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Abstract 

High-precision measurements of violations of fundamental symmetries in 
atoms are a very effective means of testing the standard model of elementary 
particles and searching for new physics beyond it. Such studies complement 
measurements at high energies. We review the recent progress in atomic parity 
nonconservation and atomic electric dipole moments (time reversal symmetry 
violation), with a particular focus on the atomic theory required to interpret 
the measurements. 
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I. INTRODUCTION 



The Glashow-Weinberg-Salam SU(3) x SU(2) x U(l) standard model of elementary 
particles [1] has enjoyed 30 years of undisputed success. It has been tested in physical 
processes covering a range in momentum transfer exceeding ten orders of magnitude. It 
correctly predicted the existence of new particles such as the neutral Z boson. However, the 
standard model fails to provide a deep explanation for the physics that it describes. For 
example, why are there three generations of fermions? What determines their masses and 
the masses of gauge bosons? What is the origin of CP violation? The Higgs boson (which 
gives masses to the particles in the standard model) has not yet been found. The standard 
model is unable to explain Big Bang baryogenesis which is believed to arise as a consequence 
of CP violation. 

It is widely believed that the standard model is a low-energy manifestation of a more 
complete theory (perhaps one that unifies the four forces). Many well-motivated extensions 
to the standard model have been proposed, such as super symmetric, technicolor, and left- 
right symmetric models, and these give predictions for physical phenomena that differ from 
those of the standard model. 

Some searches for new physics beyond the standard model are performed at high-energy 
and medium-energy particle colliders where new processes or particles would be seen directly. 
However, a very sensitive probe can be carried out at low energies through precision studies 
of quantities that can be described by the standard model. The new physics is manifested 
indirectly through a deviation of the measured values from the standard model predictions. 
The atomic physics tests that are the subject of this review lie in this second category. These 
tests exploit the fact that low-energy phenomena are especially sensitive to new physics that 
is manifested in the violations of fundamental symmetries, in particular P (parity) and T 
(time-reversal), that occur in the weak interaction. The deviations from the standard model, 
or the effects themselves, may be very small. To this end, exquisitely precise measurements 
and calculations are required. 

More than twenty years ago atomic experiments played an important role in the verifi- 
cation of the standard model. While the first evidence for neutral weak currents (existence 
of the neutral Z boson) was discovered in neutrino scattering [2], the fact that neutral cur- 
rents violate parity was first established in atomic experiments [3] and only later observed in 
high-energy electron scattering [4]. Now atomic physics plays a major role in the search for 
possible physics beyond the standard model. Precision atomic and high-energy experiments 
have different sensitivities to models of new physics and so they provide complementary 
tests. In fact the energies probed in atomic measurements exceed those currently accessible 
at high-energy facilities. For example, the most precise measurement of parity nonconser- 
vation (PNC) in the cesium atom sets a lower bound on an extra Z boson popular in many 
extensions of the standard model that is tighter than the bound set directly at the Teva- 
tron (see Section V). Also, the null measurements of electric dipole moments (EDMs) in 
atoms (an EDM is a P- and T-violating quantity) place severe restrictions on new sources 
of CP- violation which arise naturally in models beyond the standard model such as super- 
symmetry. (Assuming CPT invariance, CP-violation is accompanied by T-violation.) Such 
limits on new physics have not been set by the detection of CP-violation in the neutral K 
[5] and B [6] mesons (see, e.g., Ref. [7] for a review of CP violation in these systems). 
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Let us note that while new physics would bring a relatively small correction to a very 
small signal in atomic parity violation, in atomic EDMs the standard model value is sup- 
pressed and is many orders of magnitude below the value expected from new theories. There- 
fore, detection of an EDM would be unambiguous evidence of new physics. 

This review is motivated by the great progress that has been made recently in both 
the measurements and calculations of violations of fundamental symmetries in atoms. This 
includes the discovery of the nuclear anapole moment (an electromagnetic multipole that 
violates parity) [8], the measurement of the parity violating electron- nucleon interaction in 
cesium to 0.35% accuracy [8], the improvement in the accuracy (to 0.5%) of the atomic 
theory required to interpret the cesium measurement [9] , and greatly improved limits on the 
atomic [10] and electron [11] electric dipole moments. 

The aim of this review is to describe the theory of parity and time-reversal violation in 
atoms and explain how atomic experiments are used to test the standard model of elementary 
particles and search for new physics beyond it. We track the recent progress in the field. 
In particular, we clarify the situation in atomic parity violation in cesium: it is now firmly 
established that the cesium measurement [8] is in excellent agreement with the standard 
model; see Section V. 

The structure of the review is the following. Broadly, it is divided into two parts. The 
first part, Section II to Section VII, is devoted to parity violation in atoms. The second 
part, Section VIII to Section X, is concerned with atomic electric dipole moments. 

In Section II the sources of parity violation, and the standard model predictions, are 
described. In Section III a summary of the measurements of parity violation in atoms 
is given, with particular emphasis on the measurements with cesium. Also the atomic 
calculations are summarized. In Section IV we present a detailed description of the methods 
for high-precision atomic structure calculations applicable to atoms with a single valence 
electron. The methods are applied to parity violation in cesium in Section V and the value 
for the weak nuclear charge is extracted and compared with the standard model prediction. 
A discussion of the new physics constraints is also presented. In Section VI a brief description 
for the method of atomic structure calculations for atoms with more than one valence electron 
is given, and the thallium PNC work is discussed. A brief discussion of the prospects for 
measuring PNC along a chain of isotopes is also presented. Then in Section VII work on 
the anapole moment is reviewed. 

A description of electric dipole moments in atoms is given in Section VIII, with a sum- 
mary of all the measurements and a discussion of the P, T- violating sources at different 
energy scales. Then in Section IX a review of P, T- violating nuclear moments is given. In 
Section X a summary of the best limits on P, T- violating parameters can be found. 

Concluding remarks are presented in Section XI. 

For a general introduction to atomic P- violation and P, T- violation we refer the reader 
to the excellent books by Khriplovich [12] and Khriplovich and Lamoreaux [13]. 

II. MANIFESTATIONS AND SOURCES OF PARITY VIOLATION IN ATOMS 

Parity nonconservation (PNC) in atoms arises largely due to the exchange of Z°-bosons 
between atomic electrons and the nucleus. The weak electron-nucleus interaction violating 
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parity, but conserving time-reversal, is given by the following product of axial vector (A) 
and vector (V) currents: 



h = "7= E \C 1N e^ 5 eN^N + C 2N e^eN^ l5 N] . (1) 



G 

N 

Here G = 1.027 x 10 _5 / m p is the Fermi weak constant, AT is a nucleon wave function, and 
the sum runs over all protons p and neutrons n in the nucleus. The Dirac matrices are 
defined as 

7o= (o-°/)' 7i =(-°^o)' 75= (-°/"o J )' (2) 

and a = 2s are the Pauli spin matrices. The coefficients C\n and C 2 n give different weights 
to the contributions of protons and neutrons to the parity violating interaction. To lowest 
order in the electroweak interaction, 

C lp = l/2(l -4 sin 2 6 W ) ^0.04 , C ln = -1/2 , 

C 2p = -C 2n = l/2(l -4sin 2 9 w )g A ~ 0.05 , (3) 

where qa ~ 1.26. The Weinberg angle 6w is a free parameter; experimentally it is sin 2 Qw ~ 
0.23. The suppression of the coefficients C\ v and C 2 n due to the small factor (1 — 4sin 2 ^) 
makes \C ln \ about 10 times larger than C\ p and |C 2 at|- 

There is a contribution to atomic parity violation arising due to Z° exchange between 
electrons. However, this effect is negligibly small for heavy atoms [14-16]. It is suppressed 
by a factor (1 — 4 sin 2 9)K/(QwR(Z)) compared to the dominant electron-nucleon parity 
violating interaction, where K is a numerical factor that decreases with Z and R(Z) is a 
relativistic factor that increases with Z [15]. For 133 Cs 6S* — 7S, K pa 2 and R(Z) = 2.8 
and so the suppression factor is ~ 0.04% of the dominant amplitude [15]. This number was 
confirmed in [16]. We will consider this interaction no further. 



1. The nuclear spin-independent electron-nucleon interaction; the nuclear weak charge 

Approximating the nucleons as non-relativistic, the time-like component of the interac- 
tion (A e , Vn) is given by the nuclear spin-independent Hamiltonian (see, e.g., [12]) 

G 

h w = -^l^[ZC lp p p {r) + NC ln p n (r)\ , (4) 

Z and A^ are the number of protons and neutrons. This is an effective single-electron 
operator. The proton and neutron densities are normalized to unity, / p ntP d 3 r = 1. Assuming 
that these densities coincide, p p = p n = p, this interaction reduces to 

G 

hw = T^Qwp(rh 5 , (5) 

where Qw is the nuclear weak charge. The nuclear weak charge Q w is very close to the 
neutron number. To lowest order in the electroweak interaction, it is 



7 



Q w = -N + Z(l - 4 sin 2 6 W ) « -iV 



(6) 



This value for Q w is modified by radiative corrections. The prediction of the standard 
electroweak model for the value of the nuclear weak charge Qw in cesium is [17] 

Q^(lf Cs) = -73.10 ± 0.03 . (7) 

The nuclear weak charge Qw is protected from strong-interaction effects by conservation 
of the nuclear vector current. The clean extraction of the weak couplings of the quarks 
from atomic measurements makes this a powerful method of testing the standard model and 
searching for new physics beyond it. 

The nuclear spin independent effects arising from the nuclear weak charge give the largest 
contribution to parity violation in heavy atoms compared to other mechanisms. However, 
note that the weak interaction (5) does not always "work". This interaction can only mix 
states with the same electron angular momentum (it is a scalar). Nuclear spin-dependent 
mechanisms (see below), which produce much smaller effects in atoms, can change electron 
angular momentum and so can contribute exclusively to certain transitions in atoms and 
dominate parity violation in molecules. 



2. Nuclear spin- dependent contributions to atomic parity violation; the nuclear anapole moment 

Using the non-relativistic approximation for the nucleons, the nuclear spin-dependent 
interaction due to neutral weak currents is (see, e.g., [12]) 

G 

h-NC = C 2N a ■ cr N p(r) , (8) 

V2 N 

where ccj = 7o7i- This term arises from the space-like component of the (V e , An) coupling. 
Averaging this interaction over the nuclear state with angular momentum / in the single- 
particle approximation gives 

f j G K — 1/2 t / \ ,„> 

NC = ~7f^ 2 /(/ + i) a ' ^ ' ^ ' 

where K — (I + 1/2) (— 1) /+1 / 2 ^' and k 2 = — C 2 . There are two reasons for the suppression 
of this contribution to parity violating effects in atoms. First, unlike the spin-independent 
effects [Eq. (5)], the nucleons do not contribute coherently; in the nuclear shell model only 
the unpaired nucleon which carries nuclear spin / makes a contribution. Second, the factor 
C 2 oc (1 — Asm 2 9w) is small in the standard model. 

There is another contribution to nuclear spin-dependent PNC in atoms arising from 
neutral currents: the "usual" weak interaction due to the nuclear weak charge, hw, perturbed 
by the hyperfine interaction [18]. In the single-particle approximation this interaction can 
be written as [18,19] 



G a I 

with 



i t isr a • 1 . . . 

h Q = ^ Q —p(r) , (10) 
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«g = ~Qw^- = 2.5 x 10- 4 A 2 / 3 ^ , (11) 
o Trip rift 

Rn = To^L 1 ^ 3 is the nuclear radius, r = 1.2 fm, and //jy is the magnetic moment of the 
nucleus in nuclear magnetons. For 133 Cs, /i^ = 2.58 and kq = 0.017. 

However, the neutral currents are not the dominant source of parity violating spin- 
dependent effects in heavy atoms. It is the nuclear anapole moment K a that gives the largest 
effects [20]. This moment arises due to parity violation inside the nucleus, and manifests 
itself in atoms through the usual electromagnetic interaction with atomic electrons. The 
Hamiltonian describing the interaction between the nuclear anapole moment and an electron 



is 1 



k = ^_*_ a ■ Mr) ■ (12) 

The anapole moment n a increases with atomic number, n a oc A 2 / 3 . This is the reason it leads 
to larger parity violating effects in heavy atoms compared to other nuclear spin-dependent 
mechanisms. In heavy atoms n a ~ a A 2 ! 3 ~ 0.1 — 1 [20,21]. (Note that the interaction 
(10,11) also increases as A 2 ^ 3 , however the numerical coefficient is very small.) 

The spin-dependent contributions [Eqs. (9,10,12)] have the same form and produce the 
same effects in atoms. We will continue our discussion of the nuclear anapole moment and 
of nuclear spin-dependent effects in atoms in Section VII. 



A. Simple calculation of the weak interaction in atoms induced by the nuclear weak 

charge; the Z 3 enhancement 

In 1974 the Bouchiats showed that parity violating effects in atoms increase with the 
nuclear charge Z faster than Z 3 [14,22]. This result was the incentive for studies of parity 
violation in heavy atoms. 

Let us briefly point out where the factor of Z 3 originates. Taking the non-relativistic 
limit of the electron wave functions and considering the approximation of infinite boson 
exchange mass, the Hamiltonian (5) reduces to 

h w = 7j^( CT ' P^M + 53 ( r )°- -p)Qw , (13) 

where to, cr, p are the electron mass, spin, and momentum. The weak Hamiltonian hw mixes 
electron states of opposite parity and the same angular momentum (it is a scalar). It is a 
local operator, so we need only consider the mixing of s and p 1 / 2 states. The matrix element 
(pi/2\hw\s), with non-relativistic single-particle s and pi/2 electron states, is proportional 
to Z 2 Qw- One factor of Z here comes from the probability for the valence electron to be 
at the nucleus, and the other from the operator p which, near the nucleus (unscreened by 



In fact, the distribution of the anapole magnetic vector potential is different from the nuclear 
density. However, the corrections produced by this difference are small; see Section VII. 
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atomic electrons), is proportional to Z. The nuclear weak charge \Qw\ ~ N ~ Z. (See 
[14,22,12] for more details.) It should be remembered that relativistic effects are important, 



since Dirac wave functions diverge at r = 0, ipj oc r 7_1 , 7 = y(j + 1/2) 2 — Z 2 a 2 . Taking 
into account the relativistic nature of the wave functions brings in a relativistic factor R(Z) 
which increases with the nuclear charge Z. The factor R 10 when Z = 80. 
As a consequence, the parity nonconserving effects in atoms increase as 



that is, faster than Z 3 . 

III. MEASUREMENTS AND CALCULATIONS OF PARITY VIOLATION IN 

ATOMS 

An account of the dramatic story of the search for parity violation in atoms can be 
found in the book [12]. Below we will briefly discuss how parity violation in atoms is mani- 
fested, which experiments have yielded non-zero signals of parity violation, what quantity is 
measured in the atomic experiments, and what is required to interpret the measurements. 

Parity violation in atoms produces a spin helix, and this helix interacts differently with 
right- and left-polarized light (see, e.g., Ref. [12]). The polarization plane of linearly polar- 
ized light will therefore be rotated in passing through an atomic vapour. 

The weak interaction mixes states of opposite parity (parity violation), e.g., \p) + (3\s). 
Therefore, an Ml transition in atoms will have a component originating from an El transition 
between states of the same nominal parity, E PNC , e.g. pi/ 2 — Pz/2- The rotation angle per 
absorption length in such a transition is proportional to the ratio lm(EpNc)/Ml. While it 
may appear that it is more rewarding to study Ml transitions that are highly forbidden, 
where there is a larger rotation angle, the ordinary Ml transitions are in fact more convenient 
for experimental investigation since the angle per unit length m ImEp^c Ml (see, e.g., [12]). 

In measurements of parity violation in highly-forbidden Ml transitions, an electric field 
e is applied to open up the forbidden transition. The Ml transition then contains a Stark- 
induced El component Pstark which the parity violating amplitude interferes with. In such 
experiments the ratio lm(EpNc)/P is measured, where (3 is the vector transition polarizabil- 
ity ^stark ~ Pe. 

Atomic many-body theory is required to calculate the parity nonconserving El transition 
amplitude E PNC . This is expressed in terms of the fundamental P-odd parameters like 
the nuclear weak charge Qw- Interpretation of the measurements in terms of the P-odd 
parameters also requires a determination of Ml or j3. 



Zel'dovich was the first to propose optical rotation experiments in atoms [23]. Unfortu- 
nately, he only considered hydrogen where PNC effects are small. Optical rotation experi- 
ments in Tl, Pb, and Bi were proposed by Khriplovich [24], Sandars [25], and Sorede and 
Fortson [26]. These proposals followed those by the Bouchiats to measure PNC in highly 
forbidden transitions in Cs and Tl [22,14]. 




{pi/2\h w \s) oc R(Z)Z 2 Q W , 



(14) 



A. Summary of measurements 
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The first signal of parity violation in atoms was seen in 1978 at Novosibirsk in an optical 
rotation experiment with bismuth [3]. Now atomic PNC has been measured in bismuth, lead, 
thallium, and cesium. PNC effects were measured by optical rotation in the following atoms 
and transitions: in 209 Bi in the transition 6s 2 Qp 3 4 S 3 / 2 — Qs 2 6p 3 2 D 5 / 2 by the Novosibirsk [3], 
Moscow [27], and Oxford [28,29] groups and in the transition 6s 2 6p 3 A S 3 / 2 - 6s 2 6p 3 2 D 3/2 
by the Seattle [30] and Oxford [31,32] groups; in 6s 2 6p 2 3 P - 6s 2 6p 2 3 Pi in 208 Pb at Seattle 
[33,34] and Oxford [35]; and in the transition 6s 2 6p 2 P±/ 2 — Qs 2 6p 2 P 3 / 2 in natural Tl (70.5% 
205 T1 and 29.5% 203 T1) at Oxford [36,37] and Seattle [38]. The highest accuracy that has 
been reached in each case is: 9% for 209 Bi 4 S 3/2 - 2 D 5/2 [29], 2% for 209 Bi 4 S 3/2 - 2 D 3/2 [32], 
1% for 208 Pb [34], and 1% for Tl [38] . 

The Stark-PNC interference method was used to measure PNC in the highly-forbidden 
Ml transitions: 6s 2 S 1/2 - 7s 2 S 1/2 in 133 Cs at Paris [39-42] and Boulder [43,44,8] and 
Qs 2 Qp 2 P 1/2 - Qs 2 7p 2 P 1/2 in 203 > 205 T1 at Berkeley [45,46]. In the most precise Tl Stark-PNC 
experiment [46] an accuracy of 20% was reached. In 1997, PNC in Cs was measured with 
an accuracy of 0.35% [8] - an accuracy unprecedented in measurements of PNC in atoms. 

Results of atomic PNC measurements accurate to sub-5% are listed in Table I. 

Several PNC experiments in rare-earth atoms have been prompted by the possibility of 
enhancement of the PNC effects due to the presence of anomously close levels of opposite 
parity [47]. Another attractive feature of rare earth atoms is their abundance of stable 
isotopes. Taking ratios of measurements of PNC in different isotopes of the same element 
removes from the interpretation the dependence on atomic theory [47]; see Section VI. Null 
measurements of PNC have been reported for Ml transitions in the ground state configu- 
ration 4/ 6 6s 2 of samarium at Oxford [48,49] and for the 4/ 10 6s 2 J = 8 - 4/ 10 5d6s J = 10 
transition in dysprosium at Berkeley [50]. The upper limits were smaller than expected by 
theory. 

For a recent review of measurements of atomic PNC, we refer the reader to [51]; for a 
review of the early measurements, see, e.g., [52]. For comprehensive reviews, please see the 
book [12] and the more recent review [53]. 

B. Summary of calculations 

The interpretation of the PNC measurements is limited by atomic structure calculations. 
The theoretical uncertainty for thallium is at the level of 2.5-3% for the transition 6P\/ 2 — 
QP3/2 [54,55], and is worse for the transition QPi/ 2 — 7P\/ 2 at 6% [54] and for lead (8%) [56] 
and bismuth (12% for the 876 nm transition i S 3 / 2 — 2 D 3 / 2 and about 70% for the 648 nm 
transition 4 S 3 / 2 — 2 D^j 2 ) [56,57]. The sizeable error in the calculation for the Bi 648 nm 
transition arises because there is a strong cancellation of the zeroth order contribution by the 
first-order correlation corrections, with the amplitude then being comprised largely of the 
contributions of higher-order correlations [56,57]. Cesium is the simplest atom of interest 
in PNC experiments, it has one electron above compact, closed shells. The precision of 
the atomic calculations for Cs is 0.5% [9] (see also calculations accurate to better than 1%, 
[57,16,58]). For references to earlier calculations for the above atoms and transitions, see, 
e.g., the book [12]. 

In Table II we present the values of the most precise calculations for the PNC amplitudes 
corresponding to those atoms and transitions in which high-precision measurements (< 5% 
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error) have been performed (Table I). 



C. Cesium 

Because of the extraordinary precision that has been achieved in measurements of cesium, 
and the clean interpretation of the measurements (compared to other heavy atoms), in this 
review we concentrate mainly on parity violation in cesium. The high precision of the 
nuclear weak charge extracted from cesium has made this system important in low-energy 
tests of the standard model and has made it one of the most sensitive probes of new physics. 
Measurements of parity violation in cesium have also opened up a new window from which 
parity violation within the nucleus (the nuclear anapole moment; see Section VII) can be 
studied. 

Below we list the measurements and calculations for cesium that have been performed 
over the years, culminating in a 0.35% measurement and 0.5% calculation. 

1. Measurements 

Measurements of parity violation in the highly forbidden 6S — 7S transition in Cs were 
first suggested and considered in detail in the landmark works of the Bouchiats [14,22]. Mea- 
surements have been performed independently by the Paris group [39-42] and the Boulder 
group [43,44,8]. The results of the Cs PNC experiments are summarized in Table III. 

The Paris result in the first row is the average [41] of their (revised) results for the 
measurements of PNC in the transitions 65^=4 — 7Sp=A [39] and 6Sp =3 — 7Sp=A [40]. (The 
nuclear angular momentum of 133 Cs 1 = 7/2 and the electron angular momentum J = 1/2, 
so the total angular momentum of the atom is F — 3, 4). The Paris group have very recently 
performed a new measurement of PNC in Cs (last row) using a novel approach, chiral optical 
gain [42]. 

Each of the Boulder results [43,44,8] cited in the table is an average of PNC in the 
hyperfine transitions 65^=4 — 7Sp=3 and 6Sp =3 — 7Sp=A- The accuracy of the latest result 
is 0.35%, several times more precise than the best measurements of parity violation in other 
atoms. 

The PNC nuclear spin-independent component, arising from the nuclear weak charge, 
makes the same contribution to all hyperfine transitions. So averaging the PNC amplitudes 
over the hyperfine transitions gives the contribution from the nuclear weak charge. 

PNC in atoms dependent on the nuclear spin was detected for the first (and only) time 
in Ref. [8] where it appeared as a difference in the PNC amplitude in different hyperfine 
transitions. The dominant mechanism for nuclear spin dependent effects in atoms, the 
nuclear anapole moment, is the subject of Section VII. 

2. Calculations 

Numerous calculations of the Cs QS — 7S E PNC amplitude have been performed over the 
years. These calculations are summarized in Table IV. The many-body calculations [57,16], 
accurate to 1%, performed more than ten years ago represented a significant step forward 
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for atomic many-body theory and parity violation in atoms. At the time, these calculations 
were unmatched by the PNC measurements which were accurate to 2%. The method of 
calculation used in Ref. [57] is the subject of Sections IV, V. The method used in Ref. [16] 
is based on the popular coupled-cluster method, and we refer the interested reader to this 
work for details. 

In the last ten years a series of new measurements have been performed for quantities 
used to test the accuracy of the atomic calculations [57,16], such as electric dipole transition 
amplitudes (see [77]). The new measurements are in agreement with the calculations, re- 
solving a previous discrepancy between theory and experiment. This inspired Bennett and 
Wieman [77] to claim that the atomic theory is accurate to 0.4% rather than 1% claimed 
by theorists. Since then, a number of previously unaccounted for contributions to the PNC 
amplitude have been discovered, the Breit interaction and more recently the strong-field 
radiative corrections, that enter above the 0.4% level, but below 1% (see Section V). 

A re-calculation of the work [57], with some further improvements, was performed re- 
cently, with a full analysis of the accuracy of the PNC amplitude. This work, Ref. [9], 
represents the most accurate (0.5%) calculation to date. It is described in detail in Section 
V. 

The result of [9] differs from [57,16] by only ~ 0.1% if Breit, vacuum polarization, 
and neutron distribution corrections are excluded. One may interpret this as grounds for 
asserting that the many-body calculations [57,16,58,9] have an accuracy of 0.5% in agreement 
with the conclusion of [77]. 

IV. METHOD FOR HIGH-PRECISION ATOMIC STRUCTURE 
CALCULATIONS IN HEAVY ALKALI-METAL ATOMS 

In this section we describe methods that can be used to obtain high accuracy in calcula- 
tions involving many-electron atoms with a single valence electron. These are the methods 
that have been used to obtain the most precise calculation of parity nonconservation in Cs. 
They were originally developed in works [65,78,79,54] and applied to the calculation of PNC 
in Cs in Ref. [57]. In [57] it was claimed that the atomic theory is accurate to 1%. A com- 
plete re-calculation of PNC in Cs using this method, with a new analysis of the accuracy, 
indicates that the error is as small as 0.5% [9]. (We refer the reader to Section V, where 
this question of accuracy is discussed in general; please also see Section V for an in-depth 
discussion of PNC in Cs.) 

In this section the method is applied to energies, electric dipole (El) transition am- 
plitudes, and hyperfine structure (hfs). A comparison of the calculated and experimental 
values gives an indication of the quality of the many-body wave functions. Note that the 
above quantities are sensitive to the wave functions at different distances from the nucleus. 
Hyperfine structure, energies, and El amplitudes are dominated by the contribution of the 
wave functions at small, intermediate, and large distances from the nucleus. We concentrate 
on calculations for Cs relevant to the 65* — 7S PNC El amplitude (see Eq. (56) and Section 
VA5). 

A brief overview of the method is presented in Section IV A. For those not interested in 
the technical details of the atomic structure calculations, Sections IV B-IV H may be omitted 
without loss of consistency. 
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A. Overview 

The calculations begin in the relativistic Hartree-Fock (RHF) approximation. The N — 1 
self-consistent RHF orbitals of the core are found (N is the total number of electrons in the 
atom), and the external electron is solved in the potential of the core electrons (the V 1 *' 1 
potential). RHF wave functions, energies, and Green's functions are obtained in this way. 

Correlation corrections to the external electron orbitals are included in second (low- 
est) order in the residual interaction (V cxajCt — V 1 ^^ 1 ), where V cxact is the exact Coulomb 
interaction between the atomic electrons. The correlations are included into the external 
electron orbitals by adding the correlation potential (the self-energy operator) to the RHF 
potential when solving for the external electron. Using the Feynman diagram technique, 
important higher-order diagrams are included into the self-energy in all orders: screening of 
the electron-electron interaction and the hole-particle interaction. The self-energy is then 
iterated using the correlation potential method. 

Interactions of the atomic electrons with external fields are calculated using the time- 
dependent Hartree-Fock (TDHF) method; this method is equivalent to the random-phase 
approximation (RPA) with exchange. Using this approach we can take into account the 
polarization of the atomic core by external fields to all orders. Then the major correla- 
tion corrections are included as corrections to electron orbitals (Brueckner orbitals). Small 
correlation corrections (structural radiation, normalization) are taken into account using 
many-body perturbation theory. 

B. Zeroth-order approximation: relativistic Hartree-Fock method 
The full Hamiltonian we wish to solve is the many-electron Dirac equation 2 

N 2 

H = Eh • Pi + (P - l)m - Ze 2 /r t ] + £ _ . (15) 

i=l i<j \ T i T j\ 

Here p is the electron momentum, ot and (3 are Dirac matrices, Ze is the nuclear charge 
and N is the number of electrons in the atom (N = 55 for cesium). This equation cannot 
be solved exactly, so some approximation scheme must be used. This is done by excluding 
the complicated Coulomb term and adding instead some averaged potential in which the 
electrons move. The Coulomb term, minus the averaged potential, can be added back into 
the equation perturbatively. 

It is well known that choosing the electrons to move in the self-consistent Hartree-Fock 
potential V N ~ 1 , in the zeroth order approximation, simplifies the calculations of higher-order 
terms (we will come to this in the next section). The single-particle relativistic Hartree-Fock 
(RHF) Hamiltonian is 

h = a ■ p + {(3 - l)m - Ze 2 /r + V N ~ l , (16) 



2 In Section VA2 we discuss the inclusion of the Breit interaction into the Hamiltonian. 
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Ho = J2i ho\ where the Hartree-Fock potential 

V N 1 = Vdir + Kxch , 

is the sum of the direct and nonlocal exchange potentials created by the (N— 1) core electrons 
n, 

„=1 J I 1 " r l| 

W(r) = -e 2 £ / ^^^ ^(r) • (18) 

The direct and exchange Hartree-Fock potentials are presented diagrammatically in Fig. 1. 
The Schrodinger equation 

h il>i = e^i , (19) 

where ipi, are single-particle wave functions and energies, is solved self-consistently for the 
N—l core electrons. The Hartree-Fock potential is then kept "frozen" and the RHF equation 
(16,19) is solved for the states of the external electron. The Hamiltonian H thus generates 
a complete orthogonal set of single-particle orbitals for the core and valence electrons [80]. 

Because we are performing calculations for heavy atoms, and we are interested in inter- 
actions that take place in the vicinity of the nucleus (the weak and hyperfme interactions), 
the finite size of the nucleus needs to be taken into account. We use the standard formula 
for the charge distribution in the nucleus 

p( r ) = £ (20) 

^ v ; l + exp[(r-c)/a]' V ; 

where po is the normalization constant found from the condition / p(r)d 3 r = 1, t = a(41n3) 
is the skin-thickness, and c is the half-density radius. We take t = 2.5 fm and c = 5.6710 fm 
(( r 2)V2 = 4. 8 o4 f m ) [81]. 

Energy levels of cesium states relevant to the 6S-7S El PNC transition are presented in 
Table V. It is seen that the RHF energies agree with experiment to 10%. 

In order to obtain more realistic wave functions, we need to take into account the effect 
of correlations between the external electron and the core. We describe the techniques used 
to calculate these correlations in the following sections. 



C. Correlation corrections and many-body perturbation theory 

The subject of this section is the inclusion of electron-electron correlations into the 
single-particle electron orbitals using many-body perturbation theory. We will see that high 
accuracy can be reached in the calculations by using the Feynman diagram technique as a 
means of including dominating classes of diagrams in all orders. 

The correlation corrections can be most accurately calculated in the case of alkali-metal 
atoms (for example, cesium). This is because the external electron has very little overlap 
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with the electrons of the tightly bound core, enabling the use of perturbation theory in the 
calculation of the residual interaction of the external electron with the core. 

The exact Hamiltonian of an atom [Eq. (15)] can be divided into two parts: the first part 
is the sum of the single-particle Hamiltonians, and the second part represents the residual 
Coulomb interaction 

N 

H = J2ho( r i) + U , (21) 
i=i 

u = Y.T-^- r Y,v N - 1 ^)- (22) 

i<j I 1 * l n i=i 

Correlation corrections to the single-particle orbitals are included perturbatively in the resid- 
ual interaction U . By calculating the wave functions in the Hartree-Fock potential V N ~ l for 
the zeroth-order approximation, the perturbation corrections are simplified. The first-order 
corrections (in the residual Coulomb interaction U) to the ionization energy vanish, since 
the diagrams in first-order in the Coulomb interaction are nothing but the Hartree-Fock ones 
(Fig. 1). So two terms in Eq. (22) cancel each other. The lowest-order corrections therefore 
correspond to those arising in second-order perturbation theory, U^ 2 \ These corrections 
are determined by the four Goldstone diagrams in Fig. 2 [80]. They can be calculated by 
direct summation over intermediate states [80] or by the "correlation potential" method 
[65] . This latter method gives higher accuracy and, along with the Feynman diagram tech- 
nique to be discussed in the following section, enables the inclusion of higher-order effects: 
electron-electron screening, the hole-particle interaction, and the nonlinear contributions of 
the correlation potential. 

The correlation potential method corresponds to adding a nonlocal correlation potential 
£ to the potential V N ~ l in the RHF equation (16) and then solving for the states of the 
external electron. The correlation potential is defined such that its average value coincides 
with the correlation correction to the energy, 

5e a = (a|E|o;) (23) 
U a = I S( ri ,r 2 ,e a )^(ri)rf 3 ri . (24) 

It is easy to write the correlation potential explicitly. For example, a part of the operator 
E(ri,r 2 ,e a ) corresponding to Fig. 2(a) is given by 

A B/ v 4 [ [js j3 ^n( r 4)^24V/3(r4)^ 7 (r2)^(r3)^(r 1 )rr 3 V TO (r3) 
S a (ri,r 2 ,e Q ) = e ^ / / d r z d A r A _ — 1 . (25) 

Note that £ is a single-electron and energy-dependent operator. By solving the RHF equa- 
tion for the states of the external electron in the field V^ -1 + £, we obtain "Brueckner" 
orbitals and energies. 3 The largest correlation corrections are included in the Brueckner 
orbitals. 



3 Note that there is a slight distinction in the definition of these Brueckner orbitals and those 
defined in, e.g., [83]. 
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See Table V for Brueckner energies of the lower states of cesium calculated in the secon- 
order correlation potential. S*- 2 - 1 = denotes the "pure" second-order correlation po- 
tential (without screening, etc.). It is seen that the inclusion of these corrections improves 
the energies significantly, from the level of 10% deviation from experiment for the RHF 
approximation to the level of 1%. 

D. All-orders summation of dominating diagrams 

We saw in the previous section that when we take into account second-order correlation 
corrections, the accuracy for energies is improved significantly beyond that for energies 
calculated in the RHF approximation. However, the corrections are overestimated. This 
overestimation is largely due to the neglect of screening in the electron-electron interaction. 

In this section we describe the calculations of three series of higher-order diagrams: 
screening of the electron-electron interaction and the hole-particle interaction, which are 
inserted into the correlation potential E; and iterations of E. With the inclusion of these 
diagrams the accuracy for energies is improved to the level of 0.1% (see Table V). 

The screening of the electron-electron interaction is a collective phenomenon and is sim- 
ilar to Debye screening in a plasma; the corresponding chain of diagrams is enhanced by a 
factor approximately equal to the number of electrons in the external closed subshell (the bp 
electrons in cesium) [78]. The importance of this effect can be understood by looking at a 
not dissimilar example in which screening effects are important, for instance, the screening 
of an external electric field in an atom. According to the Schiff theorem [84], a homogeneous 
electric field is screened by atomic electrons (and at the nucleus it is zero). (See [85] where 
a numerical calculation of an external electric field inside the atom has been performed.) 

The hole-particle interaction is enhanced by the large zero-multipolarity diagonal matrix 
elements of the Coulomb interaction [79]. The importance of this effect can be seen by 
noticing that the existence of the discrete spectrum excitations in noble gas atoms are due 
only to this interaction (see, e.g., [86]). 

The non-linear effects of the correlation potential are calculated by iterating the self- 
energy operator. These effects are enhanced by the small denominator, which is the energy 
for the excitation of an external electron (in comparison with the excitation energy of a core 
electron) [79]. 

All other diagrams of perturbation theory are proportional to powers of the small pa- 
rameter Q n d/Aei llt ~ 10~ 2 , where Q n d is a nondiagonal Coulomb integral and Ae int is a large 
energy denominator corresponding to the excitation of a core electron [79]. 

1. Screening of the electron- electron interaction 

The main correction to the correlation potential comes from the inclusion of the screening 
of the Coulomb field by the core electrons. Some examples of the lowest-order screening 
corrections are presented in Fig. 3. When screening diagrams in the lowest (third) order 
of perturbation theory are taken into account, a correction is obtained of opposite sign and 
almost the same absolute value as the corresponding second-order diagram [78] . Due to these 
strong cancellations there is a need to sum the whole chain of screening diagrams. However, 
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this task causes difficulties in standard perturbation theory as the screening diagrams in 
the correlation correction cannot be represented by a simple geometric progression due to 
the overlap of the energy denominators of different loops (such an overlap indicates a large 
number of excited electrons in the intermediate states, see e.g. Fig. 3(b,c)). This summation 
problem is solved by using the Feynman diagram technique. 

The correlation corrections to the energy in the Feynman diagram technique are pre- 
sented in Fig. 4. The Feynman Green's function is of the form 

where ip n is an occupied core electron state, ip 7 is a state outside the core. While the 
simplest way of calculating the Green's function is by direct summation over the discrete 
and continuous spectrum, there is another method in which higher numerical accuracy can 
be achieved. As is known, the radial Green's function Go for the equation without the 
nonlocal exchange interaction V exc h can be expressed in terms of the solutions xo and Xoo 
of the Schrodinger or Dirac equation that are regular at r — > and r — > oo, respectively: 
Co( r i> 7-2) ~ Xo( r <)Xoo( r >)j r < = min(r 1 ,r 2 ), r> = max(ri,r 2 ). The exchange interaction is 
taken into account by solving the matrix equation G = G + G V exc hG. The polarization 
operator (Fig. 5) is given by 

^- GriT2 (u; + e)G r2ri (e). (27) 
This integration is carried out analytically, giving 

n rir2 H = E u i{ ll I'll , 2 ^(ri)V 7 (ri)^(r 2 )Vn(r 2 ) 

71,7 v n 7' 

= »E^( r i)fe + w ) + 6 ( e " - ^)]^n(r 2 ). (28) 

n 

Using formulae (26) and (28), it is easy to perform analytical integration over u in the 
calculation of the diagrams in Fig. 4. After integration, diagram 4(a) transforms to 2(a,c) 
and diagram 4(b) transforms to 2(b,d). 

Electron-electron screening, to all orders in the Coulomb interaction, corresponds to 
the diagram chain presented in Fig. 6. To calculate this we perform summation of the 
polarization operators before carrying out the integration over uo. The whole sum of screening 
diagrams in Fig. 6 can be represented by 

ii(uj) = + iQTliuj)}- 1 . (29) 

The integration over u is performed numerically. The integration contour is rotated 90° 
from the real axis to the complex uj plane parallel to the imaginary axis (see Fig. 7) - this 
aids the numerical convergence by keeping the poles far from the integration contour. 

The all-order electron-electron screening reduces the second-order correlation corrections 
to the energies of S and P states of 133 Cs by 40%. 
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2. The hole-particle interaction 



The hole-particle interaction is presented diagrammatically in Fig. 8. This diagram 
describes the alteration of the core potential due to the excitation of the electron from the 
core to the virtual intermediate state. This electron now moves in the potential created 
by the N — 2 electrons, and no longer contributes to the Hartree-Fock potential. Denoting 
Vo as the zero multipolarity direct potential of the outgoing electron, the potential which 
describes the excited and core states simultaneously is [79] 

V = V N - 1 - (1 -P)Vo(l -P) , (30) 
where P is the projection operator on the core orbitals, 

N-l 

P=H \n){n\ . (31) 

n=l 

The projection operator P is introduced into the potential to make the excited states orthog- 
onal to the core states. It is easily seen that for the occupied orbitals (V) = (l/^ 1 ), while 
for the excited orbitals (V) = (V N - 1 ) - (V ). Strictly one should also make subtractions for 
higher multipolarities and for the exchange interaction as well, however these contributions 
are relatively small and are therefore safe to ignore [79]. 

To obtain high accuracy, the hole-particle interaction in the polarization operator needs 
to be taken into account in all orders (see Fig. 9). This is achieved by calculating the Green's 
function in the potential (30) and then using it in the expression for the polarization operator 
(28). The screened polarization operator, with hole-particle interaction included, is found 
by using the Green's function in Eq. (29). 

The Coulomb interaction, with screening and the hole-particle interaction included in all 
orders, is calculated from the matrix equation [79] 

Q = Q- iQtxQ. (32) 

This is depicted diagrammatically in Fig. 10. 

The infinite series of diagrams representing the screening and hole-particle interaction can 
now be included into the correlation potential. This is done by introducing the renormalized 
Coulomb interaction (Fig. 10) and the polarization operator (Fig. 9) into the second-order 
diagrams according to Fig. 11. 

The screened second-order correlation corrections to the energies of S and P states of 
cesium are increased by 30% when the hole-particle interaction is taken into account in all 
orders. 

3. Chaining of the self-energy 

The accuracy of the calculations can be further improved by taking into account the 
nonlinear contributions of the correlation potential £ (see Fig. 12). The chaining of the 
correlation potential (Fig. 11) to all orders is calculated by adding £ to the Hartree-Fock 
potential, V N ~ 1 ) and solving the equation 
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(ho + t- e)i) = 



(33) 



iteratively for the states of the external electron. The inclusion of £ into the Schrodinger 
equation is what we call the "correlation potential method" and the resulting orbitals and 
energies "Brueckner" orbitals and "Brueckner" energies (see Section IV C). 

Iterations of the correlation potential £ increase the contributions of S (with screening 
and hole-particle interaction) to the energies of S and P states of cesium by about 10%. 

The final results for the energies are listed in Table V. The inclusion of the three series 
of higher-order diagrams improves the accuracy of the calculations of the energies to the 
level of 0.1%. 

E. Other low-order correlation diagrams 

Third-order diagrams for the interaction of a hole and particle in the polarization loop 
with an external electron are depicted in Fig. 13. These are not taken into account in 
the method described above. However, these diagrams are of opposite sign and cancel each 
other almost exactly [79]: the small and almost constant potential of a distant external 
electron practically does not influence the wave functions of the core and excited electrons 
in the loop; it shifts the energies of the core and excited electrons by the same amount. This 
cancellation was proved in the work [87] by direct calculation. 

Also, correlation corrections to the external electron energy arising from the inclusion 
of the self-energy into orbitals belonging to closed electron shells, depicted in Fig. 14, are 
small and can be safely omitted [65]. 

F. Empirical fitting of the energies 

The calculations of the external electron wave functions can be refined by placing co- 
efficients before the self-energy operator, £ — > /£, such that the energies are reproduced 
exactly. This can be considered as a way of including higher-order diagrams not explicitly 
included in the calculations. Comparison of quantities calculated with and without fitting 
can be used to test the stability of the wave functions and to estimate the contribution of 
unaccounted diagrams. 

G. Asymptotic form of the correlation potential 

At large distances, the correlation potential £ approaches the local polarization potential 

[88], 

E r _oo « -ae 2 /2r 4 , (34) 

where a is the polarizability of the core. This explains the universal behaviour of the 
correlation corrections to the energies of states of the external electron. 
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H. Interaction with external fields 



In this section we will describe the procedures used to achieve high accuracy in the 
calculations of the interactions between atomic electrons and external fields (in particular, we 
will present calculations for El transition amplitudes and hyperfine structure (hfs) constants, 
arising from the interaction of the atomic electrons with the electric field of the photon and 
with the magnetic field of the nucleus, respectively). 

We need to calculate the effect of the external field on the wave functions of the core 
electrons (core polarization) and then take into account the effect of this polarization and 
the external field on the valence electron. This is achieved by using the time- dependent 
Hartree-Fock (TDHF) method. We describe this method in Section IV H 1 and apply it to 
the calculation of El transition amplitudes and hfs constants in Sections IV H 2 and IV H 3, 
respectively. 

The dominant correlation corrections correspond to those diagrams in which the interac- 
tions occur in the external lines of the self-energy operator ( "Brueckner-type corrections"). 
These diagrams are presented in Fig. 15; they are enhanced by the small energy denominator 
e cxt corresponding to the excitation of the external electron in the intermediate states. The 
Brueckner-type corrections are calculated in a similar way to the correlation corrections to 
energy (see Sections IVH1, IVH2, IVH3). In Section IVH4 we describe the calculations 
of the remaining second-order corrections, i.e., "structural radiation" and normalization of 
states. Structural radiation diagrams are presented in Fig. 16; in these diagrams the exter- 
nal fields occur in the internal lines, and so their contributions are small due to the large 
energy denominators ei nt corresponding to the excitation of the core electrons. The rela- 
tive suppression of the contributions of structural radiation compared to Brueckner-type 
corrections is e ext /e int ~ 1/10. 

This section is largely based on the works [54,88]. 



1. Time- dependent Hartree-Fock method 

The polarization effects are taken into account by using the time-dependent Hartree-Fock 
method (see, e.g., [54,88] and references therein). We will consider how the RHF equations 
are modified in the presence of a time-dependent field 

h ext = (fe-^ + pe^). (35) 

We can assume that the time-dependent single-particle orbitals are then given by 

^k = ^k + Xke~ iut + ri k e iut ), (36) 
with corresponding eigenvalues ( "quasi-energies" ) 

~e k = €k + Se^e^ + 5e k 2) e lu)t ; (37) 

Xk-, Vk and b~e k l \ 8e^ are corrections to the RHF wave functions ipk and energies e k , respec- 
tively, induced by h ext . The Schrodinger equation 
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* _ . _ d ~ 

{h (ip) + h ext )^k = ifo^k, (38) 
can now be used to obtain equations for the corrections Xk and i] k , 

{h - e k - u)xk = -(/ + 8V ext )ip k + Se k ip k 

(h - e k + uj)r] k = -(/ f + 5V^ t )ip k + 5e k tp k , (39) 

where we take into account terms up to first order; the energy shift 5e k = 5e k ^ = 5e k 2 ^ = 
(i^klf + 8Vext\ipk)- The RHF Hamiltonian h(ip) corresponds to ho with the Hartree-Fock 
potential V N ~ l calculated with the new core wave functions ip, and SV cxt is the difference 
between the potential found in the external field and the RHF potential, 

N-l , i3„ 

5V cxt ^(r) = e 2 J2 '-Mi^Mn) + ^(r 1 ) Xn (r 1 )Wr) - (40) 

n =i J |r-ri|LV ) 

-(^(ri)^n(r) +^(r 1 )x„(r))^(r 1 )' . 

Eqs. (39), (40) should be solved self-consistently for the (N — 1) core electrons. The wave 
function of the external electron is then found in the field of the frozen core (the V^ -1 
approximation). So in the same way as in the RHF case, here we can find a complete set of 
ortho normal TDHF orbitals ip k with quasi-energy e k . 

The external electron transition amplitude Mp a from state \a) to state \f3) induced by 
the field h cxt can be found by comparison of the amplitude obtained from Eqs. (36) and 
(39), 

(MM = Wx«e-'> = <Mf+i%!M e -« (41) 

e« - tf3 + LU 

with conventional time-dependent perturbation theory, 

4> a = + Mp \ #e~^, (42) 

where only the resonant term (uj — e a ) is considered. Comparing Eqs. (42) and (41) 
gives [54,88] 

Mfr = (if>p\f + 6V ext \1> a ) . (43) 

This formula corresponds to the well-known random-phase approximation (RPA) with ex- 
change (see, e.g., [86]. When the orbitals a and (3 are calculated in the potential V N ~ 1 , the 
transition amplitude corresponds to the RHF value. 

Using the TDHF procedure described above, core polarization is included in all orders 
of perturbation theory. This is equilavent to summation of the diagram series presented in 
Fig. 17. 

The Brueckner-type correlation corrections are calculated in the same way as the cor- 
rections to energies. Brueckner, instead of RHF, orbitals are used for a and (5 in Eq. (43). 



22 



This is equivalent to calculating the diagrams presented in Fig. 15. However, it must be 
noted that by using this technique we neglect the diagrams in which the interaction takes 
place in the internal lines (see Fig. 16), although these diagrams give only small corrections; 
we will deal with these contributions in Section IV H 4. Also in Section IV H 4 we look at 
another second-order correction: the normalization of the many-body states. 

2. El transition amplitudes 

The Hamiltonian of the electron interaction with the electric field 

E(t) = E {e- iuJt + e iujt ) (44) 

of an electromagnetic wave depends on the choice of gauge. In "length" form fi — er • E 
and in "velocity" form /„ = — ie(a ■ E /u;), where h E i = fe~ tuJt + Pe wt . 

It is known that in TDHF calculations the amplitude (43) is gauge invariant (see, e.g., 
[86,89,67]. Comparing results obtained from the two forms for the amplitudes (43) provides 
a test of the numerical calculation. In [88], the length and velocity forms were shown to 
give the same results for the El transition amplitudes when the correlation corrections 
(correlation potential, frequency shift in the velocity form operator, structural radiation, 
and normalization of states) are taken into account. Calculations using the dipole operator 
in length form are more stable than those obtained in velocity form (see, e.g., [90,91,88]). 
For this reason we consider the calculations of El transition amplitudes in length form. 

As seen from Eq. (43), inclusion of the core polarization into the El transition amplitude 
is reduced to the addition of the operator 8Ve\- Because the El operator can only mix 
opposite parity states, there is no energy shift, and so when calculating the wave function 
corrections Xk and rjk from Eq. (39), Set = 0. 

The results for El transition amplitudes calculated in length form between the lower 
states of cesium are presented in Table VI. 

Another test of the accuracy and self-consistency of the TDHF equations is the value of 
an external electric field at the nucleus. According to the Schiff theorem [84], an external 
static electric field (u = in TDHF) at the nucleus in a neutral atom is shielded completely 
by atomic electrons. For an atom with charge Zi the total static electric field E tot at the 
nucleus is [85] 

E tot (0) = E + (E e (0)> = E Z t /Z , (45) 

where E is the external field and E e is the induced electron field. TDHF equations reproduce 
this result correctly. The oscillations of the electric field inside the atom are quite complex. 
We refer the reader to [85] for a plot of the electric field inside Tl + . 4 



4 Note that in the paper [85] the figures and the captions have been switched. 
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3. Hyperfine structure constants 



The hyperfine interaction between a relativistic electron and a point nucleus is given by 

h h f = ea-A, (46) 

where 

A = ^ (47) 

is the vector potential created by the nuclear magnetic moment fi and a. is the Dirac matrix. 
When performing the calculation of the hyperfine structure of heavy atoms, it is important 
to take into account the finite size of the nucleus. Using a simple model in which the nucleus 
represents a uniformly magnetized ball, the hyperfine interaction is given by 




h hf = en ■ f (r) , f (r) = _ < J ' ' ; ' ™ (48) 



where we take the magnetic radius r m = l.lA 1//3 fm (A is the mass number of the nucleus). 
Note that while the distribution of currents in the nucleus (produced by unpaired nucleons) 
is very complex, the hfs is only weakly dependent on its form (see, e.g., [99]). 

Corrections to the RHF wave functions induced by the hyperfine interaction are calcu- 
lated using Eq. (39). There is no time- variation in the hyperfine interaction, and so we set 
uj — 0. The hfs equations are then 

(h - e)<ty = -(h hf + 8V hf )il> + <W (49) 
6e={if>\h hf + 6V hf \if>) . (50) 

The hyperfine interaction does not alter the direct contribution to the core potential, so 
fiVhf — Vexch(ip) — V C xch(ip)- We use the V N ~ l approximation to calculate the complete set 
of orbitals. The external electron correction Se determines the atomic hyperfine structure. 

Hyperfine structure results for cesium are presented in Table VII. The correlations in- 
crease the density of the external electron in the nuclear vicinity by about 30%. An accurate 
inclusion of the correlations is therefore very important when considering interactions sin- 
gular on the nucleus, for example the PNC weak interaction. 



4- Structural radiation and normalization of states 

"Structural radiation" is the term we use for correlation corrections with the external 
field in the internal lines. (We will make a further distinction when we are dealing with the 
PNC El amplitude Ep^c, since here there are two fields. We call correlation corrections 
with the weak interaction in the internal lines "weak correlation potential"; see Section 
VA.) Examples of diagrams which represent structural radiation are presented in Fig. 16. 
In second-order perturbation theory in the residual interaction there is another contribution 
which arises due to the change of the normalization of the single-particle wave functions due 
to admixture with many-particle states. 
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For the El transition amplitudes in length form, the following approximate formula is 
used to calculate the structural radiation [88], 

1 f)Y f)Y 

M str = --(/3|D^ + ^D|a), (51) 

where D = d + 5V. The derivation can be found in Ref. [88]. We refer the interested 
reader to Refs. [88,54] for the structural radiation in the velocity gauge. The normalization 
contribution for El transitions has the form [88,54] 

M norm = -((3\T>\ a )(((3\^\(3) + (a\^\a)) . (52) 

The structural radiation and normalization contributions to the hyperfine structure and 
El transition amplitudes of low-lying states of cesium are small. Their combined contribu- 
tion, for both hfs and El amplitudes, usually lies in the range 0.1 — 1.0%. 



V. HIGH-PRECISION CALCULATION OF PARITY VIOLATION IN CESIUM 
AND EXTRACTION OF THE NUCLEAR WEAK CHARGE 

In this section the method described in the preceding section is applied to the QS — 7S 
parity violating amplitude in cesium and the value for the nuclear weak charge is extracted 
from the measurement of Wood et al. [8]. The value for this amplitude has been the source 
of much confusion recently. It has jumped around erratically in the last few years, and 
finally it has stabilized. The origin of this instability is described below. 

In 1999, experimentalists Bennett and Wieman argued that the standard model is in 
contradiction with atomic experiments by 2.5 a [77]. This conclusion was mainly based 
on their analysis of the accuracy of atomic structure calculations which were published in 
Ref. [57] in 1989 and in Ref. [16] in 1990,1992. The point is that the new measurements 
of electromagnetic amplitudes in atoms have demonstrated that the accuracy of the atomic 
calculations is much better than it seemed to be ten years ago. Indeed, all disagreements 
between theory and experiments were resolved in favour of theory. Based on this, Bennett 
and Wieman reduced the theoretical error from the 1% claimed by theorists [57,16] to 0.4% 
and came to the conclusion that there may be new physics beyond the standard model. 

Particle physicists put forward new physics possibilities, such as an extra Z boson in the 
weak interaction, leptoquarks, and composite fermions; see, e.g., Refs. [103]. 

However, the 1% error placed on the atomic calculations [57,16] was based not only on 
the comparison of measured and calculated quantities such as El transition amplitudes. 
It was based also on an allowance for unnaccounted contributions. It was soon after dis- 
covered that the Breit contribution is larger than expected (-0.6%) [104] (Section VA2). 
Then it was suggested that strong-field radiative corrections could make a contribution as 
large as Breit [105]. The Uehling contribution has been found to give a contribution of 
0.4% [106,107] (Section VA4). And just recently it has been established that self-energy 
and vertex contributions are ~ —0.8% [108,109] (Section VA4). A recent comprehensive 
calculation of atomic PNC is accurate to 0.5% [9]. The results of this calculation, and the 
discussion of accuracy, are presented in the following sections. The final value for parity 
violation in cesium is in agreement with the standard model and tightly constrains possible 
new physics (Section VC). 
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A. High-precision calculations of parity violation in cesium 



The weak interaction H w = Y^f h l w [Eq. (5)] 5 mixes atomic wave functions of opposite 
parity and leads to small opposite-parity admixtures in atomic states ^, ^ = ^ + 8^. This 
gives rise to El transitions between states of the same nominal parity. The parity violating 
QS — 7S El transition amplitude in Cs is 

E PNC = {fS\H E1 \QS) = (5(7S)\H E1 \6S) + (7S\H E1 \5(6S)) . (53) 

Calculations of PNC El amplitudes can be performed using the following approaches: 
from a mixed-states approach, in which there is a small opposite-parity admixture in each 
state [Eq. (53)]; or from a sum-over-states approach, in which the amplitude [Eq. (53)] 
is broken down into contributions arising from opposite-parity admixtures and a direct 
summation over the intermediate states is performed [Eq. (54)]. 

In the sum-over-states approach, the Cs QS — 7S PNC El transition amplitude is written 
in terms of a sum over intermediate, many-particle states NPi/ 2 

^ l (7S\H E1 \NP 1/2 )(NP 1/2 \H w \QS) (7S\H W \NP 1/2 )(NP 1/2 \H E1 \QS) ] 
Epnc = 2^ p p + p v • ^ 

If one neglects configuration mixing, this sum can be represented in terms of single-particle 
states; in this case, the sum also runs over core states (corresponding to many-particle states 
with a single core excitation), 

^ r(7s\h E1 \np 1/2 )(np 1/2 \h w \6s) (7s\h w \np 1/2 ) (np 1/2 \h E1 \6s)i 
&pnc — 2^ [ 1 J • 

There are three dominating contributions to this sum: 

(7s\h E1 \6p 1/2 )(Qp 1/2 \h w \Qs) (7s|/i w |6p 1/2 )(6p 1/2 |/i B i|6s) 

&PNC = 1 

| (7s\h E i\7p 1/2 ) (7p l/2 \h w \Qs) + 

e 6s — e 7p 1/2 

= -1.908 + 1.493 + 1.352 + ... = 0.937 + ... . 

The numbers are from the work [16] where the sum-over-states method was used; here we 
just demonstrate that these terms dominate. An advantage of the sum-over-states approach 
is that experimental values for the energies and El transition amplitudes can be explicitly 
included into the sum. This was the procedure for some of the early calculations of PNC in 
Cs (see, e.g., [71]). 

In Refs. [57,16,58,9] PNC calculations were performed in the mixed-states approach, and 
in Ref. [16] a calculation was carried out in the sum-over-states approach also. Here we refer 
to the most precise calculations, < 1% accuracy. 



5 It is seen from Eqs. (3,4), by inserting the coefficients Cuv, that the density p(r) is essentially 
the (poorly understood) neutron density in the nucleus. In the calculations, p(r) will be taken 
equal to the charge density, Eq. (20), and then in Section VA3 we will consider the effect on the 
PNC El amplitude as a result of correcting for p(r). 
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1. Mixed-states calculation 



In the TDHF method (Section IV HI, Eq. (36)), a single-electron wave function in 
external weak and El fields is 

i/j = ?fj + S^fj + Xe~ iuJt + Ye itjjt + 8Xe~ iujt + 5Ye iu}t , (57) 

where ipo is the unperturbed state, 5ip is the correction due to the weak interaction acting 
alone, X and Y are corrections due to the photon field acting alone, and 5X and SY are 
corrections due to both fields acting simultaneously. These corrections are found by solving 
self-consistently the system of the TDHF equations for the core states 



(ho - e)Stp = -(h w + SV w )ip, (58) 

(ho-t- u)X = -(h E1 + 5V E1 )ij, (59) 

(h - e + u)Y = -(h E1 + 5V E1 )^, (60) 

(h -e- uj)5X = -5V E1 5ij - 5V W X - 5V Elw i>, (61) 

(h -e + u)5Y = -5V E1 S^ - 5V W Y - 5V Elw ij, (62) 



where 5Vw and 5V E \ are corrections to the core potential due to the weak and El interactions, 
respectively, and 8V E \w is the correction to the core potential due to the simultaneous action 
of the weak field and the electric field of the photon. 

The TDHF contribution to E PNC between the states 6S and 7S is given by 

E T pnc F = (Mh E1 + 5V E1 \5^ 6s ) + (^ 7s \h w + 5V w \X 6s ) + (^ 7s \5V Elw \^ s ) . (63) 

The corrections 5ipQ S and X Qs are found by solving the equations (58-59) in the field of the 
frozen core (of course, the amplitude (63) can instead be expressed in terms of corrections 
to tp 7s ). 

Now we need to include the correlation corrections to the PNC El amplitude. In the 
previous sections (Sections IVC,IVD,IVH4) we have discussed two types of corrections: 
the dominant Brueckner-type corrections, represented by diagrams in which the external 
field appears in the external electron line (see Fig. 18); and structural radiation, in which 
the external field acts on an internal electron line. In the case of PNC El amplitudes, in 
order to distinguish between structural radiation diagrams with different fields, we refer to 
diagrams with the weak interaction attached to the internal electron line as "weak correlation 
potential" diagrams. Structural radiation and the weak correlation potential diagrams are 
presented in Fig. 19. 

We will consider first the dominating Brueckner-type corrections to the El PNC ampli- 
tude. Remember that the correlation potential is energy- dependent, S = S(e). This means 
that the £ operators for the 6s and 7s states are different. We should consider the proper 
energy-dependence at least in first-order in E (higher-order corrections are small and the 
proper energy-dependence is not important for them). The first-order in S correction to 
E PNC is presented diagrammatically in Fig. 18. We can write this as 

{tlj 7s \± s (e 7s )\5X 6s ) + (5^p 7s \± p (e 7s )\X 6s ) + (5Y 7s \± s (e 6s )\^p 6s ) + (Y 7s \± p (e 6s )\5^ 6s ) . (64) 

The non-linear in S contribution to the Brueckner-type correction is found using the 
correlation potential method (Section IV C): the all-orders in £ contribution is calculated 
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and from this the first-order contribution, found in the same method, is subtracted. The 
all-orders term is calculated using external electron orbitals, and corrections to these orbitals 
induced by the weak interaction and the photon field, found in the potential V N ^ + E. The 
PNC El amplitude is then calculated, using these new orbitals, in the same way as in the 
usual time-dependent Hartree-Fock method. The all-orders contribution to E PNC is 

^all-orders = EpNC ^N-l + £ } _ EpNC {V N - 1 ) . (65) 

The first-order in E contribution is found by placing a small coefficient a before the corre- 
lation potential, E — > aE. When a < 1, the linear in E contribution to E P nc dominates. 
Its extrapolation to a = 1 gives the first-order in E contribution. So the non-linear in E 
contribution to E PNC is [110] 



7-inon— lin 
^PNC 



E PNC {V N - 1 + E) - E PNC {V N - l )\ - -\E PNC {V N -' + Q E) - E PNC {V N - X ) 



(66) 

To complete the calculation of corrections second-order in the residual Coulomb inter- 
action the weak correlation potential, structural radiation, and normalization contributions 
to the PNC amplitude must be included. 

The weak correlation potential is calculated by direct summation over intermediate 
states. See Section IV H 4 for the approximate form for structural radiation in length form, 
and for the form for the normalization of the many-body states. Due to parity violation 
there is an opposite-parity correction to the orbitals a and (3, a = a + 5a and j3 = (3 + 5/3, 
and to the correlation potential E, E = E + 5E. 

Structural radiation is then given by 

i _ fly fly 

M str = -i(/?|D^ + ^D|a>. (67) 

There are two contributions to structural radiation for the PNC El amplitude: one in which 
the electromagnetic vertex is parity conserving, the weak interaction included in the external 
lines: 

i _ fly fly 

M Fig . 19(b) = --(/3|D— + — D|a>. (68) 

(see diagram (b) Fig. 19); and the other in which the weak interaction is included in the 
electromagnetic vertex (we call this structural radiation and not weak correlation potential): 



M Fig . 19(c) = -2<0|D— + ^ D l«>- ( 69 ) 



(see diagram (c) of Fig. 19). Note that in each case the amplitude first-order in the weak 
interaction is considered. 

The normalization contribution is 

Mnorm = ^ 0\ D| &) (((3\ ^ \(3) + («%!«>)■ (70) 
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The results of the calculation [9] for the 6S — 7S PNC amplitude are presented in Ta- 
ble VIII. Taking into account all corrections discussed in this section, the following value is 
obtained for the QS — 7S PNC amplitude in cesium 

E PNC = 0.9078 x 10- n iea B (-Q w /N) . (71) 

This corresponds to "Subtotal" of Table VIII. This is in agreement with the 1989 result 
[57]. Notice the stability of the PNC amplitude. The time- dependent Hartree-Fock value 
gives a contribution to the total amplitude of about 98%. The point is that there is a strong 
cancellation of the correlation corrections. 

The mixed-states approach has also been performed in [16] and [58] to determine the 
PNC amplitude in cesium. However, in these works the screening of the electron-electron 
interaction was included in a simplified way. In [16] empirical screening factors were placed 
before the second-order correlation corrections to fit the experimental values of energies. 
Kozlov et al. [58] introduced screening factors based on average screening factors calculated 
for the Coulomb integrals between valence electron states. The results obtained by these 
groups (without the Breit interaction, i.e., corresponding to the Subtotal of Table VIII) are 
0.904 [16] and 0.905 [58]. 6 As a check, a pure second-order (i.e., using E( 2 )) calculation with 
energy-fitting was also performed in [9] (in the same way as [16]), and the result 0.904 was 
reproduced. 

Contributions of the Breit interaction, the neutron distribution, and radiative corrections 
to E PN c are considered in the following sections. 



2. Inclusion of the Breit interaction 
The Breit interaction is a two-particle operator 

tt e2 V- gj - a J + (gj • n ij)( a j ■ n ij) {70 ^ 
-"Breit = -77 2^ i i > K' Z ) 

a are Dirac matrices, n = (r, — rj)/|r? — Tj\. It gives magnetic (Gaunt) and retardation 
corrections to the Coulomb interaction. A few years ago it was thought that the correction 
to Epnc arising due to inclusion of the Breit interaction in the Hamiltonian (15) is small 
(safely smaller than 1%). In the work [57] the Breit interaction was neglected, and in [16] it 
was only partially calculated. (Remember that these works claimed an accuracy of 1%.) The 
huge improvement in the experimental precision of the cesium PNC measurement in 1997 
[8] and the claim of Bennett and Wieman in 1999 [77] that the theoretical accuracy is 0.4% 
prompted theorists to revisit their calculations. Naturally this also involves a consideration 
of previously neglected contributions which, while at the 1% level could be neglected, are 



6 The numbers differ from those presented in Table II due to the Breit interaction. In [16] a value 
for Breit of —0.2% of the PNC amplitude was included (this value was underestimated), while in 
[58] the magnetic (Gaunt) part of the Breit interaction was included and calculated to be —0.4%. 
See Section V A 2. 
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significant at the 0.4% level. Derevianko [104] calculated the contribution of the Breit 
interaction to Ep^c and found that it is larger than had been expected. Its contribution to 
E PNC is -0.6%. This result has been confirmed by subsequent calculations [111,58,9]. 

3. Neutron distribution 

The weak Hamiltonian Eq. (5) was used to obtain the result Eq. (71) with p(r) taken 
to be the charge density, parametrized according to Eq. (20). However, as we mentioned in 
a footnote at the beginning of Section VA, the weak interaction is sensitive to the distri- 
bution of neutrons in the nucleus. Here we look at the effect of correcting for the neutron 
distribution. 

For the neutron density the two-parameter Fermi model (20) is used. The result of Ref. 
[112] was used in [9] for the difference Ar np = 0.13(4) fm in the root-mean-square radii of 
the neutrons (r 2 ) 1 / 2 and protons (r 2 ) 1 / 2 . Three cases which correspond to the same value 
of (r 2 ) were considered: (i) c n = c p , a n > a p ; (ii) c n > c p , a n > a p ; and (iii) c n > c p , 
o-n = a P (using the relation (r 2 ) |c 2 + ^7r 2 a^) . It is found that Ep^c shifts from —0.18% 
to —0.21% when moving from the extreme c n = c p to the extreme a n = a p . Therefore, Ep^c 
changes by about —0.2% (—0.0018) due to consideration of the neutron distribution. This 
is in agreement with Derevianko's estimate, —0.19(8)% [113]. 

4- Strong-field QED radiative corrections 

It was noted in Ref. [105] that corrections to the PNC amplitude due to vacuum po- 
larization by the strong Coulomb field of the nucleus could be comparable in size to the 
Breit correction. This has been confirmed by calculations, the strong-field radiative correc- 
tions associated with the Uehling potential (vacuum polarization) increase Ep^c by 0.4% 
[107,106,9]. 

In Ref. [114] it was pointed out that the self-energy correction can give a larger contribu- 
tion to 133 Cs PNC with opposite sign (~ —0.65%). The self-energy and vertex corrections 
were first calculated in [108] and found to be —0.73(20)% for 133 Cs. The relation between 
the PNC correction and radiative corrections to finite nuclear size energy shifts was used in 
this work. This result was confirmed in direct analytical calculations using Za expansion 
[115,109,116] and by all-orders in Za numerical calculations of the PNC matrix element of 
the 2S\/2 — 2Pi/2 transition in hydrogenic ions performed in [117]. 

Note that corrections occur at very small distances (r < 1/m) where the nuclear Coulomb 
field is not screened and the electron energy is negligible. Therefore, the relative radiative 
corrections to weak matrix elements in neutral atoms like Cs are approximately the same as 
for the 2<Si/2 — 2Py 2 transition in hydrogenic ions. 

There is good agreement between the different calculations for all values of Z; see [118] 
and the review [119]. For the strong-field self-energy and vertex contribution to PNC in 
133 Cs we will quote the value —0.8%, which is the average value of [108,116] corresponding 
also to the value obtained in [117]. 

Above we discussed the radiative corrections to the weak matrix elements. However, the 
sum-over-states expression for the PNC amplitude contains also energy denominators and 
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El electromagnetic amplitudes. It was shown in [114,9] that for Cs the corrections to the 
energies (—0.3%) and El matrix elements (+0.3%) cancel. 

The contributions of strong-field radiative corrections to E PNC of cesium are listed in 
Table VIII. 

5. Tests of accuracy 

There are two main methods used to estimate the accuracy of the PNC amplitude E PNC - 
(i) root-mean-square (rms) deviation of the calculated energy intervals, El amplitudes, and 
hyperfine structure constants from the accurate experimental values; (ii) influence of fitting 
of energies and hyperfine structure constants on the PNC amplitude. 

The PNC amplitude can be expressed as a sum over intermediate states (see beginning 
of Section VA). Notice that there are three dominating contributions to the 6S — 7S PNC 
amplitude in Cs; see Eq. (56). Each term in the sum is a product of El transition amplitudes, 
weak matrix elements, and energy denominators. Therefore, this amplitude is sensitive to 
the electron wave functions at all distances. (The weak matrix elements, energies, and El 
amplitudes are sensitive to the wave functions at small, intermediate, and large distances 
from the nucleus, respectively.) While mixed-states calculations of PNC amplitudes do not 
involve a direct summation over intermediate states, it is instructive to analyze the accuracy 
of the weak matrix elements, energy intervals, and El transition amplitudes which contribute 
to Eq. (56) calculated using the same method as that used to calculate Ep^c- The accuracy 
of these quantities is determined by comparing the calculated values with experiment. Note 
that we cannot directly compare weak matrix elements with experiment. However, like the 
weak matrix elements, hyperfine structure is determined by the electron wave functions in 
the vicinity of the nucleus, and this is known very accurately. 

In Section IV we presented calculations of the energies, El transition amplitudes, and 
hyperfine structure constants relevant to Cs QS — 7S E PNC - The states that have been 
considered in these calculations are those relavant to E PNC [in the sum (56)]. 

The calculated removal energies are presented in Table V. The Hartree-Fock values 
deviate from experiment by 10%. Including the second-order correlation corrections 
reduces the error to ~ 1%. When screening and the hole-particle interaction are included 
into S( 2 ) in all orders, the energies improve, 0.2 — 0.3%. The rms deviation between the 
calculated and experimental energy intervals e 6s — e 6pi/2 , e 7s — £q Pi/21 and e 6s — e 7pi/2 is 0.3%. 

We mentioned at the end of Section IV D that the experimental values for energies can 
be fitted exactly by placing a coefficient before the correlation potential £. The stability of 
the amplitude Ep^c (as well as the El amplitudes and hfs constants) with fitting gives us 
an indication of the size of omitted contributions. Note that the accuracy for the energies 
is already very high and the remaining discrepancy with experiment is of the same order of 
magnitude as the Breit and radiative corrections. Therefore, generally speaking, we should 
not expect that fitting of the energy will always improve the results for amplitudes and 
hyperfine structure. In fact, as we will see below, some values do improve while others do 
not. The overall accuracy, however, remains at the same level. 

Below we present results for Ep^c obtained in three different approximations: with 
unfitted £, and with E^ 2 ) and £ fitted with coefficients to reproduce experimental removal 
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energies. First, we analyse the El transition amplitudes and hfs constants calculated in 
these approximations. 

The relevant El transition amplitudes (radial integrals) are presented in Table VI. These 
are calculated with the energy-fitted "bare" correlation potential E( 2 ) and the (unfitted and 
fitted) "dressed" potential E. Structural radiation and normalization contributions are also 
included. The rms deviations of the calculated El amplitudes from experiment are the 
following: without energy fitting, the rms deviation is 0.1%; fitting the energy gives a rms 
deviation of 0.2% for E( 2 ) and 0.3% for the complete E. Note, these correspond to the 
deviations between the calculations and the central points of the measurements. The errors 
associated with the measurements are in fact comparable to this difference. So it is unclear 
if the theory is limited to this precision or is in fact much better. Regardless, the uncertainty 
in the theoretical accuracy remains the same. 

The hyperfine structure constants calculated in different approximations are presented in 
Table VII. Corrections due to the Breit interaction, structural radiation, and normalization 
are included. The rms deviation of the calculated hfs values from experiment using the 
unfitted E is 0.5%. With fitting, the rms deviation in the pure second-order approximation is 
0.3%; with higher orders it is 0.4%. The point is to estimate the accuracy of the s—pi/2 weak 

matrix elements. It seems reasonable then to use the square-root formula, y / hfs(s)hfs(pi/2)- 
Notice that by using this approach the deviation is smaller. Without energy fitting, the rms 
deviation is 0.5%. With fitting, the rms deviation in the second-order calculation (E^ 2 ^) is 
0.2% and in the full calculation (E) it is 0.3%. 

From the above consideration it is seen that the rms deviation for the relevant parameters 
is 0.5% or better. Note that from this analysis the error for a sum-over-states calculation 
of E PNC would be larger than this, as the errors for the energies, hfs constants, and El 
amplitudes contribute to each of the three terms in Eq. (56). However, in the mixed-states 
approach, the errors do not add in this way. 

We now consider calculations of the PNC amplitude performed in [9] in different approx- 
imations (with unfitted E, and with energy-fitted E^ 2 ^ and E). The spread of the results can 
be used to estimate the error. The results are listed in Table IX. It can be seen that the PNC 
amplitude is very stable. The PNC amplitude is much more stable than hyperfine structure. 
This can be explained by the much smaller correlation corrections to E PNC (~ 2% for E PNC 
and ~ 30% for hfs; compare Table VIII with Table VII). One can say that this small value 
of the correlation correction is a result of cancellation of different terms in (64) but each 
term is not small (see Table VIII). However, this cancellation has a regular behaviour. The 
stability of E PNC may be compared to the stability of the usual electromagnetic amplitudes 
where the error is very small (even without fitting). 8 



7 Note that fitting is an empirical method to estimate screening corrections (which were 
accurately calculated in £). Agreement between results with fitted X^ 2 ** and ab initio S shows that 
the fitting procedure is a reasonable way to estimate omitted diagrams. 

8 Note that different methods also give different signs of the errors for hfs. This is one more 
argument that the true value of Ep^c is somewhere in the interval between the results of different 
calculations in Table IX. 
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In [9] the fitting of hyperfine structure was also considered, using different coefficients 
before each E. The first-order in E correlation correction (64) changes by about 10%. It 
was found that the PNC amplitude changes by about 0.4%. 

It is also instructive to look at the spread of E PNC obtained in different schemes. The 
result of the work [9] (the number we present here) is in excellent agreement with the earlier 
result [57] while the calculation scheme is significantly different. The only other calculation 
of the Epnc in Cs which is as complete as [57,9] is that of Blundell et al. [16]. Their result 
in the all-orders sum-over-states approach is 0.909 (without Breit) and is very close to the 
value of 0.908 (corresponding to "Subtotal" of Table VIII). 

A note on the sum-over-states procedure. The authors of reference [16] include single, 
double, and selected triple excitations into their wave functions. Note, however, that even if 
wave functions of 65', 7S, and intermediate NP states are calculated exactly (i.e., with all 
configuration mixing included) there are still some missed contributions in this approach. 
Consider, e.g., the intermediate state 6P = 5p 6 6p. It contains an admixture of states 
5p 5 nsQd: QP = 5p 6 6p + a5p 5 ns6d + .... This mixed state is included into the sum (54). 
However, the sum (54) must include all many-body states of opposite parity. This means 
that the state 5p 5 ns6d = 5p 5 ns6d — a5p 6 6p + ... should also be included into the sum. 
Such contributions to Ep^c have never been estimated directly within the sum-over-states 
approach. However, they are included into the mixed-states calculations [57,16,58,9]. 

It is important to note that the omitted higher-order many-body corrections are different 
in the sum-over-states [16] and mixed-states [57,9] calculations. This may be considered as 
an argument that the omitted many-body corrections in both calculations are small. Of 
course, here it is assumed that the omitted many-body corrections to both values (which, in 
principle, are completely different) do not "conspire" to give exactly the same magnitude. 

A comparison of calculations of Ep^c in second-order with fitting of the energies is also 
useful in determining the accuracy of the calculations of E PNC - (Remember that this value 
is in agreement with results of similar calculations performed in [16,58]; see Section VAl.) 
One can see that replacing the all-order E by its very rough second-order (with fitting) 
approximation changes Ep^c by less than 0.4% only. On the other hand, if the higher 
orders are included accurately, the difference between the two very different approaches is 
0.1% only. 

The maximum deviation obtained in the above analysis is 0.5%. This is the error claimed 
in the E PNC calculation [9]. 

B. The vector transition polarizability 

The determination of the nuclear weak charge from the Stark-PNC interference measure- 
ments also requires knowledge of the vector transition polarizability (3. This can be found 
in a number of ways: 

(i) from a direct calculation of f3. (3 can be expressed as a sum over intermediate states 
and experimental El transition amplitudes and energies can be used [14] (see also [16,120]). 
However, this calculation is unstable due to strong cancellations of different terms in the 
sum (see Ref. [120]). These cancellations are explained by the fact that (3 is proportional 
to the spin-orbit interaction, therefore for zero spin-orbit interaction the sum for (3 must be 
zero; 



33 



(ii) from the measurement of the ratio of the off-diagonal hyperfme amplitude to the 
vector transition polarizability, Mhf s /j3 [121]. f3 is then extracted from the ratio using a 
theoretical determination of M h f s ; 

(hi) from the measurement of the ratio of the scalar to vector polarizabilities, a/ (3. a can 
be calculated accurately using experimental values for El transition amplitudes and energies 
in the sum-over-states approach (the calculation of a is much more stable than that of (3 
[120]). 

There are currently two very precise determinations of (3. One was obtained from 
the analysis [122] (calculation of M h f s ) of the measurement [77] of the ratio M h f s /{3, 
(3 = 26.957(51)af 3 , and another is from an analysis [9] (semi-empirical calculation of a; 
see [120] for details, where a similar calculation was performed) of the measurement [123] 
of the ratio a/ (3 using the most accurate experimental data for El transition amplitudes 
including the recent measurements of Ref. [124], (3 = 27.15(ll)a 3 3 . An average of these 
values gives 

(3 = 26.99(5)4 • (73) 



C. The final value for the Cs nuclear weak charge Qw and implications 

Combining the measurement [8] 

- Im(£ T g) = 1.5935(56)^ 
p cm 

with the calculated value (see Table VIII) 

E PNC = 0.897(1 ± 0.5%) x lO~ n iea B (-Q w /N) (75) 

(from the calculation [9] with the averaged value —0.8% of works [108,109] for the self-energy 
and vertex radiative corrections) and the averaged value for (3 [Eq. (73)], gives 

Q w = -72.74(29) exp (36) theor (76) 

for the value of the nuclear weak charge for 133 Cs. The difference between this value and 
that predicted by the standard model, Q^^fCs) = -73.19 ± 0.13 [125], 9 is 

AQ- W = Q W - Q S w = 0.45(48) , (77) 

adding the errors in quadrature. 

Let us briefly consider the constraints on physics beyond the standard model set by the 
nuclear weak charge of cesium [Eq. (77)]. New physics corrections to the standard model 
are divided into two groups, those that originate through vacuum polarization corrections to 



(74) 



9 We use this value rather than the Particle Data Group value, Eq. (7), since we use the new 
physics analysis of Ref. [125]. 
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gauge boson propagators and those that originate from all other mechanisms (such as new 
physics arising through vertex and self-energy diagrams and new tree-level physics). The 
former are termed "oblique" and the latter "direct" . 

A formalism for oblique corrections has been devised by Peskin and Takeuchi [126] in 
terms of weak isospin conserving and weak isospin breaking parameters S and T, respectively. 
Atomic parity violation (on a single isotope) is unique among other electroweak probes of 
new physics in its almost exclusive dependence on the parameter S. Indeed, the dependence 
of the nuclear weak charge on the parameter T cancels almost exactly. For 133 Cs [125] 

AQ$ liquc = -0.8005 - 0.007T . (78) 

The standard model value corresponds to S = T = (no new physics) at values m t = 
174.3 GeV for the top quark mass and M H = 100 GeV for the Higgs boson mass. 10 The 
constraint on S from PNC in ^Cs [comparing Eqs. (77) and (78)] is 

S = -0.56(60) , (79) 

adding the errors in quadrature. 

In terms of direct new physics, a positive AQ^ W [Eq. (77)] could be indicative of an 
extra Z boson in the weak interaction. A lower bound for the Z x boson mass predicted 
in 5*0(10) theories can be obtained from the deviation of the measured weak charge from 
theory, according to [128] 

^Qw\ ee « 0.4(2iV + Z){M w /M Zx f . (80) 

To one standard deviation, the lower bound on the mass M z% from parity violation in ,j| 3 Cs 
is 

M Zx > 750 GeV . (81) 

A lower bound of about 600 GeV has been obtained from a direct search at the Tevatron 
[129]. 

For a discussion of atomic physics sensitivities to new physics, see, e.g., [130,131,125] 
and references therein. For a recent analysis of electroweak tests of the standard model, 
including atomic parity violation, see, e.g., [125,132,133]; an earlier review on this topic 
[134] is also very informative. 

D. Ongoing/future studies of PNC in atoms with a single valence electron 

A PNC measurement in cesium is continuing at Paris [42], where ~ 1% precision is 
expected to be reached. The possibility of performing PNC measurements on the single 
trapped ions Ba + and Ra + is being considered at Seattle [135,136]. Preliminary atomic 



10 The result is actually not very sensitive to Mh which is currently Mh > 114.4 GeV at the 95% 
confidence level [127]. 
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PNC calculations [137] indicate that the accuracy for calculations of Ba + could compete 
with that of cesium. An experiment to measure PNC in francium was discussed at Stony 
Brook [138,139]. Francium is a heavier analog of cesium, and correspondingly the PNC 
effect is an order of magnitude (18 times [110]) larger than that for cesium. However, there 
are no stable isotopes in francium. The group uses laser trapping to suspend the atoms. 

A new facility (TRRiP) is under development at the Kernfysisch Versneller Instituut, 
Groningen. One aim is to measure parity violation in radioactive atoms and ions, attractive 
candidates being Ra + and Fr. See Ref. [140]. 

Note that the atomic theory is limited by the error associated with calculations of corre- 
lation corrections. It's expected that the relative corrections due to correlations are the same 
in analogous atoms, e.g., cesium and rubidium. Thus, taking ratios of PNC measurements 
in analogous atoms possibly provides a way to circumvent the troublesome correlations in 
the determination of the ratio of the nuclear weak charges. [Another way to exclude atomic 
theory is through isotope ratios, see Section VI B.] Light atoms by themselves may provide 
some advantage from the point of view of atomic theory. Breit, neutron skin, and QED 
radiative corrections for rubidium are much smaller than for cesium. Therefore, atomic 
PNC calculations for rubidium could, in principle, reach a precision better than for cesium, 
however the amplitude is much smaller [70,54], about 6 times [54]. 

VI. ATOMS WITH SEVERAL ELECTRONS IN UNFILLED SHELLS 

There are some advantages in measuring PNC effects in heavier atoms (where the PNC 
signal is expected to be larger) with relatively complicated electronic structure. However, the 
interpretation of the measurements in such a case is strongly impeded by the poor knowledge 
of the atomic wave functions. Even for Tl, which has a relatively simple structure, the atomic 
theory has an error of 2.5-3% [54,55]. 

The method we discussed in Sections IV and V is applicable for heavy alkali-metal atoms. 
However, this method may not be accurate for atoms with more than one electron in unfilled 
shells. This is because the correlations between the electrons in the unfilled shells may be 
important. 

The most effective many-body method for atoms with several electrons above closed 
shells was developed in the works [141,142] (see also [56]). It is a combination of many-body 
perturbation theory and configuration interaction methods (it is called "MBPT+CI"). An 
effective Hamiltonian is constructed for the valence electrons using many-body perturbation 
theory for the interaction of the valence electrons with the core. In this way the correlations 
between the external electrons and the core are taken into account using MBPT while the 
correlations between the external electrons are calculated using the CI method. 

A detailed description of the MBPT+CI method is beyond the scope of this review. We 
refer the interested reader to the works [141,142]. 

A. Parity nonconservation in thallium 

High-precision measurements of parity violation in the 6P1/2— QP3/2 transition in thallium 
have been performed by the Oxford and Seattle groups, 
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(82) 



We will use the latter, most precise measurement to extract the value for Q w . A recent 
measurement of the ratio E2/M1 indicates that the Oxford result should be rescaled to 
1Z = —15.34 ± 0.45 [143], reducing the discrepancy between the two measurements (82). 

Let us move on now to the calculations. The electronic configuration for the ground 
state of thallium is 6s 2 6pi/ 2 . This system can be treated as a one-electron or a three- 
electron system above closed shells. However, the Qp electron is energetically close to the 6s 
electrons, and correlations between 6s and 6p are significant. This means that the method 
described in Sections IV,V will not be as effective for thallium as for cesium. This method 
was employed in the work [54] to the transition 6pi/ 2 — 6p 3 / 2 , yielding the following result: 

E PNC = (-2.70 ± 0.08) x 10- 10 iea B (-Q w /N) . (83) 

This result includes core polarization of El and weak fields calculated in the TDHF method 
and all second order correlation corrections (higher-order diagrams - electron-electron 
screening and the hole-particle interaction - are not included). Finally, empirical fitting 
of the energies is used to take into account some missed higher-order correlations. 
Consideration of the Breit interaction, Eq. (72), gives a correction [144] 

AEf^/E PNC = -0.0098 . (84) 

The size of radiative corrections to Tl PNC is the following. The Uehling correction is 0.94% 
[145] and we take the self-energy and vertex corrections to be —1.51% (average of —1.61% 
from [108] and -1.41% from [109]), 

AEgftc/EpNc = -0.0057 . (85) 
The correction for the neutron distribution is small (with a relatively large error) [55], 

AE™%™/E PNC = -0.003 . (86) 

With these corrections, Eq. (83) becomes 

E PNC = (-2.65 ± 0.08) x 10' 10 iea B (-Q w /N) . (87) 

(the errors associated with the corrections to Eq. (83) are well below 1%.) 

More recently the method MBPT+CI was applied in Ref. [55], with the result 

E PNC = (-2.72 ± 0.07) x l0~ l0 iea B (-Q w /N) . (88) 

This method is well-suited for a system like thallium, where there are three electrons above 
a compact core, with the correlations between the three external electrons treated non- 
perturbatively. In this work the Gaunt interaction (the dominant, magnetic part of the 
Breit interaction, Eq. (72)) was included self-consistently at every stage of the calculation. 
Inclusion of radiative corrections shifts the value (88) slightly, E PNC = (—2.70 ± 0.07) x 
lO~ w iea B (— Q w /N). The difference between this value and Eq. (87) is within the assigned 
errors. 
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We need a value for the 6P1/2 — QP3/2 magnetic dipole transition amplitude, Ml. Calcu- 
lations of this quantity are very stable. In Ref. [55] the values Ml = 1.692 x 1CT 3 a.u. and 
Ml = 1.694 x 1CT 3 a.u. were calculated, the former using the MBPT+CI method and the 
latter in the MBPT method (treating thallium as a single-electron atom). 

The Seattle measurement (82), the PNC calculation (87), and the Ml amplitude calcu- 
lated in [55] lead to the following value for the nuclear weak charge of thallium, 

Q^( 205 T1) = -116.3(1.3) exp (3.5) theor , (89) 

in agreement with the standard model value [17] 

Q^ M ( 205 T1) = -116.67(5) (90) 



B. A method to exclude the error from atomic theory: isotope ratios and the 

neutron distribution 

In Ref. [47] it was pointed out that atomic theory (and hence the large associated error) 
can be excluded by taking ratios of atomic PNC measurements along an isotope chain. 
Indeed, from a simple non-relativistic consideration, the measured parity violating amplitude 
can be expressed as 

Apnc = iQw , (91) 

where the atomic theory is included into £. It is seen that in this approximation taking 
ratios of measurements of different isotopes with neutron numbers N and N' 

Apnc(N') = Q W (N') 

Ap N c(N) Q W (N) { ] 

removes the dependence on £, assuming that the atomic structure does not change signifi- 
cantly from isotope to isotope. However, a consideration of relativistic effects (the variation 
of the electron wave functions inside the nucleus) changes this simple formulation. 

It was shown in Ref. [146] that while the atomic structure cancels in the isotope ratios, 
there is an enhanced sensitivity to the neutron distribution p n (r), and the uncertainty as- 
sociated with p n (r) may limit the extraction of new physics. The dependence of the PNC 
amplitude on nuclear structure effects can be incorporated through the correction Q nnc [146], 

Apnc = Z(Qw + Q n w) , (93) 

where 

Q n w = -N(q n - 1) + Z(l -Asm 2 9w )(q p - 1) (94) 

and 

q n = J Pn(r)f(r)d 3 r , q p = J p p (r)f(r)d 3 r . (95) 
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The variation of the electron wave functions inside the nucleus is given by f(r), which is 
normalized to /(0) = 1. The isotope ratio 



where Aq n = q n — q' n . It is seen that R is sensitive, in particular, to the difference in the 
neutron distributions, and it is the error associated with this difference, S(Aq n ), that could 
limit the interpretation of the ratio R of the PNC measurements in heavy atoms in terms 
of new physics. 

The isotope ratios are sensitive to a combination of new physics that is different from 
that probed by measurements of parity violation in a single isotope. While measurements 
of atomic parity violation in a single isotope are sensitive to the weak isospin conserving 
parameter S and to new tree level physics, the isotope ratios are sensitive to both S and T 
through sin 2 9w as well as to new tree level physics. The sensitivity of the isotope ratios to 
oblique new physics has been investigated in [128,147,148], and the sensitivity to direct new 
physics in [131,148]. 

Nuclear structure corrections to the weak charge distribution have been calculated for 
different isotopes of lead [148,149], cesium [150,149,151], barium [149], and ytterbium [149]. 
The calculations are model-dependent, and it is seen from the spread in the results that 
at this time our knowledge of the neutron distribution would preclude a determination 
of new physics competitive with other probes. The problem is that there is a very limited 
amount of empirical information on neutron distributions. In Ref. [152] new data on neutron 
distributions from experiments with antiprotonic atoms was used to re-analyze the impact 
of the nuclear structure uncertainties on the extraction of new physics. It was found that the 
errors associated with nuclear structure are slightly reduced compared with those associated 
with nuclear calculations. It was suggested that isotope ratios of PNC measurements with 
atoms having Z < 50 may be more effective at searching for new physics, since the nuclear 
structure uncertainty in the extraction of (direct) new physics increases roughly as Z 8 ^ 3 and 
the required experimental uncertainty is less strict for lighter atoms [152]. Since the PNC 
signal oc Z 3 , this would only be feasible if there is an accidental enhancement of the PNC 
effect (for example, very close levels of opposite parity). 

The high sensitivity of isotope ratios of PNC measurements to variations in the neutron 
distribution could be used to determine the nuclear structure and test nuclear models [146]. 
Parity violating electron scattering measurements appear to be very promising for extracting 
information on the neutron distribution (see, e.g, [153]). It appears that such measurements 
could be used to significantly reduce the nuclear structure uncertainties to the level where 
the new physics sensitivity of PNC isotope ratios is significant. 



A thallium PNC measurement is continuing at Seattle [154] in which sub-1% accuracy 
is expected to be reached. Experimental studies of rare-earth elements dysprosium and yt- 
terbium are in progress at Berkeley [155]. Studies with dysprosium are continuing despite 
the small upper limit found for the PNC amplitude [50]. Dysprosium appeared to be par- 
ticularly attractive for PNC studies due to the presence of a state (Af 9 5d 2 Qs J = 10) that 




(96) 



C. Ongoing/future studies of PNC in complex atoms 
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is nearly degenerate with the state 4f 10 5dQs J = 10 with opposite parity and the same 
angular momentum. However, the weak interaction does not mix the dominant electronic 
configurations of the nearly degenerate states, a non-zero weak matrix element shows up 
only through configuration mixing and core polarization which makes it very small. An 
improvement in the statistical sensitivity of a few orders of magnitude is anticipated in the 
current search [50,156]. If the effect is of the same order of magnitude as the current upper 
limit, then dysprosium will provide an important test of the standard model and allow the 
determination of the nuclear anapole moment [50,156]. A measurement of PNC in ytterbium 
is in progress also at Berkeley [157,158]. In ytterbium there are seven stable isotopes in the 
range A = 168 — 176 with two of these isotopes having non-zero nuclear spin (possibility 
to observe the nuclear anapole moment). The atomic theory in this case should be more 
reliable than in the other rare-earth atoms, since here there is a closed 4/-shell. 

PNC studies with samarium appear to have no advantages over other systems. Prelimi- 
nary work showed that PNC optical rotation on the forbidden Ml transition 4/ 6 6s 2 7 F — 
Af 6 5dQs 7 Gi, which could have benefited from having close (~ 10 cm -1 ) opposite parity 
states, is not feasible [159]. Moreover, it has recently been found, through a term reassign- 
ment of several of the odd-parity states in Sm [160], that an estimate of the weak matrix 
element made in the original proposal for the experiment [161] should be reduced by about 
40 times. The null PNC results of the Oxford group reported in Refs. [48,49] for Ml tran- 
sitions between the ground state configuration show that the possible enhancement due to 
the relatively small energy denominators (~ 200 cm -1 ) is not realized due to very small El 
and weak matrix elements. 

There are some interesting proposals for new experiments where large enhancements are 
expected. E.g., atomic calculations [162,163] indicate that parity violation in the electronic 
transitions 7s 2 1 S , — 7sQd 3 Di of the radium isotopes 223 Ra and 225 Ra are enhanced due 
to a very small interval between the states 7s 7p 3 P\ and 7s6d 3 D\. The parity violating 
amplitude is calculated to be 100 times larger than that of Cs. (Note that the electronic 
structure of radium is relatively simple and so it is expected that the PNC amplitude will 
not suffer from large suppressions that appeared in the PNC amplitudes of the rare earth 
atoms.) 

VII. THE NUCLEAR ANAPOLE MOMENT AND MEASUREMENTS OF P-ODD 
NUCLEAR FORCES IN ATOMIC EXPERIMENTS 

A. The anapole moment 

The notion of the anapole moment was introduced by Zel'dovich [164] just after the 
discovery of parity violation. He noted that a particle may have a parity violating elec- 
tromagnetic form factor, in addition to the usual electric and magnetic form factors. The 
first realistic example, the anapole moment of the nucleus, was considered in Ref. [20] and 
calculated in Ref. [21]. In these works it was also demonstrated that atomic and molecular 
experiments could detect anapole moments. Subsequently, a number of experiments were 
performed in Paris, Boulder, Oxford, and Seattle [40,43,44,37,38] and some limits on the 
magnitude of the anapole moment were established. However, it was only recently that a 
nuclear anapole moment was unambiguously detected - in 1997 a group in Boulder measured 
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a nuclear anapole moment in 133 Cs (using atomic experiments) to an accuracy of 14% [8] 
(see also the recent review [177]). This is the first observation of an electromagnetic moment 
that violates fundamental discrete symmetries. 

Besides the usual electric and magnetic monopole, dipole and quadrupole (and so on) 
moments, there are also other electromagnetic multipole moments, which are not usually 
dealt with in multipole moment expansions as they give rise to contact, rather than long- 
range, potentials. The anapole moment is such a moment. It obeys time reversal invariance 
but violates parity conservation and charge conjugation invariance (i.e., T-even and P- and 
C-odd). The anapole moment arises out of an expansion of the vector potential as a series in 
P _1 , where R is the distance from the center of the charge distribution (see, e.g., [165,12]). 
The part of the vector potential that is due to the anapole moment is 11 



j is the electromagnetic current density. Eq. (98) can be taken as the definition of the 
anapole moment. Notice the contact form of the potential - this is true for any T-even, 
P-odd moment; see, e.g., [20,12] for a proof. 

What kind of current distribution corresponds to an anapole moment? Such a current 
distribution is shown in Fig. 20. This will give a nonzero anapole moment, because the 
places at which the current is pointing downwards are further from the center than the 
places where the current is pointing upwards. Since the current in Eq. (98) is weighted by a 
factor of r 2 this current distribution will produce an anapole moment pointing perpendicular 
to the plane of the doughnut. A magnetic field is produced inside the current distribution, 
as shown in the figure. 

The expression for the anapole moment (98) contains the current vector j, which changes 
its sign under reflection of co-ordinates. The anapole moment is directed along the nuclear 
spin I: (a) = — vr(r 2 j) = al/I. However, the spin I does not change its sign under co- 
ordinate reflection. The different behavior of the right and left hand sides of the relation 
(r 2 j) oc I under reflection of co-ordinates means that the existence of the anapole moment 
violates parity, i.e., symmetry under the reflection of co-ordinates (but it does not violate 
time reversal invariance). 



We now turn to the question of how an anapole moment can actually be produced. 
A T-even, P-odd moment like the anapole can only arise if there is some kind of P-odd 
force present; for this the weak interaction is needed. The two-body P-odd nucleon-nucleon 
interaction, in the contact limit, has the form (see, e.g., [12]) 



In the gauge V • A = there is a long-range term in the anapole vector-potential which may 
be removed by a gauge transformation. 




(97) 



where 




(98) 



B. Origin of the nuclear anapole moment 
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G 1 

^ ab = V2 2m ({ ( 9ab(7a ~ 9ba(Tb) • (p a - Pfc) , 5(r„ - r fe ) } 

+g' ab \(T a x <r 6 ] • V5(r„ - r 6 )) , (99) 

where { , } is an anticommutator, G is the Fermi constant of the weak interaction, m is the 
nucleon mass, and cr, p, and r are the spins, momenta, and coordinates of the nucleons a and 
b. The dimensionless constants g a b, gba, g' a b §i ye the strength of the weak interaction between 
nucleons. The effective one-body P-odd weak interaction between an unpaired nucleon and 
the nuclear core can then be obtained, 

W = ^^9[<t ■ PP(r) + P(r)<r • p] , (100) 

where p{r) is the number density of core nucleons and 

Z N . , 

g = jg av + -jg an • (ioi) 

Here a = p,n denotes the unpaired nucleon. For the 133 Cs atom the unpaired nucleon is a 
proton and so g = g p . 

The interaction (100) perturbs the wave function of the unpaired nucleon, resulting in 
the mixing of opposite parity states: ip = ipo + Sip, where ipo = |0) is the unperturbed wave 
function and Sip = J2 n \ n ){ n \W\0)(E — En)' 1 . An approximate analytical solution for the 
perturbed Schrodinger equation (H + W)ip = Eip (which assumes that the nuclear density 
is constant) gives (see, e.g., [21,166]) 

^ = e i6(T - T ^ , (102) 

where 9 = — gGpj\[2. What this means is that the spin (s = |cr) of the unperturbed 
wave function will be rotated around the vector r by an angle of 29r. If, for example, 
the unperturbed wave function was in a spin up state, the spin at different points for the 
perturbed wave function will be as shown in Fig. 21. Thus we have a spin helix, with a 
definite chirality, i.e., right- or left-handedness (see, e.g., [12]). This means that the parity 
symmetry has been broken. Now consider the current and magnetic field produced by such 
a spin helix. The electromagnetic current of the unpaired nucleon in the non-relativistic 
limit has the form 

j = -|^ f VV - (V^)V] + |£v x (^W), (103) 

where q — (1) for a neutron (proton) and p is the nucleon magnetic moment in nuclear 
magnetons. The first term comes from the orbital motion of the nucleon (convection or 
orbital current), while the second term is a magnetic moment current term, which produces 
the dominating contribution. The current distribution and the magnetic field produced by 
the wave function tp of Eq. (102) have been calculated, e.g., in Ref. [167]. 

Using the expression for the electromagnetic current (103) in Eq. (98), the operator of 
the anapole moment, a (a = (^|a|^>)) can be written as 

a = (ne/m)\fi(r x cr) - (g/2)(pr 2 + r 2 p)], (104) 
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where r and p are the position and momentum operators of the nucleon. The domi- 
nant contribution to the nuclear anapole comes from the first, spin term and so we can 
express the anapole moment operator in terms of the magnetic dipole moment operator 
M = tr(e J u)/(2m) as a « 2?r(r x M). 

The anapole moment is usually described by a dimensionless parameter, K a , defined by 
the following equation 

1 G KI . . 

^eV^WTTf- (105) 

where K = (7+|)(— iy +1 / 2 - 1 and / is the orbital angular momentum of the external nucleon. 

In Ref. [21] an approximate analytical result for n a (in terms of the parameter g p ) was 
obtained by using the wave function (102) to calculate the mean value of the anapole moment 
operator (104). The result is 

Ka = ^ — AV3 9 P ( 106 ) 

10 mr 

= 0.08ft for 133 Cs , (107) 

where a = 1/137, fi is the magnetic moment of the external nucleon in nuclear magnetons, 
and tq = 1.2 fm. It is the dependence on A 2 / 3 that makes the anapole moment the dominant 
nuclear spin-dependent effect in heavy atoms. 

Single-particle calculations of K a for Cs in the Woods-Saxon and harmonic oscillator 
potentials have been performed in Refs. [21,168,19,174] (note that in Refs. [168,19] the effects 
of configuration mixing were also included semi-empirically) and many-body calculations in 
Refs. [169-172]. The results of these calculations are presented in Table X; a compilation of 
these results in terms of DDH "best values" (see Section VII E) can also be found in Ref. 
[173]. Single-particle calculations for the anapole moment are remarkably stable with the 
choice of potential [174]. We will quote the value 

n s a p = 0.06ft, (108) 

for 133 Cs, obtained from numerical calculations in the Woods-Saxon potential with spin-orbit 
interaction [21,174]. We will leave the discussion of many-body effects to Section VII E. 



C. Parity violating effects in atoms dependent on the nuclear spin 

The nuclear anapole moment interacts with an atom's electrons due to its magnetic field. 
The interaction is [using Eqs. (97) and (105)] 

h a = eat ■ A = eac ■ a5 3 (r) = -^=J^-^ K J 3 ( r ), (109) 

where A is the anapole vector-potential and a. is the relativistic velocity operator (Dirac 
matrices). The magnetic field corresponding to the anapole is localized inside the nucleus. 
Therefore the interaction with atomic electrons occurs only if the electron wave functions 
penetrate the nucleus. 
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The anapole moment produces parity violating nuclear spin-dependent effects in atoms. 
In heavy atoms it is the dominant mechanism producing such effects. However, its ef- 
fects are indistinguishable from the smaller parity violating neutral current effects, and 
these need to be accounted for in the extraction of re a from atomic measurements (we men- 
tioned these effects briefly in Section II 2). One contribution to parity violating nuclear 
spin-dependent effects arises from Z° exchange in an electron-nucleus interaction, with an 
axial- vector Z°-nucleus coupling and a vector Z°-electron coupling [Eq. (9)]; see Ref. [175]. 
Another contribution arises from the perturbation of the nuclear spin-independent contri- 
bution, corresponding to the nuclear weak charge, by the hyperfine interaction (Eq. (10)); 
see [18,19,176]. The total effective nuclear spin-dependent interaction can be expressed in 
the form 12 

h '" = T2 K WTT) a - Ip{rh (uo) 



where 



if -1/2 J + l nm 

K = K a — K 2 + — -g~ K Q i ( 1U ) 



andK = (/ + l/2)(-l) /+1 / 2 - / . 



D. Measurement of nuclear spin-dependent effects in cesium and extraction of the 

nuclear anapole moment 

The anapole moment can be detected by observing the amplitudes of transitions between 
atomic levels that violate parity. In the case of the Boulder experiment [8], it was an El 
transition between the 6S and 7S states of the cesium atom. The nuclear spin-dependent 
P-odd effects can be separated from the dominant effect (the weak interaction between the 
electron and the weak charge of the nucleus; see Sections III,V) by observing the dependence 
of the parity violating effects in two different hyperfine transitions. The nuclear spin has 
different relative orientations in the different hyperfine states. For 133 Cs, F = 3 and F — 4, 
where F is the total angular momentum of the atom (F = J + I, where J is the electron's 
angular momentum). By observing the transitions 6SV=4 — * 7Sp=3 and 65^=3 — > 7Sp = 4, 
the anapole moment can be detected. In the Boulder experiment [8] these two amplitudes 
were measured, and they were found to be significantly different, indicating the presence of 
an anapole moment. 



1. Atomic calculations and extraction of k 



Atomic many-body calculations are required to extract the effective constant re from 
experiment. Atomic calculations of nuclear spin-dependent effects in cesium have been per- 



Note that different definitions for re, re a , K2, kq are used in different works. For example, in 
[172,177,176] the effective Hamiltonian h J e ^ = ^^ot ■ Ip( r ), re = re a + K2 + kq. Also different 
notations for the constants are used in these works. 
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formed in [175,178,179,168,16,176]. In [180] the value for k was extracted from [8] by taking 
the ratio of the nuclear spin dependent PNC amplitude to the main spin independent PNC 
amplitude, using the spin dependent calculation [179] and the spin independent calculation 
[54]. The calculations [179,54] were performed in the relativistic Hartree-Fock approxima- 
tion, with the effects of core polarization included using the time-dependent Hartree-Fock 
method and correlations included through the use of Brueckner orbitals (see Section IV). 
These calculations were performed using the same method and computer codes, so the the- 
oretical errors should cancel in the ratio. The value obtained in [180] is 

K ( 133 Cs) = 0.442(63) . (112) 

In the recent work [176] the value k( 133 Cs) = 0.462(63) was obtained. The calculations were 
performed in zeroth-order in the relativistic Hartree-Fock approximation, and core polariza- 
tion was included using the random phase approximation. This result was extracted directly 
from the spin-dependent component measured in [8], with (3 from [77]. The 4% difference 
between the two values for k( 133 Cs) is explained by correlation corrections: correlation cor- 
rections (Brueckner orbitals) are included in [179], while they are not considered in [176]. 
At the TDHF level, the results coincide. 



2. Extraction of n a 

The value of k contains three contributions (see Eq. (111)). The contributions to k from 
the nuclear spin-dependent neutral current k 2 and from the combined neutral weak charge 
and hyperfme interaction kq must be subtracted to give the anapole constant K a . In the 
single-particle approximation, k 2 = — C 2p ~ —0.05 [see Section II, Eq. (9)]. Large-basis 
nuclear shell-model calculations for 133 Cs give k 2 = —0.063 [172]. Atomic calculations for 
the combined weak charge and hyperfme interaction have been performed in [18,19,176]. The 
recent calculation [176] is the most complete: it was performed in third-order perturbation 
theory, with core polarization taken into account in the random phase approximation. In 
the other calculations [18,19] the single-particle form for the operator [Eq. (10)] was used. 
The result of Johnson et al. [176] for 133 Cs is kq = 0.017. Results of calculations for k 2 and 
kq are compiled in Table XL With k 2 from [172] and k,q from [176], the anapole moment 
constant is 

Ka = 0.368(63) . (113) 

The interaction (109) is for a point-like nucleus. However, a real nucleus has a finite size. 
In the works [178,179,168,16,176] the atomic calculations were performed by replacing 5 3 (r) 
with the nuclear density p(r). A more accurate treatment of finite nuclear size effects was 
given in [167,180]. For Cs the correction is small. The corrected value for n a [Eq. (113)] is 

n a = 0.362(62) . (114) 
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E. The nuclear anapole moment and parity violating nuclear forces 



The weak potential is usually parametrized in terms of a one-boson exchange model. The 
proton-nucleus and neutron- nucleus constants, g p and g n , can be expressed in terms of the 
following combination of meson-nucleon parity nonconserving interaction constants [21,181] 
(using the notation of Ref. [182]): 



g p = 2.0 x 10 5 W P 



g n = 2.0 x 10 5 W P 



176^A - 19.5/i° - 4.7/iJ + 1.3h 2 p - 11.3(/i° + h l J 



W, 



118-f U ~ 18.9/i° + %Ah\ - 1.3/i? - 12.8(/i° - hi) 



W, 



(115) 
(116) 



f n = h\ and the /i's are weak meson-nucleon coupling constants — the subscript denotes the 
type of meson involved and the superscript indicates whether it is an isoscalar, isovector or 
isotensor interaction (0, 1, or 2). These are the effective constants obtained in the contact 
limit of the nucleon-nucleon PNC interaction [Eq. (99)], with short-range nucleon-nucleon 
repulsion and long-range effects taken into account through the parameters W p and W^. 
Following [21,180], we take W p = 0.4 and W n = 0.16 using calculations of PNC for neutron 
and proton scattering on 4 He [183]. 

The standard reference values for weak meson-nucleon couplings are those of Desplan- 
ques, Donoghue, and Holstein [182] [we list the DDH "best values" for the weak couplings 
and the effective coupling constants (Eqs. (115,116,101); Refs. [21,184,174]) in terms of these 
values in Table XII for easy reference] . However, there are large uncertainties in the possible 
values (known as the "DDH reasonable ranges"). It is therefore important to determine the 
weak coupling constants experimentally. Information on these constants can be obtained 
from the measurement of the Cs anapole moment. 



1. The cesium result and comparison with other experiments 

The measurement of the anapole moment K a (114) and the single-particle calculation 
(108) give the following value for the weak interaction constant between an unpaired proton 
and the nuclear core [180], 

^ 5P = 6±l(exp.) . (117) 

(Note that only the experimental error is included here; there is also a theoretical error from 
the nuclear calculation of n a (108).) There is a lot of uncertainty about the value for f n . 
Following [180], we use DDH best values for h p and h u to obtain a value for f w from the 
measurement of the cesium anapole moment. Using Eqs. (117), (115) it is found that 

fS p = ( g sP _ 2) x 1.8 x 10~ 7 = [7 ± 2 (exp.)] x 10~ 7 . (118) 

(The contribution of p and uj to g p is 2.) This result (118) agrees with QCD calculations 
[185,186] which give f n = 5-6 x 10~ 7 and is in agreement with the DDH best value f n = 
4.6 x 10~ 7 (note that the DDH "reasonable range" for f w is 0.0 - 11.4 x 10~ 7 ). 
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However, there is a serious discrepancy between the weak meson-nucleon couplings ex- 
tracted from the Cs anapole moment and those extracted from other experiments. The 
following experiments are thought to give reliable information on the weak meson-nucleon 
couplings (that is, their interpretation is not hindered too much by nuclear structure un- 
certainties) [187,172]: the longitudinal analyzing power for pp scattering at 13.6 MeV [188], 
45 MeV [189], and 221 MeV [190] and pa scattering at 46 MeV [191], the circular polariza- 
tion of 7-rays emitted from the 1081 keV state in 18 F [192], and the asymmetry of 7-rays 
emitted in the decay of the state 110 keV in polarized 19 F [193]. These experiments depend 
on different combinations of the coupling constants. They are consistent and favor a small 
value of /„., roughly lying between —1 x 10 -7 and 1 x 10~ 7 (see [172,177]), in contradiction 
with the Cs anapole result. 

The value for /„- extracted above for Cs was obtained in the single particle model. Let us 
briefly discuss how this value changes when many-body effects are included. Core polariza- 
tion effects have been calculated in Refs. [170,171] using the random phase approximation. 
In Ref. [171] it was found that for 133 Cs core polarization and pairing effects decrease the 
single-particle value for K a by almost a factor of two. These calculations therefore increase 
the discrepancy of the cesium result with other experiments and with the DDH best value. 
Large-basis shell model calculations have been performed in Refs. [169,172]. These calcula- 
tions also reduce the size of K a , although not to the extent of the calculation [171]. 

In Table X results of single-particle and many-body calculations of n a for 133 Cs are 
presented. It appears that consideration of many-body effects tends to exacerbate the dis- 
agreement between the meson-nucleon couplings extracted from the Cs anapole moment and 
those obtained from other experiments. 

Note that there may be other significant corrections to the anapole moment calculations 
that may resolve the discrepancy between theory and experiment. For instance, consid- 
eration of the strong renormalization of the weak potential [168,194]. Another correction 
to the calculations could stem from the use of an improved set of strong meson-nucleon 
coupling constants. (In the one-boson exchange model the weak potential is formed from a 
product of weak and strong meson-nucleon couplings. The strong couplings g w = 13.45, 
g p = 2.79, g w = 8.37 and the isoscalar xs — —0.12 and isovector xv — 3.7 anoma- 
lous magnetic moments of the nucleon have been employed in the calculations; see, e.g., 
[182,183,172,195]. These values were used to arrive at Eqs. (115,116)). In the recent work 
[195] it was shown that using strong coupling constants in the one-boson exchange model 
constrained by nucleon-nucleon phase shifts leads to an enhancement of parity violating 
effects by about a factor of two. 

The cesium anapole moment is not the only case in which parity violating effects appear 
to be enhanced. Parity violation in the Mossbauer transitions in 57 Fe and 199 Sn nuclei give 
a value for the circular polarization of 7-rays at least four orders of magnitude larger than 
calculations (see Ref. [184]). Statistical methods were applied in Refs. [196,197] to study 
parity violating effects in polarized neutron scattering from compound nuclei. It was found 
in Ref. [197] that the experiments give parity violating effects 1.7 — 3 times larger than those 
estimated using the DDH values for weak meson-nucleon couplings. 

The Cs anapole measurement is also inconsistent with the limit on the anapole moment 
of Tl [38], 

re(Tl) = -0.22 ±0.30 . (119) 
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In fact, the Tl anapole moment itself produces problems. Its central point is of opposite 
sign to that predicted by theory. 13 Results of calculations of K a for Tl are presented in 
Table XIII (these results in terms of DDH "best values" can also be found in [173]). Using 
DDH best values for the weak couplings, the most complete many-body calculations of the 
anapole moment of Tl give « a (Tl) = 0.24 within RPA [171] and « (T1) = 0.24 [199] and 
K a (T\) = 0.10 [172] in shell model calculations. 

Note that the anapole moment of Tl extracted above (119) was obtained from the nu- 
clear spin-dependent component of the PNC optical rotation using a single-particle atomic 
calculation [198]. In a recent many-body calculation [200], in which the MBPT+CI method 
was implemented (see Section VI for a brief description and references), the central point for 
the Tl anapole moment is interpreted as K a = —0.26; that is, the disagreement with theory 
is more obvious. 

F. Ongoing/future studies of nuclear anapole moments 

It is clear that further experimental investigation is required to resolve these inconsisten- 
cies. An improved measurement of the Tl anapole moment is important, as is an anapole 
measurement involving a nucleus with an unpaired neutron (Cs and Tl have unpaired pro- 
tons). This latter experiment would give information on g n which depends on a different 
combination of weak meson-nucleon couplings that are roughly perpendicular to the current 
anapole measurements (compare Eqs. (115,116)) and so would provide a very important 
cross-check. 

We have already discussed in Sections V D, VI C the PNC experiments underway. Nuclear 
spin-dependent effects in those atoms with non-zero nuclear spin will be measured. 

Enhancement of anapole moment effects in atoms can occur due to close levels of opposite 
parity. For example, calculations [162,163] show that the PNC amplitudes in the transitions 
7s 2 % - 7s6d 3 D 1 and 7s 2 % - 7s6d 3 D 2 in 223 Ra and 225 Ra can be more than 1000 
times larger than the nuclear spin-dependent transition amplitudes in Cs. Moreover, in the 
transition 7s 2 1 5'o — 7s6d 3 D 2 there is no background nuclear spin-independent amplitude, 
since the large change in angular momentum AJ = 2 forbids it. Note that measurements 
with 223 ' 225 Ra, sensitive to the neutron constant g n , would provide a sought-after cross-check 
on the parity violating meson-nucleon couplings. 

Work towards measurements of anapole moments in diatomic molecules is underway 
[201]. Diatomic molecules are attractive for anapole moment searches because the nuclear 
spin-dependent effects, but not the (usually dominating) nuclear spin-independent ones, are 
enhanced compared to atoms due to the presence of close rotational levels of opposite parity 
[202,203]; see also the review [204]. 



13 Note that while the Oxford result [37], «(T1) = 0.23 ± 1.20 [using the single-particle atomic 
calculation [198] used to obtain (119)], has a central point that is positive, its error is four times 
larger than the Seattle value [38]. The Oxford result has been misquoted in the literature. 
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VIII. ELECTRIC DIPOLE MOMENTS: MANIFESTATION OF TIME 
REVERSAL VIOLATION IN ATOMS 



Violation of CP symmetry (combined symmetry of charge conjugation, C, and parity) 
was discovered in 1964 in the decays of the neutral K mesons [205]. It is incorporated 
into the standard electroweak model as a single complex phase in the quark mixing matrix 
(Kobayashi-Maskawa mechanism [206]). For a long time K mesons remained the only system 
in which CP violation had been observed, until just recently, in 2001, when the collaborations 
BaBar and Belle detected it in the neutral B mesons [207]. The CP violation seen there is 
consistent with the standard model predictions. 

However, a striking problem arises from cosmology. Sakharov proposed that CP viola- 
tion, present at the time of the Big Bang, is a necessary ingredient in the asymmetry of 
matter and antimatter [208]. It is well-known that standard model CP violation is insuffi- 
cient to generate the level of matter-antimatter asymmetry in the Universe. Understanding 
the origin of CP violation and searching for possible new sources is therefore a very inter- 
esting and fundamental problem. 

If CPT is a good symmetry, as it is in gauge theories, then CP violation is accompanied 
by T (time-reversal) violation. So far there has been no (undisputed) direct observation of 
T violation and its detection is of interest in its own right. Also, detection of T-violation 
may shed light on the origin of CP violation. The measurement of a permanent electric 
dipole moment (EDM) of, e.g., a neutron, atom, or molecule would be direct evidence of 
T-violation. This can be seen very simply: an EDM, d, of a non-degenerate quantum system 
is directed along the total angular momentum F (the only vector specifying the system), and 
such a correlation can only occur if both T and P are violated (d and F behave differently 
on reversal of r and t) . 

T violating electric dipole moments have an extraordinary sensitivity to new models of 
CP violation. The reason is that the standard model CP violation, appearing through just 
a single phase, is highly suppressed in such flavour-conserving phenomena. In particular, 
the standard model prediction for the neutron EDM is d n (SM) ~ 10~ 34 e cm [209], with the 
proton EDM the same order of magnitude, and the electron EDM is even smaller d e (SM) < 
10~ 38 e cm [210]. The EDMs appear only in high-order loops (3-loop diagrams for the neutron 
and 3- to 4-loop diagrams for the electron). The new sources of CP violation appearing in 
new theories generate EDMs that are many orders of magnitude larger than the standard 
model values and the size of these EDMs are in reach of current experiments. Already 
the parameter space of popular theories such as supersymmetry, multi-Higgs models, and 
left-right symmetric models are strongly restricted by current measurements. 

A CP-violating phase also appears in quantum chromodynamics (QCD) and it can be 
constrained (or detected) from EDM measurements. It defines the strength of the P- and 
T-violating term in the QCD Lagrangian containing G^G^ (this is analogous to F^F^ 
in electrodynamics which reduces to E • B, clearly a P, T-violating correlation; see, e.g, 
Ref. [12]). The "unnatural" smallness of this phase (Oqcd ~ 10~ 10 , compare this to the 
Kobayashi-Maskawa angle 5km ~ 1) is the famous "strong CP problem". One solution 
to this problem involves the introduction of a new Goldstone boson known as the "axion" 
[211], and this particle has become a popular candidate for dark matter. However, after many 
experimental searches, there is no evidence of its existence. (For a review of experimental 
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searches for axions, see, e.g., [212].) 



A. Atomic EDMs 

In this section we will consider the contributions to an atomic EDM arising from various 
P, T-odd mechanisms. Contributions to an atomic EDM arise from (i) the sum of intrinsic 
EDMs of the atomic constituents (averaged over the atomic state) and (ii) the mixing of 
opposite parity wave functions due to a P, T-odd interaction Hpt- 

The atomic EDM induced in an atomic state K due to admixture with the opposite- 
parity wave functions M has the form 

d _ = 2E (MKj^> = d ^ (F/F) (l2Q) 

where D = — e J2i r i is the electric dipole operator, Hpt is the P, T-odd operator that mixes 
K with the set of wave functions M, and F is the total angular momentum of the atom 
corresponding to the state K. 

We will use the following notations for the electron wave functions: 

*<«> = (V- ( &0 • (121) 

where flji m is a spherical spinor, n = R/P, and f(R) and g(R) are radial functions (see, 
e.g., [12]). 



1. Electronic enhancement mechanisms 

What are the mechanisms that lead to enhancement of parity and time invariance vio- 
lating effects in atoms? Like the P-odd, T-even effects we discussed in the first part of this 
review, the effects of P, T- violation increase in atoms with: 

(i) a high nuclear charge Z. P, T-odd effects in atoms increase rapidly with Z, [213,214]. 
Therefore, it is best to search for atomic EDMs in heavy atoms. 

(ii) close levels of opposite parity. From Eq. (120) we can see that if H PT mixes oppo- 
site parity levels with energies E\ ps E 2 , then the atomic EDM induced will be enhanced, 
c^atom oc l/(Pi — E2) . This enhancement occurs, for example, in rare-earth atoms, where 
there are anomalously close levels of opposite parity. 



B. Enhancement of T-odd effects in polar diatomic molecules 

Sandars was the first to point out the enhancement of T-odd effects in polar diatomic 
molecules compared to atoms in his proposal to measure the proton EDM in TIF [215]. 
The idea is to polarize the molecules along an external electric field, thereby aligning the 
enormous intramolecular field, leading to an effective enhancement of the external field by 
several orders of magnitude. This enhancement mechanism is related to that arising from 
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the presence of close rotational levels of opposite parity [203] . The two cases correspond to 
the respective cases of strong and weak interactions with the electric field compared to the 
rotational level spacing. See [204] for a review. 

C. Limits on neutron, atomic, and molecular EDMs 

A typical EDM experiment is performed in parallel electric and magnetic fields. The 
corresponding Hamiltonian is 

h = -fi-B-d E . (122) 

A linear Stark shift is measured by observing the change in frequency when the electric field 
is reversed (this is a measure of the P, T-odd correlation E • B). 

To date, permanent EDMs in neutrons, atoms, and molecules have escaped detection. 
Null measurements of EDMs have been obtained for the paramagnetic atoms 85 Rb [216,217], 
133 Cs [218-220], 205 T1 [221-224,11], and 129 Xe in the metastable state 5p 5 Qs 3 P 2 [225] and, 
very recently, the molecule YbF in the ground state [226] . An experiment aimed to measure 
the EDM of the ion Fe 3+ in the solid state was carried out a long time ago [227]. Interest to 
carry out measurements in solids, in particular in gadolinium gallium garnet and gadolinium 
iron garnet, has been sparked recently [228,229] (see Section X). In diamagnetic systems, 
EDM experiments have been performed for atomic 129 Xe [230,231] and 199 Hg [232,233,10] in 
the ground states and for the polar molecule TIF [234-238]. 

Measurements of EDMs in paramagnetic systems (that is, with non-zero total electron 
angular momentum) are most sensitive to leptonic sources of P, T- violation, in particular 
the electron EDM, while measurements of EDMs in diamagnetic systems (zero total electron 
angular momentum) are most sensitive to P, T-odd mechanisms in the hadronic sector. 
Both are sensitive to P, T-odd semi-leptonic processes (the electron-nucleon interaction), 
the former to those involving the electron spin, and the latter to those involving the nuclear 
spin. 

The most precise measurement of an atomic EDM in a paramagnetic system has been 
obtained for 205 T1 [11], 

rf( 205 Tl) = -(4.0 ± 4.3) x 10" 25 e cm . (123) 

With the use of atomic calculations, described in Section VIII D 2, this measurement can 
be expressed in terms of a limit on the EDM of the electron, for which it currently sets the 
tightest constraint. 

The most precise measurement of an atomic EDM has been carried out with 199 Hg [10], 

d( 199 Hg) = -(1.06 ± 0.49 ± 0.40) x 10~ 28 e cm . (124) 

This measurement sets the best limits on a number of P, T- violating mechanisms in the 
hadronic sector. The nucleus of 199 Hg has an unpaired neutron. A limit on the neutron 
EDM extracted from (124) is competitive with the best direct neutron EDM measurements 14 



14 In line with the comments of Ref. [239], we present the most recent 1999 result of the ILL group, 
rather than the final value cited in Ref. [240] which is the average of their 1999 and 1990 results. 
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performed at the Institut Laue-Langevin (ILL) and the Petersburg Nuclear Physics Institute 
(PNPI) 



d 



a, 



(1.9 ±5.4) x 1(T 26 e cm 
(2.6 ±4.0 ±1.6) x 10~ 26 



e cm 



ILL [240] 
PNPI [241] ' 



(125) 



As we will see in Section X, at the most fundamental scale the limits on P, T- violating 
parameters from the measurement (124) are more strict than those from the neutron mea- 
surements (125) due to the presence of a nucleon-nucleon P, T- violating interaction that 
induces P, T- violating nuclear moments more efficiently than does an intrinsic EDM of a 
nucleon. 

It is interesting that although the 199 Hg nucleus has an unpaired neutron, the measure- 
ment (124) sets a constraint on the proton EDM (it contributes due to configuration mixing) 
that is tighter than that from the measurement of the EDM of the TIF molecule (Tl has an 
unpaired proton) used in the past to set the proton EDM upper limit. In TIF the P, T-odd 
interaction is given by H = —hdcr ■ A, where cr is the spin operator of the Tl nucleus, A 
is the unit vector along the internuclear axis, and h is Planck's constant. The most precise 
limit on the P, T-odd coupling constant d in TIF is [238] 



The coupling constant d can be expressed, e.g., in terms of permanent EDMs of the proton 
and electron and P, T- violating interactions; such calculations for TIF have been performed 
in [215,242-244,165,18,245-247]; see, e.g, Ref. [238] for a general overview. The limit on the 
proton EDM following from the calculations of Ref. [247] is presented in Table XIX. 

A more detailed comparison of limits on P, T-violating parameters will be postponed 
until Section X. In the following sections we discuss the different mechanisms that induce 
atomic EDMs. 



An external electric field acting on a neutral atom consisting of non-relativistic point-like 
charged particles with EDMs, interacting via electrostatic forces, is screened exactly at each 
particle [248,84]; see also [12]. 15 This occurs due to polarization of the atomic electrons by 
the external field. An atomic EDM cannot be induced in such a case. However, as shown 
by Schiff [84], if magnetic or finite-size effects are taken into account, there is incomplete 
shielding, and so atomic EDMs can in principle be measured. 

An atomic EDM can be induced from the following P, T-odd mechanisms: 

(i) an intrinsic EDM of an electron - the electron EDM can interact with the atomic field 



15 Actually, for moving particles like electrons the average value of the screened field is zero. An 
explanation is the following: a neutral atom is not accelerated by the homogeneous external field. 
If a charged particle inside the atom is not accelerated the electric field acting on this particle is 
zero. Therefore, the EDM of the particle has nothing to interact with, d • (E) = 0. 



d= -(0.13 ±0.22) mHz . 



(126) 



D. Mechanisms that induce atomic EDMs 
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producing an atomic EDM that is many times larger than the single electron EDM [213]; 

(ii) a P, T-odd electron-nucleon interaction [249,242]; 

(iii) an intrinsic EDM of an external nucleon [215,250]; 

(iv) a P, T-odd nucleon-nucleon interaction - this interaction can induce P, T-odd nuclear 
moments that can greatly exceed the moments of single nucleons [251,165]. 

The latter two mechanisms (iii) , (iv) can be grouped together at the nuclear scale since 
they both produce P, T-odd nuclear multipole moments. 

One may consider other exotic T-violating mechanisms such as dyon 16 vacuum polariza- 
tion [252]. 

There are in fact other mechanisms that can lead to an atomic EDM that are not P, T- 
violating at the fundamental level. For example, T-odd, P-even interactions can have P, T- 
odd effects due to P-odd radiative corrections (see, e.g., [13]). 

In Fig. 22 we present a flowdiagram (slightly modified from the review Ref. [253]) show- 
ing the CP-violating mechanisms at different energy scales that induce neutron, atomic, 
and molecular EDMs. Read from left to right, it is seen clearly which measurements con- 
strain which CP- violating parameters at smaller distances (and which popular CP- violating 
models). (Read in the other direction, it is seen which small-distance mechanisms induce 
large-distance CP- violating effects.) Calculations are required to relate the CP- violating 
parameters at different scales. Solid lines indicate the parameters that are most strongly 
constrained (induced) by the parameters to the left (right), while dashed lines indicate a 
weaker constraint (inducing mechanism). 

In this review we focus primarily on atomic EDMs induced by P, T-odd nuclear mo- 
ments that originate from P, T-odd nuclear forces, since there have been several recent 
developments in this area. We will first discuss atomic EDMs induced by the P, T-violating 
electron-nucleon interaction in Section VIII D 1 and by the electron EDM in Section VIII D 2. 
For the sake of completion, we briefly introduce nuclear P, T-violating moments in Section 
VIII D 3 before a comprehensive overview of these moments is given in Section IX. 

1. The F ', T '-violating electron-nucleon interaction 

The P, T-violating electron-nucleus interaction has the following form (see, e.g. [13]): 
G 

N 



h = ]T \C s N p NNe lb e + C PS N l5 Nee + C^N^a^Nea^e] . (127) 

V 2 AT 



The real, dimensionless constants Cff ', C PS ', and C^ give the strength of the scalar- 
pseudoscalar, pseudoscalar-scalar, and tensor P, T-odd electron-nucleon interactions for the 
nucleon N. Upper limits on the constants C SP , C PS , and C T can be obtained from mea- 
surements of atomic EDMs. 

In the limit of an infinitely heavy nucleon the following form for the electron-nucleus 
interaction is obtained [13] 



16 A dyon is a particle with both electric and magnetic charges. 
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h SP > T = h SP + h p = i^=5(r) [{ZC S / + NC s n p ) lol5 + 2{C T p £ <r p + C T n £ a n ) ■ 7] . (128) 

V ^ p n 

In this approximation the term containing C PS vanishes. 

Notice the similarity between this expression (128) and Eqs. (4), (8). However, here the 
matrix element is real (the factor % is placed here to make the operator Hermitian; another 
factor of % arises due to the mixing of the upper and lower components of the electron wave 
function (121) due to 75 and 7$). Therefore the interaction (128) mixes atomic states of 
opposite parity and induces static electric dipole moments in atoms. 

Like their P-odd, T-even analogues (4) and (8), the nuclear spin-independent term in 
(128) receives coherent contributions from the nucleons inside the nucleus, whereas the 
smaller second term, dependent on the nuclear spin, arises from the unpaired nucleons. 

Each of the interactions H SP and H T can induce EDMs in paramagnetic atoms. However, 
only H T can open up the closed electron shells of diamagnetic atoms and thus induce an 
EDM; this term is dependent on the nuclear spin, while the scalar-pseudoscalar term is not. 
Therefore, the interaction H SP cannot by itself contribute to EDMs in diamagnetic atoms. 

However, by allowing for the hyperfine interaction, measurements of EDMs of diamag- 
netic atoms can place limits on C SP which are, in fact, as competitive as those obtained 
from experiments with paramagnetic atoms. The atomic EDM induced by the P, T-odd 
scalar-pseudoscalar electron-nucleus interaction H SP along with the hyperfine interaction 
Hhf arises in the third order of perturbation theory [18], 

datom = E ^ ' ' ~ ' Z j T7^ + Permutations . (129) 

mn I A) - An MA) ~ &n) 

The matrix element of an effective operator constructed from H SP and Hhf, 



(g«'|B)(B|g w ) + (g v |n)(n|g"-) 

^eflf = p p > 

n A — A 

has the form (pi/2|-H'Jf'l s i/2) °^ (ZCp P + NC^ p )j ■ I. This matrix element can be related to 
that of (pi/2\H T \si/2) oc j • (CpJ^p&p + C« En^n), the brackets ( ) denote averaging over 
the nuclear state. The correspondence between rf at om(C 5P ) and rf at om(C T ) [18,254,13], 

(fqf P + ^Oy-1.9x^ , (131) 

where /x is the nuclear magnetic moment in nuclear magnetons, can be used to obtain the 
sensitivity of c4tom to C SP from calculations of d atom (C T ) for diamagnetic atoms. 

Calculations of atomic EDMs induced by H SP and H T are presented in Table XIV. 

The current best limits for C SP have been obtained from the 205 T1 and 199 Hg measure- 
ments, Eqs. (123) and (124), 



(0.40C- + 0.60O - { ^sfo,) x 10- Hg • < 132 ' 
and for C T from Hg, 
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C T n = -(5.3 ±2.5 ±2.0) x 1(T 9 . (133) 

Here we have used the simple shell model of the nucleus, (cr n ) = —(1/3)1/7. 

Let us now consider the contribution to an atomic EDM arising from the pseudoscalar- 
scalar component of the electron-nucleus interaction (127). In the lowest non- vanishing 
approximation in m" 1 , the Hamiltonian of the electron- nucleus interaction reduces to the 
form [13] 

h PS = -^^-K S X>p + C P n S £ OW(r)7o • (134) 

Again, the matrix element of this interaction (p\/2\H PS \si/2) oc j • (C PS : J2 P < 7 p + C PS : I] n c n ) 
can be related to that of H T . The correspondence is [18,12,13] 

A 1 / 3 

C PS <-> 4.6 x 10 3 ^— C T . (135) 

Zj 

See Table XIV for calculations of the sensitivities of atomic EDMs to C PS . It is seen 
that this interaction induces EDMs much less efficiently than H SP and H T . 

It is the 199 Hg EDM measurement (124) that currently places the tightest constraint on 

C PS , 

C PS = -(1.8 ±0.8 ±0.7) x 1(T 6 . (136) 



2. The electron EDM 

The best limits on the electron electric dipole moment are derived from measurements of 
atomic EDMs. Salpeter first noted the possibility of an enhancement of the electron EDM 
in atoms through consideration of the metastable 2s state in hydrogen [258]. 

As noted above, when magnetic effects are considered, the screening of the EDMs in an 
atom (Schiff theorem) is lifted [84]. Sandars [213] pointed out that due to relativistic mag- 
netic effects, the atomic EDM induced in heavy atoms can be strongly enhanced compared 
to the electron EDM inducing it. The value of the atomic EDM compared to the electron 
EDM is expressed through an enhancement factor 

K = d atom /d e . (137) 

The enhancement factor K increases with nuclear charge Z faster than Z z . 

We will look briefly at how an electron EDM induces an EDM of an atom as a whole; 
for a more detailed consideration we refer the reader to the works [259,214,12,13,260]. As 
we mentioned at the beginning of Section VIII A, there are two types of contributions to 
an atomic EDM arising from constituent EDMs. For the case of electron EDMs these 
are the following: (i) the sum of the intrinsic EDMs of the electrons (0|4 E 4 =i 7o^|0); 
and (ii) the admixture of opposite parity atomic states due to the pseudoscalar interaction 
H e = —d e Yh=\ 7o^* ' E mt ; see Eq. (120). Here 7 and £ = 7o7s7 are Dirac matrices defined 
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in Eq. (2), |0) is the unperturbed state of the atom (it is an eigenstate of the P, T-even 
Hamiltonian with no external electric field), and E int is the internal atomic electric field. 

Let us consider for a moment the Stark shift generated by the presence of the electron 
EDMs. It has the form AE = (0| - 4E*=i7o s * • Ei |0), where |0) is an eigenstate of the 
P, T-even Hamiltonian H which includes the external electric field and E = E int + E cxt is 
the total electric field. (When the external field is treated perturbatively and only terms 
first-order in this field are considered, this reduces to the linear Stark shift generated by 
the two contributions to the atomic EDM mentioned above.) It is convenient to break up 
the pseudoscalar interaction 7 £ • E = E ■ E + (70 — 1)£ • E, since the first term on the 
right-hand-side does not contribute to a linear Stark shift, S • E = (l/e)[S • V, H}. 

The enhancement factor for the atom can then be written as 

KWJ) = (o| XK - de;|o> + 2eE <^£M^LgM^kiW _ m 

i=l M M 

where J is the electron angular momentum of the state |0). The operator 

^°- i ) s =(o-L) < i39 > 

mixes the lower components of the wave functions (121). So it is seen that this is a purely 
relativistic effect. In heavy atoms the first term in Eq. (138) is small compared to the second 
and so can be omitted in the calculations. 

The first analytical formulation of the enhancement factor was performed in Ref. [214]. 
The following expression was obtained for alkaline atoms, in terms of quantities that can be 
determined experimentally, 

^ = V 4(Za) 3 r 0m hc 

m (J + 1K7(4 7 2 - l){w m ?' 2 {E m -E Q )- { U> 

The sum over m is taken over the excited states of the external electron, z/ , u m are effective 



principal quantum numbers for the ground and excited states, 7 = y (J + 1/2) 2 — (Za) 2 , 
and r 0m is the electric dipole radial integral. Taking into account mixing with only the 
nearest level, and assuming the values r 0m = 5a#, u Q = u m = 2, and E m — E = (l/10)Ry, 
we obtain: 

1*1 ~ ' V+Su+i)' * ' (141) 

where R is a relativistic enhancement factor that increases with Z and is 1.2 for Rb and 2.8 
for Fr in the ground states (R tends to unity when J is large). 

This simple formula (141) illustrates the dependence of the enhancement factor on Z 3 
and on the angular momentum J and can be used to obtain order of magnitude estimates. 
It is seen that K is large for high Z and low J. 

For atoms with more complex configurations these formulae are not applicable and nu- 
merical calculations of the enhancement factor are required. 

As we mentioned earlier, experiments with paramagnetic atoms and molecules are most 
sensitive to the electron EDM. The current best limit on the electron EDM comes from the 
Tl measurement (123). The sensitivity of atomic thallium to the electron EDM is [261] 
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d( 205 Tl) = -5854 (142) 

and accordingly the measurement (123) of the electron EDM is 

d e = (6.9 ± 7.4) x KT 28 e cm . (143) 

There is some sensitivity of diamagnetic systems to the electron EDM [262], al- 
though this sensitivity is very weak. The dominant contribution appears in third-order 
perturbation theory due to consideration of the hyperfine interaction, Eq. (129) with 
H SP replaced with the interaction H e [18]. (The second-order "bare" contribution, 
2 l]ra(0| — d e Yh=i 7o^*l n ) ( n I 1 0) / (-Eo — E n ), is significantly smaller and can be neglected 
[18].) Another contribution [18] comes from the direct interaction of the nuclear mag- 
netic field B (arising from the nuclear magnetic moment) with the electron EDM [258], 
Hpt = —id e Y^i = il % • B\ This latter interaction contributes in second-order, Eq. (120). 
The following relation between d e and the P, T-odd tensor electron-nucleon interaction has 
been obtained [18,12,13], 

^'-i^T^r-v^^' (i44) 

R is a relativistic enhancement factor. Using the result of Ref. [256] for the calculation of 
^atom(C T ) for 199 Hg and the above relation, the enhancement factor fT( 199 Hg) = —0.014 was 
found in Ref. [18]. A TDHF self-consistent calculation performed in [263] yielded the result 

d( 199 Hg) = 1.16 x 10" 2 d e . (145) 

Due to huge polarization corrections, it is of opposite sign, and the same order of magnitude, 
as the lowest-order result, and its value may change with inclusion of correlation corrections 
[263]. 

While the enhancement factors are small for diamagnetic systems, the extraordinary 
precision that has been achieved in the Hg measurement (124) makes the corresponding 
measurement of the electron EDM, d e — — (9.1 ± 4.2 ± 3.4) x 10~ 27 e cm, comparable with 
those from the best paramagnetic EDM measurements. 

In Table XV we list enhancement factors for both paramagnetic and diamagnetic atoms 
of experimental interest. 

Paramagnetic polar diatomic molecules are attractive for electron EDM studies, in par- 
ticular those with electron states 2 Si/2 and 2 U 1 / 2 (see, e.g., [203,264] and the review [204]). 
In [203] an analytical estimate for the energy shift in such molecules (Q — 1/2) was made, 



4(7o - 1)S • E tot \P) = K d a ■ A , |« d | ~ 77J73 ^ 2 , (146) 



Z 3 a 2 ed e 
7(4 7 2 - 1)4 



where \/3) are molecular orbitals built up from atomic orbitals mixed by the strong internal 
electric field. Molecules with electron ground state 2 Si/ 2 include BaF, YbF, HgF, PbF. The 
same estimate (146) is valid for the metastable a(l) 3 S state in PbO with which an EDM 
experiment is in progress [270]; more refined calculations can be found in Refs. [271,272]. 
The recent measurement of the electron EDM in YbF yielded the result [226] 
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d e = (-0.2 ±3.2) x 1(T 26 e cm 



(147) 



(calculations for the effective electric fields have been performed in Refs. [273]). This is 
the first measurement of an EDM in a paramagnetic molecule, and while the limit on the 
electron EDM is not as impressive as that from Tl or even Hg, it is limited only by statistics. 

3. P, T -violating nuclear moments 

Atomic EDMs can be induced if the nucleus possesses P, T-odd nuclear moments. These 
moments arise at the nucleon scale due to a P, T- violating interaction between nucleons or 
due to intrinsic nucleon EDMs. The induced nuclear moments can be electric or magnetic. 
For example, the following moments violate parity and time-reversal invariance: electric 
dipole, magnetic quadrupole, electric octupole. For the electric case, the interaction Hamil- 
tonian that mixes opposite parity electron states, and induces an EDM of the atom, is of 
the form hpr = — ey?, where ip is the electrostatic potential of the nucleus corresponding to 
a P, T-odd charge distribution. In fact, due to Schiff's theorem [84], there is an additional 
screening term which we will look at in Section IX. In the magnetic case Schiff's theorem is 
not valid, and the interaction Hamiltonian of a relativistic electron with the vector potential 
A corresponding to a P, T-odd current distribution is simply hpx = ecx ■ A. See formula 

( 12 °)- 

The operator hpr has electronic and nuclear components. While the overall operator hpr 
is a scalar, the electronic and nuclear operators can be of any (equal) rank. Accordingly, the 
triangle rule for addition of angular momenta imposes restrictions on the angular momenta 
of the electron and nuclear states for non-zero matrix elements. For example, the electron 
interaction with the nuclear magnetic quadrupole moment cannot mix s and pi/2 electron 
states, since we must have \ji — jz\ < 2 < ji + j 2 . Similarly, a static magnetic quadrupole 
moment of the nucleus cannot arise in nuclei with total angular momentum / < 1. 

We leave a detailed consideration of P, T-odd nuclear moments for the next section. 

IX. P, T- VIOLATING NUCLEAR MOMENTS AND THE ATOMIC EDMS THEY 

INDUCE 

This section is devoted to a consideration of the P, T-odd nuclear moments that can 
induce atomic EDMs. In Sections IX A, IX B we look at the form of the P, T- violating 
electric and magnetic moments. In Section IX C we discuss how P, T- violating nuclear 
moments are induced by P, T- violating mechanisms at the nucleon scale. Enhancement 
mechanisms for nuclear moments are reviewed in Section IX D. Finally, in Section IX E we 
look at calculations of atomic EDMs induced by nuclear moments. 

A. Electric moments; the Schiff moment 

When considering the P, T-odd electric moments of the nucleus, we must take note of an 
important screening phenomenon - the Schiff theorem (we mentioned this at the beginning 
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of Section VIII D). The electron screening is taken into account by using the following 
(screened) electrostatic potential of the nucleus (for a derivation, see, e.g., [274]): 

* e P( r ) j3„ , 1 ,a tt\ f P( r ) A 



^ R) = J^ dr + z<- d ' v) Jwr7r- (148) 

where p(r) is the nuclear charge density, / p{r)d 3 r = Z, and 

d = J erp(r)d 3 r = dl/I (149) 

is the P, T-odd nuclear EDM. 17 

We are interested in the contributions of the first and second terms to ip that are first 
order in the P, T-odd interaction. The first term on the right-hand-side of Eq. (148) is 
P, T-odd if the charge density is distorted due to a P, T-odd interaction. The density in the 
second term can be considered spherical, since the nuclear EDM is P, T-violating (it arises 
due to a P, T-violating component of the density in Eq. (149)). 

If we consider the nucleus to be point-like, then we can perform a multipole expansion of 
the potential (148) in terms of r/R. According to Schiff's theorem, the nuclear electrostatic 
potential is screened by atomic electrons such that the dominant nuclear P, T-odd moment, 
the nuclear EDM, of a point-like nucleus cannot generate an atomic EDM. It is easy to see 
this from Eq. (148): 

- J ep(r)(r • V-^)d 3 r + i(d • V)-^ J p{r)d 3 r = . (150) 

The first non-zero P, T-odd term in Eq. (148) is then 

V? {3) = - \ J ep{r)r a r p r y d 3 rV a V p V^ + ^(d • V)V Q V /3 -^ J p{v)r a r p d 3 r . (151) 

The first term r a r^r 7 on the right-hand-side of the equation is a reducible rank-3 tensor, 
while the second r a rp is a reducible rank-2 tensor. Separating the trace, 



r a r e ,r 1 - -r 2 (r a 6fa + r p 5 ai + r 7 <5 Q/3 ) + ^r 2 (r a 5 Pl + r p 5 ai + r 7 <5 a/3 ) (152) 



1 



1 



r a rp - -r 5 a/ 3 + -r 5 aj3 , (153) 



3 



it is seen that tp^ is comprised of a rank-3 octupole potential <Poctupoic and a rank-1 "Schiff" 
potential y?Schifr, 



y? (3) = y^ctupoic + (pschifr , (154) 



17 The screening term appears as a result of a unitary transformation of the Hamiltonian which 
does not change the linear Stark shift. For exact atomic wave functions the result for the atomic 
EDM must be the same with and without the screening term. However, in real (approximate) 
calculations inclusion of the screening term is a must. 
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where 

V^octupoic = -^0 Q /3 7 V a V/3V 7 4 + -TTzQapid • V)V Q V /? 4 (155) 
o it e Jx 

V^schis = 4ttS • V5(R) , (156) 

and we have used V 2 (l/-R) = — 47r5(R). S is the P, T-odd nuclear Schiff moment, O a ^ is 
the P, T-odd nuclear electric octupole moment, and Q a p = J ep{v){r a rp — ^r 2 5 a p)d 3 r is the 
P, T-even nuclear electric quadrupole moment. The second term in Eq. (155) is small since 
only protons in the external shell contribute to Q a p and there is a factor -|. The nuclear 
octupole and Schiff moments are given by 

f 1 

O aM = J ep(r) [r Q r^r 7 - -r 2 (r a S Pl + r p S aj + r 7 <5 a/3 )j d 3 r , (157) 
Y ep(r)rr 2 d 3 r - ^d^ ^ p(r)r 2 d 3 r] = SI/ 1 . (158) 



s = ± 

10 



Because the octupole moment carries 3 units of angular momentum it can only 
arise in nuclei with spin I > 3/2, whereas the Schiff moment can arise in nuclei with spin 
I > 1/2 (due to the triangle rule for the addition of angular momenta). Of the atomic EDM 
measurements performed so far, only Cs (I = 7/2) has a nuclear spin large enough to have 
a static octupole moment; all other nuclei have spin / = 1/2. 

However, for these moments to induce an atomic EDM they must satisfy electronic 
angular momentum requirements. Due to the higher rank of the octupole moment it mixes 
electronic states of higher angular momentum than the Schiff moment. This means that the 
atomic EDM induced by the octupole moment is smaller than that induced by the Schiff 
moment because the wave functions of electrons with higher angular momentum penetrate 
the vicinity of the nucleus less due to the greater centrifugal barrier. The lowest value for 
the angular momentum of the electrons that can induce an atomic EDM due to mixing by 
the octupole moment is j = 3/2. [The conditions imposed on the allowed electronic angular 
momentum for mixing by the octupole moment is \j\ — j'2] < 3 < ji + j 2 and that allowed 
by the electric dipole mixing is \j± — j 2 ] < 1 < ji + 32, so the conditions for inducing an 
atomic EDM are \ji — j 2 ] < 1 and ji + j 2 > 3.] This means that s states cannot contribute 
to the EDM produced by the electric octupole moment, so that in fact the octupole moment 
of the Cs nucleus cannot induce an atomic EDM in Cs in the ground state, as this state 
corresponds to a configuration with a single electron in an s-state above closed shells. Static 
nuclear octupole moments and the atomic EDMs they induce have been considered in detail 
in Ref. [275]. The EDMs they induce in atoms are very small, so we will consider them no 
further. 

It is therefore obvious that the Schiff moment is the only P, T-odd moment that induces 
an EDM in atoms with closed electron subshells such as Xe and Hg. (The nuclear magnetic 
quadrupole moment, which will be discussed in the next section, cannot induce an atomic 
EDM in systems with zero electron angular momentum, since there is no magnetic field 
of the electrons for the MQM to interact with. The same conclusion also follows from the 
triangle rule applied to Eq. (120).) In fact, all the atoms for which EDM measurements have 
been performed ( 133 Cs (I = 7/2, J = 1/2), 205 T1 (I = 1/2, J = 1/2), 129 Xe 3 P 2 (I = 1/2, 

J = 2), 129 Xe 1 S (I = 1/2, J = 0), 199 Hg (I = 1/2, J = 0)) can have contributions from 
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the nuclear Schiff moment, however it is only Cs that can have an EDM arising due to a 
static magnetic quadrupole moment. 

Let us consider the form of the atomic EDM (120) induced by the interaction of electrons 
with the Schiff moment (156). The contact interaction H PT = — e J2i V^schifr mixes s- and 
p-wave electron orbitals and produces EDMs in atoms. The expression (156) is consistently 
defined for non-relativistic electrons. Using integration by parts, it is seen that the matrix 
element (s\ — e</?s c hiff|p) is finite, 

(s\ - ey^schiff b) = 47reS • (Viplifj p ) R= o = constant . (159) 

However, atomic electrons near the nucleus are ultra-relativistic, the ratio of the kinetic 
or potential energy to mc 2 in heavy atoms is about 100. For the solution of the Dirac 
equation, (V^|V ; p)r^o ^ oo for a point-like nucleus. Usually this problem is solved by 
a cut-off of the electron wave functions at the nuclear surface. However, even inside the 
nucleus Vip\ip p varies significantly, ~ Z 2 a 2 , where a is the fine-structure constant, Z is the 
nuclear charge. In Hg (Z = 80), Z 2 a 2 = 0.34. 

A more accurate treatment requires the calculation of a new nuclear characteristic which 
has been termed the local dipole moment (LDM) [276]. This moment takes into account 
relativistic corrections to the nuclear Schiff moment which originate from the electron wave 
functions. So in the non- relativistic limit, Za -> 0, the LDM L = S. For 199 Hg, L w 
S(l — 0.8Z 2 a 2 ) pa 0.75S. When considering the interaction of atomic electrons with the 
LDM it is defined as placed at the center of the nucleus, that is the electrostatic potential is 

(p(R) = 4ttL • V5(R) . (160) 

See Ref. [276] for the explicit form for L. 

It is more convenient to use a real electric field distribution produced by a P, T-odd 
perturbation. In [276] it was shown (by considering several nuclear models) that the natural 
generalization of the Schiff moment potential for a finite-size nucleus is 

^( R ) = -^^ n ( R - R n) , (161) 

where is the nuclear radius and n(R—R]y) is a smooth function which is 1 for R < Rjy — 5 
and for R > R N + 5; n(R — R N ) can be taken as proportional to the nuclear density. This 
form for the electrostatic potential has no singularities and is suitable for relativistic atomic 
calculations. 



1. The P,T-odd electric field distribution in nuclei created by the nuclear Schiff moment 

From the new form for the electrostatic potential Eq. (161) it can easily be seen that 
the Schiff moment gives rise to a constant electric field inside the nucleus (see Fig. 23), 
E = — Vy?. The correlation between the electric field and the nuclear spin, E oc I, is 
naturally P, T-odd. This electric field polarizes atomic electrons, producing an EDM of the 
atom. 
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B. Magnetic moments; the magnetic quadrupole moment 



In the gauge V • A = the vector potential produced by a steady current is 

a(R) = / |^r^ 3 r > ( 162 ) 

where j is the vector current density. The lowest-order term in the multipole expansion of 
Eq. (162) is the P, T-even magnetic dipole moment, and we have no interest in this. The 
lowest-order P, T-odd moment arises in second-order and is the rank-2 magnetic quadrupole 
moment (MQM); it appears alongside the P-odd, T-even anapole moment a (see Section 
VII) [165], 

4 2) = \j hr S r a d 3 rV 5 V a ^ (163) 
1 lnl 

— (8 ia a 5 - 5 Sa a 1 ) - -e^M^J V«jV a — (164) 

= A a +A MQM^ (165) 

where the anapole moment is given by Eq. (97) and the MQM is 

M a/3 = - J (r a e Kv + rpe a ^ v )j^r v d 3 r . (166) 

The P, T-odd component of the current density j will give rise to a non-zero magnetic 
quadrupole moment. Its form is specific to the P, T-odd mechanism creating it. For instance, 
if we consider that it is produced by an external nucleon perturbed by P, T-odd nuclear 
forces, then we can use the current (103), 

e r r 2 i 

M a(3 = — J [3fi(r a ap + rpa a - -5 a(3 cr • r) + 2q(rjp + rpl a )\p(r)d 3 r . (167) 

The magnetic quadrupole moment can now be calculated using a P, T-odd perturbed density 
P(r). 

A general expression for the MQM can be constructed in terms of the total angular 
momentum of the system I, 



M, 



M 



a/3 



21(21-1] 



I a Ip + I p I a -^I(I + l)S aP ] . (168) 



The quantity M is conventionally referred to as the MQM and is defined as the maximum 
projection of M a p on the nuclear axis, M = M zz . It is easily seen from a comparison of 
Eqs. (166), (168) that the magnetic quadrupole moment violates parity and time- reversal 
invariance. 

The interaction of electrons with a nuclear MQM induces an atomic EDM typically an 
order of magnitude larger than that induced by the nuclear Schiff moment [250,165]. The 
ratio of the s — p electronic matrix elements is [165] 

<^,- AMQ ">> ~ ioM-/3^ , (W9) 

(s|y?schifr|p) Rs 
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Rm, Rs are relativistic enhancement factors for the magnetic quadrupole and Schiff mo- 
ments, respectively, R — > 1 as Za — > 0. It is seen that for atoms with light nuclei the 
contribution of the MQM dominates. The relativistic factor Rs grows faster than Rm with 
increase of Z. For example, at Z — 80, Rs = 7 [pi/2],5 [P3/2] while Rm = 1-8; explicit 
formulae can be found in Ref. [165]. (In the square brackets the angular momentum of the 
p electron state is specified; for the MQM, p = ps/2-) At Z = 80, A = 200, the ratio (169) 
reaches ~ I. 

1. The spin hedgehog 

In a spherically symmetric system the P, T-odd interaction induces a "spin-hedgehog" 
whereby the spin density is proportional to the radial vector, cr oc r [277,278]. The P, T- 
odd nucleon-nucleon interaction (leading to the perturbed wave functions (177), see below) 
produces the following distributions of the spins for protons and neutrons in the nucleus, 

<r p (r) = £ P V Pp (r) , <r n (r) = £ n Vp n (r) , (170) 

the unperturbed nuclear density p = TJ IV'I 2 - This collective spin distribution, however, 
produces no current (j(r) = /iV x cr(r) oc V x Vp(r) = 0) and hence no magnetic field 
[278]. One may think that because the spin-hedgehog has no magnetic field it produces no 
effects. This is not the case. The spin-dependent part of the strong interaction is sensitive 
to this spin structure. It reduces the constants of the P, T-odd nucleon-nucleon interaction: 
for the case of an external proton interacting with the spin-hedgehog, f] p — > r] p /1.5, while for 
a neutron, i] n — > i] n /1.8 [278,194]. A distorted spin hedgehog in deformed nuclei produces a 
collective magnetic quadrupole field (see Section IX D 2). 

Note that the spin-hedgehog is not specific to nuclei. For example, the P, T-odd electron- 
nucleon interaction (Section VIII D 1) mixes atomic states of opposite parity and induces a 
spin-hedgehog of the atom [277]; see, e.g., [12,13] for details. 

C. What mechanisms induce P, T-odd nuclear moments at the nucleon scale? 

P, T-odd nuclear moments can arise due to an intrinsic EDM of an external nucleon or 
due to P, T-odd nuclear forces. The P, T-odd nuclear forces induce larger nuclear moments 
than a single nucleon EDM (Section IXC 3). This is what makes atomic experiments so 
competitive compared to neutron experiments in probing CP- violation in the hadron sector. 
As we will see in Section X, atomic experiments are more sensitive than neutron experiments 
to many underlying CP-violating mechanisms. 

1. The P,T-odd nucleon-nucleon interaction 

The P, T-odd nucleon-nucleon interaction is the dominating nuclear mechanism inducing 
atomic EDMs in diamagnetic atoms and molecules. 

The P, T-odd nucleon-nucleon interaction, to first-order in the velocities p/m, can be 
presented as [165] 
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G 1 

W ab = ^/|2^ ^ ab<Ta ~ VbaVb) ■ V a 5(r a - r 6 ) + i]' ab [cr a x cr b ] {(p a - p b ), 5{r a - r b )}) , 

(171) 

where { , } is an anticommutator, G is the Fermi constant of the weak interaction, m is 
the nucleon mass, and <x, r, and p are the spins, coordinates, and momenta of the nucleons 
a and b. The dimensionless constants 77^ and r)' ab characterize the strength of the P, T-odd 
nuclear interaction (experiments on EDMs are aimed to measure these constants). 

If we consider the P, T-odd interaction between a single unpaired nucleon and a heavy 
spherical core, then we can average the two-particle interaction (171) over the core nucleons 
to obtain the effective single-particle P, T-odd interaction between the nucleon and core 
[165], 

Here it has been assumed that the proton and neutron densities are proportional to the total 
nuclear density Pa{ v )\ the dimensionless constant 

Va = -jVap + -jVan ■ (173) 

Notice that there is only one surviving term from the P, T-odd nucleon-nucleon interaction 
(171); this is because all other terms contain the spin of the internal nucleons for which 
(a) = 0. 

The shape of the nuclear density pa and the strong potential U are known to be similar; 
we therefore take 

PA(r) = P j^U(r) . (174) 
Then Eq. (172) can be rewritten in the following form: 

W = t*.VU, ^V ,^ 7^ = -2 x 10- 21 ^ cm . (175) 

2v2m U {[)) 

Now it is easy to find the solution of the Schrodinger equation including the interaction W 
[165], 

(H + W)^ = EiP , (176) 
^ = (l + ^-V)^, (177) 

where ip is the unperturbed solution (Hip = Eip). The density arising from the wave function 
(177) is 

p = -0t,0 = ^ + £ V • (V>W) (178) 

The second term is the P, T-odd part of the density which generates the nuclear P, T-odd 
moments. 
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The electric dipole (149), Scruff (158), and magnetic quadrupole (167,168) nuclear mo- 
ments induced by the P, T-odd nucleon-nucleon interaction, through the perturbed density 
(178), are [165] 

d=-e^(q-j)t I , (179) 

5 = -» + 7TlH-H' (180) 

M q)(2I - l)t z , (181) 

m 

where q = (1) for an external neutron (proton), r 2 and r 2 x are the mean-square radii of 
the nuclear charge and external nucleon, respectively, and 

1 = 1-1/2 ' (182) 

I is the orbital angular momentum of the external nucleon. The recoil effect for the electric 
moments (the motion of the nuclear core around the center-of-mass; see Ref. [165]) has been 
taken into account (q — > q — Z/A in the expression for d). In the single-particle model, the 
recoil effect for the Schiff moment disappears due to the cancellation of its contributions to 
the first and second (screening) terms in Eq. (148). 

Note that for nuclei with an unpaired nucleon in the state S1/2 (in the simple shell model), 
such as 203 > 205 T1, the Schiff moment is reduced to the difference of two approximately equal 
terms, 

S(s l/2 ) oc (r e 2 x - r 2 ) . (183) 

In obtaining numerical values for the Schiff moment in an analytical calculation it is usually 
assumed that r 2 x = r 2 = (3/5)R%, where Rn = tqA 1 ^ 3 , t-q = 1.1 fm. Then the Schiff moment 
Eq. (183) vanishes. This cancellation makes calculations for 203 > 205 T1 unstable. 

The moments we have discussed so far are produced by a valence nucleon. In the work 
[279] it was shown that core nucleons make a contribution to P, T-odd moments that is 
comparable to that of a valence nucleon. 

Of particular interest is 199 Hg, which gives the best limit on the nuclear Schiff moment. 
In the 199 Hg nucleus the unpaired nucleon is a neutron. It doesn't contribute to the Schiff 
moment directly (see Eq. (180)). The nuclear Schiff moment arises due to the polarization of 
the protons of the core by the P, T-odd field of the external neutron. (The charge distribution 
must be distorted to give a P, T-odd correction to the charge density.) The strength of the 
P, T-odd interaction is defined by the parameter f] np . A numerical calculation of the Schiff 
moment for 199 Hg in the Woods-Saxon potential with spin-orbit interaction gives [280,279] 

S( 199 Hg) = -1.4 x 10~V P e fm 3 . (184) 

An analytical treatment of the electron relativistic corrections to the Schiff moment (that 
is, calculation of the local dipole moment) of 199 Hg shows that these corrections are small 
- they reduce the Schiff moment S only by about 25%. A many-body treatment of the 
199 Hg Schiff moment has been performed recently in the work [281]. A finite-range P, T-odd 
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nucleon-nucleon interaction was used and core polarization was calculated in the random- 
phase approximation. The result of Ref. [281] is 

S( 199 Hg) = -0.0004 gg - 0.055 gg x + 0.009 gg 2 e fm 3 , (185) 

where g = g w NN is the strong pion-nucleon coupling constant, g = g n NN are the P,T- 
violating isoscalar (i — 0), isovector (i = 1), and isotensor (i = 2) pion-nucleon couplings; 
see, e.g., Refs. [251,282] for the form of the finite-range P, T- violating interaction. Eqs. 
(184,185) can be compared by using the relation [281] i] np ~ (Gml / \/2)~ 1 g (g + g 1 — 2g 2 ) ~ 
7x 10 6 g(g + gi-2g 2 ), so Eq. (184) gives S( 199 Hg) ~ -0.09 g(g + gi-2g 2 ) efm 3 . It is seen 
that while the contribution of the isovector channel does not change much from the value in 
Eq. (184), the isoscalar channel is suppressed by two orders of magnitude and the isotensor 
channel by one order of magnitude. These corrections are due largely to the inclusion of 
core polarization [281]. 

As we mentioned earlier, of the EDM experiments performed so far, only the 133 Cs 
measurement can be interpreted in terms of a magnetic quadrupole moment of the nucleus. 
A calculation in the Woods-Saxon potential with the spin-orbit interaction gives [165] 

M( 133 Cs) = 1.7 x 10"V — fm (186) 

m 

Single-particle calculations have been performed in [283] for the external nucleon contribu- 
tion and in [284] for the core contribution. These calculations show that the MQM is very 
sensitive to the shape of the P, T-odd potential relative to the shape of the central field 
potential and to the spin-orbit potential. It was found that the core contributions arising 
from the interaction proportional to r] is comparable to that of the valence contribution, and 
that if i] ~ t]' then the core contribution is several times larger than the valence contribution. 

Results of calculations of nuclear P, T-odd moments of current interest are presented in 
Table XVI. 



2. The external nucleon EDM 

Even though the nucleon-nucleon interaction may be more effective in inducing nuclear 
P, T-odd moments and hence atomic EDMs, measurements of nucleon EDMs are interesting 
in their own right. Also, nucleon EDMs measured from atomic experiments may be compared 
to those from direct measurements. Currently, the limit on the neutron EDM from the Hg 
measurement (124) is competitive with those from direct neutron EDM searches (125); see 
Tables XIX,XX. Here we will merely quote the results of the work [250] for the nuclear 
Schiff and magnetic quadrupole moments induced by a single unpaired nucleon (neutron or 
proton) with an intrinsic EDM <i n , p , 



r 2 

cx 



<?-!,* 

10 n ' P L~V i + 

M =^(2J-l)tj. (188) 
m 

See, e.g., [250,13] for details. Notice the similarity between these expressions and those for 
the corresponding moments induced by P, T-odd nuclear forces (180), (181). In the approxi- 
mations used, there is a simple correspondence between the P, T- violating parameters. For 
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the Schiff moment, — eqC, <-> d n>p , while for the MQM, — e£(g — //)<-> d n>p . The second con- 
tribution (oc /i) to the MQM from the P, T-odd nuclear forces appears from the spin term 
in (103). This has no analog in the case of d n>p since a stationary EDM cannot induce a 
MQM. It is the orbital motion of a nucleon with an intrinsic EDM that induces the nuclear 
MQM [250]. 

In the simple shell model the unpaired neutron in the 199 Hg nucleus carries the nuclear 
spin J, and the Schiff moment is induced by the EDM of just this neutron. In this picture, the 
induced Schiff moment can be calculated using Eq. (187) and the simplifying assumptions 
r e 2 x = r\ = (3/5)R%, R N « 1.1A 1/3 fm, giving ,S( 199 Hg) « 2.2 d n fm 2 . Proton EDMs 
also contribute due to configuration mixing. It is possible to estimate their contribution 
by comparing the experimental value of the magnetic moment of 199 Hg with that of the 
simple shell model; see Ref. [257]. In this way, it is found [285] that d n can be replaced by 
(d n + O.ldp), 

,S( 199 Hg) = (2.2 d n + 0.2 d p ) fm 2 . (189) 

A numerical calculation of the nuclear Schiff moment of 199 Hg induced by neutron and 
proton EDMs has recently been performed, with core polarization accounted for in the RPA 
approximation [286]. The result is [286] 

,S( 199 Hg) = (1.9 d n + 0.2 d p ) fm 2 . (190) 

(See [286] for the discussion of uncertainty.) 



3. Comparison of the size of nuclear moments induced by the nucleon-nucleon interaction and the 

nucleon EDM 

Here we consider the enhancement of the nuclear EDM induced by P, T-odd nuclear 
forces compared to that induced by an external valence nucleon using the simple one-boson 
exchange model, following [280,279,287]. The largest contribution to the constant 77 is prob- 
ably given by the lightest 7r°-meson, 

G ^ g« NN ? NN 

9-kNN and g n NN are the constants of the strong and T-odd n meson- nucleon interactions, 

{ig-nNNU^n + g° NN pp)7r° + \/2{ig^ NN p^ b n + g~ NN pn){n') ] + ... . (192) 
A neutron EDM is induced through virtual creation of a ir~ meson [288], 

m 4tk z m n 

Here M ~ m p ~ 700 MeV is the scale at which the 7r-meson loop converges. The values of rj 
and d n are expressed in terms of different quantities, g® NN and g~ NN , respectively. However, 
for example, in the model of T- violation with the #-term, \g nNN g® NN \ = \g nNN g~ NN \ = 
O.37|0| [288]. Taking \g vNN g° NN \ ~ \g vN Ng~ NN l it is found that [165,280,279,287] 
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„ 27r(m2r 3 |f/(0)|) _1 ~ 40 , (194) 

that is, the nuclear EDM exceeds the nucleon EDM by one to two orders of magnitude. 
Similarly, P, T-odd nuclear forces generate all P, T-odd nuclear moments, such as Schiff and 
MQM moments, 10 — 100 times larger than those generated by the presence of a nucleon 
EDM [287]. 



D. Nuclear enhancement mechanisms 

So far we have considered P, T-odd nuclear moments in spherical nuclei. However, in 
non-spherical nuclei there is the possibility of enhancement due to (i) the presence of a 
low-lying level with opposite parity and the same angular momentum with respect to the 
ground state; and (ii) collective effects. 



1. Close-level enhancement 

It is known that nuclei with non-spherical symmetry have close levels of opposite parity. 
It was pointed out in Ref. [289] that due to the existence of a close level of opposite parity 
with the same angular momentum as the ground state, the nuclear EDM can be enhanced; 
calculations of enhanced EDMs and MQMs were performed in [251], Schiff moments in [165]. 
In the frozen frame, the contribution to the ^-components of the electric dipole, Schiff, and 
magnetic quadrupole nuclear moments in the ground state fl due to the close opposite parity 
state f2 is 

rp = 2 (m PT \n)(n\f\n) ^ 

En — Eq 

where T = d z , S z , M zz . The magnetic quadrupole moment in heavy stable nuclei can be 
enhanced by an order of magnitude due to the "close level" mechanism, while for the electric 
dipole and Schiff moments this enhancement hardly exceeds ~ 5 — 10 [165]. However, this 
"close-level" enhancement is not regular: it gives contributions to P, T-odd moments with 
different magnitudes and signs even in the "nearest" nuclei, the results are unstable. 18 



2. Collective enhancement 

While the "close-level" mechanism enhances the contribution of the external nucleon to 
the P, T-odd nuclear moments, there can also be a "collective" enhancement of the nuclear 
moments that occurs due to the contribution of many nucleons. In the work [278] it was 



18 The reason is explained in [165]: taking Eq. (175) for the P, T-odd interaction Hpt, it is seen 
that W = (,cr ■ 'VU oc [H,cr ■ p], where H is the single-particle Hamiltonian, has small matrix 
elements between close levels, (fJ| [H, a ■ p]\Cl) oc Eq — E n . 
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shown that a collective magnetic quadrupole moment can be produced in deformed nuclei 
by P, T-odd nuclear forces. Unlike the close-level enhancement, this collective enhancement 
is regular: in deformed nuclei about A 2 / 3 nucleons belong to open shells that contribute to 
the MQM. The P, T-odd nuclear forces create a spin hedgehog [Eq. (170)] as in the case of 
spherical nuclei (Section IX B 1), however in the deformed case there is a non-zero magnetic 
field associated with it. The MQM of a deformed nucleus (in the rotating frame) can 
basically be calculated as a summation of the single-particle MQMs (181) of all nucleons 
in the open shells. Notice, from Eqs. (181,182), that spin-orbit pairs 1 = 1 + 1/2 and 
1 = 1 — 1/2 make contributions to the collective MQM of opposite sign. A sufficiently large 
spin-orbit splitting is therefore required to avoid cancellation, and this is satisfied in nuclei. 
This collective mechanism gives an order of magnitude enhancement of the nuclear MQM 
in deformed nuclei compared to spherical nuclei. 



We will now move on to the collective P, T-odd nuclear moments produced by P, T-odd 
nuclear forces that arise in nuclei with static octupole deformation [290,274,275]. There is an 
enhancement of these collective moments, compared to single-particle nuclear moments, due 
to the collective nature of the intrinsic moments and the small energy separation between 
members of parity doublets. This enhancement can be as large as 1000 times. 

Static octupole deformation in the ground state has been demonstrated to exist in nuclei 
in the regions Ra-Th and Ba-Sm. It produces effects such as parity doublets, large dipole 
and octupole moments in the intrinsic frame of reference and enhanced electric dipole and 
octupole transitions; see the review Ref. [291]. 

While it has been shown that the Schiff and electric dipole and octupole moments are 
enhanced in nuclei with octupole deformation, we will focus our attention on the nuclear 
Schiff moment (the EDM is not of direct interest, in atoms it is screened by atomic electrons; 
also, the atomic EDM induced by the electric octupole moment is small since it does not 
mix s and p electron orbitals.) 

The mechanism generating collective P, T-odd moments is the following. In the "frozen" 
body frame collective moments can exist without any P, T- violation. However, the nucleus 
rotates, and this makes the expectation value of these moments vanish in the laboratory 
frame if there is no P, T- violation. (For example, the intrinsic Schiff moment is directed along 
the nuclear axis, Si ntr = Si n tr n j and in the laboratory frame the only possible correlation 
(n) oc I violates parity and time reversal invariance.) The P, T-odd nuclear forces mix 
rotational states of opposite parity and create an average orientation of the nuclear axis n 
along the nuclear spin I, 



3. Octupole deformation; collective Schiff moments 



(n z ) = 2a 



KM 



(196) 



1(1 + 1) 



where 



a = 



(i>-\w\^ + ) 



(197) 
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is the mixing coefficient of the opposite parity states, K — [I • n| is the absolute value of 
the projection of the nuclear spin I on the nuclear axis, M — I z , and W is the effective 
single-particle potential (172). The Schiff moment in the laboratory frame is 

S z = S intr (n z ) = Sintr-^^ • (198) 

In the "frozen" body frame the surface of an axially symmetric deformed nucleus is 
described by the following expression 

R(6) = R N (l + J2m (0)) ■ (199) 
i=i 

To keep the center-of-mass at r = we have to fix f3\ [292]: 

A = -3#E + ■ (200) 

V 4tt ^ v /( 2 / + l)(2/ + 3) 

Assuming that the distributions of the protons and neutrons are the same, the electric dipole 
moment e(r) = (since the center-of-mass of the charge distribution coincides with the 
center-of-mass) and hence there is no screening contribution to the Schiff moment. We also 
assume constant density for R < R{6). The intrinsic Schiff moment Si ntr is then [290,274] 

2°f h \/(2( + l)(2i + 3) N 2(W35' 

where the major contribution comes from P2P3, the product of the quadrupole /3 2 ~ 0.1 and 
octupole /9 3 ~ 0.1 deformations. The estimate of the Schiff moment in the laboratory frame 
gives [274] 

eV 

S ~ a,S intr ~ 0.05 e[3 2 $ZA 2/3 rirl- — ~ 700 x 10^ 8 r]dm 3 , (202) 

hi+ — hi— 

where r ~ 1.2 fm is the internucleon distance, E + — ~ 50 keV. This estimate (202) is 
about 500 times larger than the Schiff moment of a spherical nucleus like Hg (see Eq. (184)). 

See Table XVII for calculations [274] of Schiff moments in nuclei assuming static octupole 
deformation. An attractive candidate for P, T-odd studies is radium, and recently several 
laboratories around the world have considered performing EDM experiments with it. As 
well as the possibility for a large Schiff moment, the atomic EDM is large due to high Z. 
And if measurements can be performed for metastable atomic states, further enhancement 
can occur due to the presence of close opposite parity levels [162,163]. A Woods-Saxon 
calculation for the Schiff moment of 225 Ra gives [274] 

S( 225 Ra) = 300 x 10~ 8 r] np e fm 3 . (203) 

Recently, a self-consistent calculation of the nuclear Schiff moment of 225 Ra was performed, 
with core polarization taken into account [293], 

,S( 225 Ra) = -5.06 gg + 10.4 gg x - 10.1 gg 2 e fm 3 . (204) 
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This calculation was carried out in the approximation of the zero-range P, T-odd interaction. 
It was found that the Schiff moment in the rotating frame is up to twice as large as the value 
calculated in [274] (see Table XVII). However, the Schiff moment in the laboratory frame 
was found to be suppressed by between 1.5 and 3 times due to suppression of the matrix 
element of the P, T-odd interaction. [Comparison between Eqs. (203,204) can be made using 
the relation following Eq. (185); Eq. (203) then gives S( 225 Ra) ~ 20 g(g + g~i- 2g 2 ) e fm 3 .] 
Improved calculations (taking into account the finite-range of the P, T-odd interaction) are 
in progress [293]. It is seen by comparison with the calculations for 199 Hg (184,185) that the 
radium Schiff moment is several hundred times larger. 

Note that S in Eq. (202) is proportional to the squared octupole deformation parameter 
In Ref. [294] it was pointed out that in nuclei with a soft octupole vibration mode 
~ (0.1) 2 , i.e., the Schiff moments induced in nuclei with a soft octupole vibration mode 
are of the same magnitude as those induced in nuclei with static octupole deformation. In 
a recent work [295] Schiff moments of nuclei with soft octupole vibrations were calculated 
and found to have a similar enhancement as in the static case. This means that a number 
of heavy nuclei can have large collective Schiff moments. 

The effect of static octupole deformation on the size of the single-particle MQM was 
considered in [275]. It was found that generally there is no significant enhancement due to 
this mechanism. 

There are several experiments in preparation aimed to detect EDMs of heavy atoms with 
deformed nuclei; see Section X. 

E. Calculations of atomic EDMs induced by P, T- violating nuclear moments; 
intepretation of the Hg measurement in terms of hadronic parameters 

Results of atomic calculations of EDMs (of current interest) induced by nuclear Schiff 
moments are presented in Table XVIII. Of particular interest is the calculation for 199 Hg. 
A recent calculation for the atomic EDM induced by the Schiff moment yielded the result 

[285] 

d( 199 Hg) = -2.8 x 10" 17 ( g(19 f ^ 3 S) ) e cm . (205) 

This value was obtained using the new finite-size form for the Schiff potential (161) and is 
the average of two calculations: one performed in the potential V N with core polarization 
taken into account using the TDHF method and the other performed in the V N ~ 2 potential 
using the combined MBPT+CI method. The error of the result (205) is about 20%. The 
previous value for <i( 199 Hg) induced by the nuclear Schiff moment, <i( 199 Hg) = —4 (in the 
same units as Eq. (205)), was estimated in Refs. [280,279] from an atomic calculation [256] 
of the EDM induced by the tensor electron-nucleon interaction. 
From (124,205), the best limit on the Schiff moment follows, 

S( 199 Hg) = (3.8 ± 1.8 ± 1.4) x 10~ 12 e fm 3 . (206) 

As we have discussed, the Schiff moment can be induced from a number of P, T- violating 
mechanisms: due to a P, T- violating nucleon-nucleon interaction or due to a static EDM of 
an unpaired nucleon. The limit on r) np from Eqs. (206,184) is 
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^ = -(2.7 ±1.3 ±1.0) x 10 



,-4 



(207) 



The limits on neutron and proton EDMs are [Eqs. (206,189)] 



d n = (1.7 ± 0.8 ± 0.6) x 10~ 25 e cm 
d p = (17±8±6) x 10~ 25 e cm . 



(208) 
(209) 



Using instead the recent calculation (190) will change the limits on d n and d p only slightly. 
The result (185) suggests that the limits on the isoscalar and isotensor P, T-odd couplings 
may be substantially weaker than those that would follow from Eq. (207) using Eq. (184). 

X. CURRENT LIMITS ON FUNDAMENTAL P, T- VIOLATING PARAMETERS 
AND PROSPECTS FOR IMPROVEMENT 



EDM measurements have already excluded several models of CP-violation and the pa- 
rameter space of currently popular models is strongly constrained. In this review we will not 
discuss the sensitivities of the various EDM measurements to different CP- violating models; 
for such an analysis, see, e.g., the reviews [296,253,260] and the book [13]. The problem is 
that in new theories there are many free parameters, and this gives a whole range of possible 
values for CP- violating effects. Here we compare the sensitivities of different measurements 
to fundamental CP- violating interactions by placing limits on phenomeno logical parameters. 

In Tables XIX,XX we present the best limits on fundamental P, T- violating parameters 
extracted from EDM measurements in atoms, molecules, and neutrons. In previous sections 
we presented limits on hadronic and semi-leptonic CP- violating parameters at the nucleon 
scale [nucleon EDMs, Eqs. (125,208,209); nucleon-nucleon interaction, Eqs. (207); electron- 
nucleon interaction, Eqs. (132,133,136)] and in the leptonic sector directly on the electron 
EDM [Eq. (143)]. These are summarized in Table XIX. Here we constrain CP- violation at 
the quark scale from limits at the nucleon scale. Fig. 22 shows which of these parameters 
are related. 

We will begin with a discussion of limits on CP-violating parameters in the hadronic 
sector. Currently, the best limits on the neutron EDM come from direct neutron measure- 
ments (125), and the limit from mercury is not far behind (208). The best proton EDM limit 
comes from mercury (209). Constraints on nucleon EDMs constrain the following param- 
eters at the quark level: quark EDMs, chromoelectric dipole moments (CEDMs) [CEDMs 
are analogous to quark EDMs, though with the external electromagnetic field replaced by a 
gluonic field] , P, T-odd quark-quark four-fermion interactions, the QCD phase defining the 
strength of the term GG, and another parameter defining the strength of the term GGG [in 
all models of CP violation considered this parameter is not the dominating mechanism in- 
ducing neutron EDMs [302], so we will not discuss it further]. Limits on nucleon EDMs can 
also constrain P, T-odd lepton-quark four-fermion interactions which induce quark EDMs 
at the one-loop level [303]; however, the limits involving light quarks are weaker than those 
obtained at tree-level from atomic measurements. 



A. Summary of limits 
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We use the following calculations to relate the nucleon and quark parameters, with 
the corresponding limits presented in Table XX. The limit on the QCD 6 term is ar- 
rived at using the relation obtained in Ref. [301], d n = 1.2 x 10~ 16 (9e cm. To re- 
late d n to the quark CEDMs d q and quark EDMs d q we use the recent calculation 
[300], d n = (1 ± 0.5)[0.55e(J d + 0.5d u ) + 0.7(d d - 0.25«]. For the P, T-odd quark- 
quark interactions we use the phenomenological parameters [297,298,13] defined by the 
Hamiltonian h qq = ^[k s {qii^qi){q 2 q 2 ) + ^(^Ts^iX^^) + k t ±e^ ua(3 (ua ^u){da a pd) + 
kt^^ua/3(ua^ v t a u)(da a /3t a d)}, where t a are the 577(3) generators and the quarks q±, q 2 = u, d. 

From the limit on the P, T-odd nucleon-nucleon interaction 77 from mercury follows 
the best limit on the P, T-odd pion-nucleon coupling g^NN] they are related through 
V — 9nNNSInNNV^/(Gml), where the strong coupling g^NN = 13.6. The relation \g n NN\ ~ 
0.027|6>| (Ref. [288]) is used to place a constraint on the QCD 9 term, tighter than those 
obtained from direct neutron measurements. Limits on CEDMs of quarks can be obtained 
from the limit on g n NN- We use the calculation of Ref. [299], g w NN = 2(d u — dd) / (10~ 14 cm) . 
Again, this limit is better than those from the neutron. Constraints on the phenomenological 
P, T-odd quark-quark interactions defined above are obtained from the constraints on g^NN, 
with the exception of k° s with qi 7^ q 2 and k\ for which the P, T-odd neutral pion-nucleon 
vertex is insensitive [13]. 

Limits on the P, T-odd electron-nucleon parameters can be broken down into limits 
on phenomenological P, T-odd electron-quark parameters ki q , k 2q , k 3q defined according to 
K q = ^=[k lq qqei^e + k 2q \e KX ^qa KX qea^ v e + k 3q qi-f 5 qee\; see Ref. [13]. 

Until this point we have not mentioned the P, T-odd electron-electron interaction, h ee = 
-^=/cf (er/5e)(ee) . This interaction induces atomic EDMs smaller than those induced by the 
electron-nucleon interaction. It does not benefit from two enhancement factors present in 
the latter interaction: the relativistic factor arising from the outer electrons being in the 
nuclear vicinity; and the interaction with A nucleons (rather than two i^-shell electrons in 
the former interaction); see Ref. [13]. 

The limit on the electron EDM can be reduced to limits on the P, T-odd electron-quark 
interaction [303] and on the electron-electron interaction [13] that are competitive with those 
obtained at tree level. However, we do not include these into our table. 

In relating the P, T-odd parameters at the nucleon and quark levels we used the most 
recent calculations available. For references to other calculations, please see, e.g., the review 
[253] and book [13]. 

B. Ongoing/future EDM experiments in atoms, solids, and diatomic molecules 

There is a new generation of experiments in preparation aimed to measure EDMs of 
ions in solids (gadolinium gallium garnet and gadolinium iron garnet) [228,229]. Here it is 
expected that the sensitivity to T-violating effects (in particular, the electron EDM) will be 
improved by several orders of magnitude. The solid-state nature of the problem complicates 
the calculations required for interpretation of the measurements. A series of calculations 
have been performed in [304]. 

One of the primary goals of the TRRiP facility under construction at Groningen is 
to measure permanent EDMs of radioactive atoms and ions, in particular Ra (see, e.g., 
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[140]). Groups at Yale, Argonne, and Los Alamos are also considering performing EDM 
experiments with Ra and Rn. These atoms can have very large EDMs due to their high 
Z and the presence of nuclear (static/vibrational) octupole deformation (Section IX D 3). 
Also, metastable atomic states of Ra have close levels of opposite parity, and this can be 
exploited to obtain an enhanced EDM effect [162,163]. 

EDM experiments with both paramagnetic and diamagnetic diatomic molecules are un- 
derway. The first EDM experiment with paramagnetic molecules (YbF) was performed 
recently at the University of Sussex [226], and while the result gave a looser bound on the 
electron EDM than the Tl experiment [11], a substantial improvement in the result is ex- 
pected. Experiments with PbO excited to the metastable a(l) state have begun at Yale 
[270]. 

XI. CONCLUDING REMARKS 

Exciting developments in violations of fundamental symmetries are expected in the next 
few years. Improved EDM measurements underway, including the new generation of exper- 
iments in solids, diatomic molecules, and radioactive atoms, are expected to yield limits on 
electric dipole moments that are several orders of magnitude better than the current ones. 
Or perhaps an EDM will be unambiguously detected? Popular models such as supersymme- 
try will be put to the test. Improved precision tests of parity violation in atoms in a single 
isotope and in a chain of isotopes will provide crucial tests of physics beyond the standard 
model complementary to each other and to other electroweak tests. New measurements of 
the nuclear anapole moment are anticipated, and they will have important consequences for 
the theory of parity violating nuclear forces. 
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TABLES 



TABLE I. Results of atomic PNC experiments measured to better than 5%. Results of optical 
rotation experiments are given in terms of Im(EpNc /Ml); Stark-PNC experiments are given in 
terms of Im(EpNc/P)- 



Atom 


Transit inn 


Group 


Year 


Ref. 


Measurement 












lm{E PNC /Ml) 


lm(E PNC /(3) 












[io- 8 ] 


[mV/cm] 


209 Bi 




Oxford 


1991 


[32] 


-10.12(20) 




208p b 


3 Po - 3 Pi 


Seattle 


1993 


[34] 


-9.86(12) 








Oxford 


1996 


[35] 


-9.80(33) 




205 rQ 


6P 1/2 - 6P 3/2 


Oxford 


1995 


[37] 


-15.68(45) 








Seattle 


1995 


[38] 


-14.68(17) 




133 Cs 




Boulder 


1988 


[44] 




-1.576(34) 






Boulder 


1997 


[8] 




-1.5935(56) 



TABLE II. Most precise calculations of PNC amplitudes Epjyc for atoms and transitions listed 
in Table I. Units: W' u iea B (-Q w /N). 



Atom Transition Ep^c 3, Ref. 

^ 4 S 3/2 " 2 D 3/2 28(3) WW] 

208 Pb 3 Po- 3 Pi 30(2) [56] 

205 T1 6P 1/2 -6P 3/2 27.0(8) [54] 

27.2(7) [55] 

133 Cs 6S-7S 0.904(5) [9] 



a The values for the PNC amplitudes for Cs [9] and for Tl 6Pi/ 2 — 6P 3 / 2 [55] include corrections 
beyond the other calculations. In particular, for Cs the contributions of the Breit interaction 
and vacuum polarization due to the strong nuclear Coulomb field are included. For Tl the Breit 
interaction is also included. The remaining corrections for Cs and Tl are discussed in detail in 
Sections V,VI, respectively. These corrections would be inside the error bars for the other atoms 
and transitions in the table. 



TABLE III. Summary of experimental results for PNC in cesium 6S — 7S, — Im(EpNc) / P', 
units: mV/cm. 



Group 



Year 



Ref. 



Value 



Paris 

Boulder 

Boulder 

Boulder 

Paris 



1982,1984 

1985 

1988 

1997 

2003 



[39-41] 

[43] 

[44] 

[8] 
[42] 



1.52(18) 

1.65(13) 

1.576(34) 

1.5935(56) 

1.752(147) 
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TABLE IV. Summary of calculations of the PNC El amplitude for the cesium 65 — 75 tran- 
sition; units are 10 _11 ieas(— Qw/N). 



Authors 


Year 




Ref. 


Value 


Bouchiat, Bouchiat 


1974, 


1975 


[14,22] a 


1.33 


Loving, Sandars 


1975 




[59] a 


1.15 


Neuffer, Commins 


1977 




[60] a 


1.00 


Kuchiev, Sheinerman, Yahontov 


1981 




[61] a 


0.75 


Das 


1981 




[62] b 


1.06 


Bouchiat, Piketty, Pignon 


1983 




[63] a 


0.97(10) 


Dzuba, Flambaum, Silvestrov, Sushkov 


1984, 


1985 


[64,65] b 


0.88(3) 


Schafer, Miiller, Greiner, Johnson 


1984 




[66] b 


0.74 


Martensson-Pendrill 


1985 




[67] b 


0.886 


Plummer, Grant 


1985 




[68] b 


0.64 


Schafer, Miiller, Greiner 


1985 




[69] b 


0.92 


Johnson, Guo, Idrees, Sapirstein 


1985, 


1986 


[70] c 


0.754,0.876,0.856 


Johnson, Guo, Idrees, Sapirstein 


1985, 


1986 


[70] b 


0.890 


Bouchiat, Piketty 


1986 




[71] a 


0.935(20)(30) 


Dzuba, Flambaum, Silvestrov, Sushkov 


1987 




[54] b 


0.90(2) 


Johnson, Blundell, Liu, Sapirstein 


1988 




[72] b 


0.95(5) 


Parpia, Perger, Das 


1988 




[73] c 


0.879 


l^ZLLUd, JJ IdlllUdjUlIl, u U511KUV 






L° 'J 


u.yuoi v 1 


Hartley, Sandars 


1990 




[74] c 


0.904(18) 


Hartley, Lindroth, Martensson-Pendrill 


1990 




[75] b 


0.933(37) 


Blundell, Johnson, Sapirstein 


1990, 


1992 


[16] b 


0.905(9) 


Safronova, Johnson 


2000 




[76] b 


0.909(11) 


Kozlov, Porsev, Tupitsyn 


2001 




[58] bd 


0.901(9) 


Dzuba, Flambaum, Ginges 


2002 




[9] bd 


0.904(5) 



a Semi-empirical calculations. 

b Ab initio many-body calculations. 

c Combined many-body and semi-empirical calculations. 

d The difference between the values of [58,9] and previous ones is due to the inclusion of the Breit 
interaction in [58] and the Breit and strong field vacuum polarization in [9]. 



TABLE V. Removal energies for Cs in units cm . 

State RHF £W t Experiment a 

65 27954 32415 31492 31407 

75 12112 13070 12893 12871 

6P 1/2 18790 20539 20280 20228 

7P 1/2 9223 9731 9663 9641 



a Taken from [82]. 
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TABLE VI. Radial integrals of El transition amplitudes for Cs in different approximations. 
The experimental values are listed in the last column, (a.u.) 



TV a rm i t i n n 


RHF 


TDHF 


with fitting 




with fitting 


FiYnpri in pnf" 


QS - 6P 1/2 


6.464 


6.093 


5.499 


5.497 


5.509 


5.5232(91) a , 5.4979(80) b , 












5.5192(58) c , 5.512(2) d , 














5.524(5)(l) c 


7S - 6P 1/2 


5.405 


5.450 


5.198 


5.190 


5.204 


5.185(27) f 


7S-7P 1/2 


13.483 


13.376 


12.602 


12.601 


12.612 


12.625(18)8 



a Ref. [92]. 
b Ref. [93]. 

c Ref. [94] . Deduced from the van der Waals coefficient C$ . 
d Ref. [95]. Deduced from photoassociation spectroscopy. 

e Ref. [96]. Deduced from their measurement of the static dipole polarizability. 
f Ref. [97]. 
SRef. [98]. 



TABLE VII. Calculations of the hyperfine structure of Cs in different approximations. In the 
last column the experimental values are listed. Units: MHz. 



State 


RHF 


TDHF 


E( 2 ) 
with fitting 


t 


t 

with fitting 


Experiment 


6S 


1425.0 


1717.5 


2306.9 


2315.0 


2300.3 


2298.2 a 


7S 


391.6 


471.1 


544.4 


545.3 


543.8 


545.90(9) b 


6P1/2 


160.9 


200.3 


291.5 


293.6 


290.5 


291.89(8) c 


7Pl/2 


57.6 


71.2 


94.3 


94.8 


94.1 


94.35 a 



a Ref. [100]. 
b Ref. [101]. 
c Ref. [102]. 
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TABLE VIII. Contributions to the 6S — 7S Ep^c amplitude for Cs in units 
10~ n ieaB(— Qw/N). (S corresponds to the (unfitted) "dressed" self-energy operator.) 



TDHF 0.8898 
Brueckner-type correlations 

(4> 7s \± s (e 7s )\5X 6s ) 0.0773 

(d^ 7s \± p (e 7s )\X 6s ) 0.1799 

(8Y 7s \± s (e 6s )\^ 6s ) -0.0810 

(y 7s |S p (e6 S )|# 6s ) -0.1369 

Nonlinear in E correction -0.0214 

Weak correlation potential 0.0038 

Structural radiation 0.0029 

Normalization -0.0066 

Subtotal 0.9078 

Breit -0.0055 

Neutron distribution correction -0.0018 
QED radiative corrections 

Vacuum polarization (Uehling) 0.0036 

Self-energy and vertex -0.0072 

Total 0.8969 



TABLE IX. Values for Ep^c m different approximations; units 10 n iea,B(— Qw/N). 





S( 2 ) with fitting 


t 


E with fitting 


Epnc 


0.901 


0.904 


0.903 
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TABLE X. Calculations of the anapole moment n a for 133 Cs. In the last column the value 
for the anapole moment constant n a is presented, using DDH best values of the meson-nucleon 



couplings (see Table XII). 

Ref. K a x 10 2 K a {DDH) 

Single- [21] a 8^ 0.36 

particle [21] b 6g p = (34/ w - 9h° p - 2h p + l/i 2 - 5/i° - 5/£) x 10 2 0.27 

[174,170,171] c 4.9^ + 0.65^ 0.26 

[168,19] d 5.5# p 0.25 

Many- [169] c 23 x lO 2 /^ 0.11 

body [170,171] f 4Ag p + 0A5g n + 0A5g pn - 0m 9np 0.23 

[171]s 2.9<7 p + 0.18<7 n + 0.36 5pri -0.02 5np 0.15 

[172] h (26.98/^ - 7.01/i° - 1.72/i* + 0.16/i 2 - 4.48/i° - 2. 16/£) x 10 2 0.21 



a Formulae (107) and (115). 

b Woods-S axon potential with spin-orbit interaction, formulae (108) and (115). 
c Woods-Saxon potential with spin-orbit interaction; the contact and spin-orbit currents are in- 
cluded. This calculation includes some contributions beyond the single-particle approximation. 
d Harmonic oscillator potential with spin-orbit interaction. Configuration mixing is taken into 
account semi-empirically. 

e Large-basis shell-model calculations; just the pion contribution was calculated. 
f Many-body effects taken into account in RPA. 

g Many-body effects taken into account in RPA. More complete treatment than Ref. [170]. 
h Large-basis shell-model calculations. 



TABLE XL Nuclear calculations of K2 and atomic calculations of hq for 133 Cs and 203,205 T1 



Ref. 


133 Cs 


K Q 


203,205 rpj 


«0 


Eq. (3) a 


-0.05 




-0.05 




[18] b 




0.017 




0.014 


[168,19] 


-0.038 c 


0.027 


-0.027 c 


0.022 


[172] d 


-0.063 




-0.064 




[176] e 




0.017 







a Single-particle, standard model value with sin 2 9\y = 0.23. 
b Eq. (11). 

c Nuclear configuration mixing was taken into account semi-empirically. 
d Large-basis nuclear shell model calculations. 

e Atomic many-body calculations in third-order perturbation theory. Core polarization is included 



in RPA. 



89 



TABLE XII. DDH "best values" (/„-, h) in units 10 7 and effective coupling constants (g) at 
DDH values for the meson- nucleon couplings. 



fn 


h° 


h 1 
tip 


h 2 


h° 


K 


9pn 


9np 


9pp — 9nn 


9 P 


9n 


4.6 


-11.4 


-0.19 


-9.5 


-1.9 


-i.i 


6.5 


-2.2 


1.5 


4.5 


0.2 



TABLE XIII. Calculations of the anapole moment n a for 203 ' 205 T1. In the last column the value 
for the anapole moment constant K a is presented, using DDH best values of the meson-nucleon 
couplings (see Table XII). 





Ref. 


K a X 10 2 


n a (DDH) 


Single- 


[21] a 




0.48 


particle 


[21] b 


10 5p = (56/, - 16/i° - 4hl + lh 2 p - 9/tO - 9hl) x 10 2 


0.43 




[174,170,171] c 


7.8^ p + 0.85g pn 


0.41 




[168,19] d 


6.0g p 


0.27 


Many- 


[170,171]° 


7.lg p + 0.35# n + 0Mg pn - 0.06g np 


0.36 


body 


[199] f 


5.3g p + 0Ag n 


0.24 




[171]* 


A.3g p + 0.1g n + 0.64^ pn - 0.06g np 


0.24 




[172] h 


(13.98/, - 2.92/i° - 0.26^ + 0.23/i 2 - 2.Uh® - 0.98/£) x 10 2 


0.10 



a Formulae (106) and (115). 

b Woods-Saxon potential with spin-orbit interaction. Eq. (115) is used to express the effective 
coupling constant g p in terms of parity violating weak meson-nucleon couplings. 
c Woods-Saxon potential with spin-orbit interaction; the contact and spin-orbit currents are in- 
cluded. This calculation includes some contributions beyond the single-particle approximation. 
d Harmonic oscillator potential with spin-orbit interaction. Configuration mixing was taken into 
account semi-empirically. 

e Many-body effects taken into account in RPA. 

f Large-basis shell-model calculations. Just the dominant spin-current contribution was considered. 
g Many-body effects taken into account in RPA. More complete treatment than Ref. [170]. 
h Large-basis shell-model calculations. 
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TABLE XIV. Calculations of atomic EDMs induced by P, T-odd electron-nucleon interaction. 



Atom 


G^atom / C 








(1(T 18 e cm) 


(10" 2 ° e cm) 


(10~ 23 e cm) 


Cs 


0.70 a 


0.92 b 


2.2 C 




0.71 b 








0.72 d 






Tl 


-5.1 b 


0.5 b 


1.5 C 




-(7± 2) 






Xe 


-5.6 x 1(T 5 c 


0.52 f 


1.2 e 




-4.4 x 1(T 5 s 


0.41 h 


0.95 g 






0.6 b 




Hg 


-5.9 x 10" 4 c 


2.0 f 


6.0 C 






1.3 b 





a Ref. [249], semi-empirical calculation. 

b Ref. [12,13], simple analytical calculations. 

c Using formula (135) and calculations for d a tom/C T from [13]. 

d Ref. [255], TDHF calculation with correlation corrections estimated from the size of the corrections 

for the corresponding electron EDM enhancement factors in [75]. 

e Using formulae (131), (135) and the calculation <i a tom/C T from Ref. [256]. 

f Ref. [256], TDHF calculations. 

g Using formulae (131), (135) and the calculation d at om/C T from Ref. [257]. 
h Ref. [257], TDHF calculation. 
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TABLE XV. Enhancement factors for atoms of interest. 





Atom 


Enhancement factor K 
Semi-empirical 


Ab initio 


Paramagnetic 


Rb 


24 a , 24.6 b , 16.1 c , 23.7 C , 22.0 C 


24.6 C 




Cs 


1 1 qa 1 01 d 1 00b 1 00c on oc 1 Oft f)C 1 nn AC 


1 1 A Q c 1 1 4 f 1 1 Kg 




Fr 


1150 a 


910(50) h 




Tl 


-716*, -500 c , -502 c , -607 c , -562 c 


-1041 c , -301 j , -179 f , -585 k 




Xe 3 P 2 


130 1 , 120° 




Diamagnetic 


Xe 




-0.0008 111 , -0.0008 g 




Hg 




-0.014 m , 0.0116 g 



a Ref. [213]. 
b Ref. [265]. 
c Ref. [70]. 
d Ref. [266]. 
e Ref. [214]. 
f Ref. [75]. 
SRef. [263]. 
h Ref. [267]. 
'Ref. [268]. 
JRef. [269]. 
k Ref. [261]. 
! Ref. [225]. 

111 Ref. [18]. Extracted from the many-body calculation of ci a tom(C T ) [256]; see Section VIII D 2. 
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TABLE XVI. Calculations of P, T-violating moments in spherical nuclei induced by the 
P, T-odd nucleon-nucleon interaction. 

Atom 5(e fm 3 ) x 10 8 M(^- fm) x 10 7 

T29 Xe" 1.75^/ 
199 Hg -1.4^/ 

133 Cs 3.0r/ p b 1.7r/ p b 

1.6?7p, 2.Qt] p , 1.8t] p c 

-0.3ri pp - 0.2 Vpn - 2.5r? pp + \.lr( pn A 

203,205 T1 _ 2r] b 

\.lTj pp l.Af] pn 

a Calculated using wave functions and Green's functions found in the Woods-Saxon potential with 
spin-orbit interaction, Refs. [280,279]. 

b Woods-S axon potential with spin-orbit interaction, Ref. [165]. 

c Contributions of the valence nucleon carried out in different potentials: Eq. (174), harmonic 
oscillator, Woods-Saxon, respectively. Ref. [283]. 
d Core contribution, Ref. [284]. 

c This value is comprised of a valence nucleon and core contribution. The valence proton gives 
a contribution of — 0.04r/ pp — lArj pn while the core contribution is 1.23r/ pp , in the units in the 
table. The difference between the valence contributions of the works [280,279] and [165] is due to 
a difference in the potentials used; remember that in 203 ' 205 T1 the Schiff moment is very sensitive 
to the potential, see Eq. (183). 



TABLE XVII. Nuclear Schiff moments Si n t r and S (in rotating and laboratory frames, respec- 
tively) calculated in the Woods-Saxon potential. Static octupole deformation is assumed. 





223 Ra 


225 Ra 


223 Rn 


221 Fr 


223 Fr 


5 intr (e fm 3 ) a 


24 


24 


15 


21 


20 


S(i] e fm 3 ) x 10 8a 


400 


300 


1000 


43 


500 



a Numbers are from Ref. [274]. See [274] for details. 



TABLE XVIII. Ground state EDMs of diamagnetic atoms induced by nuclear Schiff moments. 
Units: 10- 17 (5/(e fm 3 ))e cm. 



129 Xe 


223 Rn 


199 Hg 


225 Ra 


0.27 a , 0.38 b 


2.0 C , 3.3 b 


-4.0 d , -2.8 b 


-7.0 C , -8.5 b 



a Ref. [257]. Relativistic Hartree-Fock calculation. 

b Ref. [285]. Average of two ab initio many-body calculations; core polarization and correlation 
corrections included. 



c Ref. [274]. Estimate found by scaling (with Z) calculations for lighter analogous atoms. Radon 
result scaled from xenon calculation [257]; radium result scaled from mercury calculation [280,279]. 
d Ref. [280,279]. Estimated from the calculation of d a tom(C T ) for mercury performed in Ref. [256]. 
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TABLE XIX. Best limits on P, T- violating parameters at the nucleon level. Signs of the central 
points are omitted. Errors are experimental. Some relevant theoretical works are presented in the 
last column. 



P, T-violating term 


Value 


System 


Exp. 


Theory 


HADRONIC 










neutron EDM d n 


(17±8±6) x 10~ 26 e cm 

(1.9 ±5.4) x 10~ 26 e cm 

(2.6 ±4.0 ±1.6) x 10~ 26 e cm 


199 Hg 

n 
n 


[10] 

[240] 

[241] 


[250,13] 


proton EDM d p 


(1.7 ±0.8 ±0.6) x 10" 24 e cm 
(17 ± 28) x 10~ 24 e cm 


199 Hg 

TIF 


[10] 
[238] 


[250,13,285] 
[215,247] 




i] np = (2.7 ± 1.3 ± 1.0) x 10~ 4 


199 Hg 


[10] 


[280,279] 


9-kNN 


g wNN = (3.0 ± 1.4 ± 1.1) x 10~ 12 


199 Hg 


[10] 


[288,13] 


SEMI-LEPTONIC 










^C SP NNei^ 5 e 


(0.40C^ P + 0.60C^ P ) = (18 ± 8 ± 7) x lO" 8 
(0.40C^ P + 0.60C^ P ) = (6 ± 6) x lO" 8 


199 Rg 

205 rp J 


[10] 
[11] 


[18,13] 
[255] 


^=C PS Ni lb Nee 


C P5 = (1.8 ± 0.8 ± 0.7) x 10~ 6 


199 Rg 


[10] 


[18,13] 




Cl = (5.3 ±2.5 ±2.0) x 10~ 9 


199 Hg 


[10] 


[256] 


LEPTONIC 










electron EDM d e 


(6.9 ± 7.4) x 10~ 28 e cm 

(9.1 ±4.2 ±3.4) x 10~ 27 e cm 


205 rQ 
199 Hg 


[11] 
[10] 


[261] 
[18] 
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TABLE XX. Best limits on P, T-violating parameters at the quark level. Signs of the central 
points are omitted. Errors are experimental. Some relevant theoretical works are presented in the 
last column. 



P, T-violating term 


Value 


System Exp. 


Theory 


HADRONIC 










^k s q\i^qxq 2 q 2 


k s = (3 ± 1 ± 1) x 10~ 7 
k s = (2 ± 5) x 1(T 5 
k s = (3 ± 4 ± 2) x 10~ 5 


199 Hg 

n 
n 


[10] 

[240] 

[241] 


[297,13] 

[297,298,13] 

[297,298,13] 


-^k c s q 1 i-f 5 t a q 1 q 2 t a q 2 

Ql = 92 


k c s = (1 ± 1 ± 1) x 10~ 6 
k c s = (1 ± 2) x 1(T 4 
fc= = (1 ± 1 ± 1) x 10~ 4 


199 Hg 

n 
n 


[10] 

[240] 

[2411 

L J 


[297,13] 

[297,298,13] 

[297,298,13] 


Qi + Q2 


k c s = (5 ± 2 ± 2) x 10~ 4 
^ = (1 ± 2) x 1(T 4 
fcf = (1 ± 1 ± 1) x 10~ 4 


199 Hg 

n 
n 


[10] 

[240] 

[241] 


[297,13] 

[297,298,13] 

[297,298,13] 


C 1 h — 1 


kt = (1 ± 1 ± 1) x 10~ 5 
fct = (1 ± 2) x 10" 5 
kt = (1 ± 1 ± 1) x 10~ 5 


199 Hg 

n 
n 


[10] 

[240] 

[241] 


[297,13] 

[297,298,13] 

[297,298,13] 


^^^i3ui(T l iut a uda al3 t a d 


fc t c = (5 ± 3 ± 2) x 10~ 5 
fcf = (1 ± 2) x 1(T 5 
k c t = (1 ± 1 ± 1) x 10~ 5 


199 Hg 

n 
n 


[10] 

[240] 

[241] 


[297,13] 

[297,298,13] 

[297,298,13] 


CEDMs d and 
EDMs d of quarks 


e(J d - d u ) = (1.5 ± 0.7 ± 0.6) x 10- 26 ecm 

e(d d + 0.5(J U ) + 1.3dd - 0.3d u 
= (3.5 ± 9.8) x 10~ 26 ecm 
= (4.7 ± 7.3 ± 2.9) x 10~ 26 ecm 


199 Hg 

n 
n 


[10] 

[240] 
[241] 


[299] 

[300] 
[300] 


QCD phase 9 


6 = (1.1 ±0.5 ±0.4) x 10- 10 

6 = (1.6 ±4.5) x 10~ 10 

6 = (2.2 ±3.3 ±1.3) x 10~ 10 


199 Hg 

n 
n 


[10] 

[240] 

[241] 


[288,13] 

[301] 

[301] 


SEMI-LEPTONIC 










^k lq qqeij 5 e 


km + k ld = (3 ± 1 ± 1) x 10- 8 
km + k ld = (1 ± 1) x 10~ 8 


199 Hg 

205 rpl 


[10] 
[11] 


[13] 
[13] 


^qi^qee 


k 3q = (2 ± 1 ± 1) x 10~ 8 


199 Hg 


[10] 


[13] 


■^qijscr^qea^e 


k 2q = (5 ± 3 ± 2) x 10~ 9 


199 Hg 


[10] 


[13] 
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FIGURES 





(a) (b) 
FIG. 1. Hartree-Fock (a) direct and exchange (b) diagrams for energies. The solid and dashed 
lines are the electron and Coulomb lines, respectively. 






(a) (b) (c) (d) 

FIG. 2. Second-order correlation diagrams for the valence electron (X operator). Dashed line 
is the Coulomb interaction between core and valence electrons. Loop is the polarization of the 
atomic core which corresponds to the virtual creation of the excited electron and a hole in the core 
shells. Here a is the state of the external electron; n, m are core states; and /3, 7 are states outside 
the core. 



'ay 



FIG. 3. Lowest order screening corrections to the diagram in Fig. 2(a). 
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FIG. 4. Correlation corrections to energy in the Feynman diagram technique. 



co+e 



FIG. 5. Polarization operator. 
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FIG. 6. Screening diagram chain for effective polarization operator. 
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FIG. 7. Rotation of the integration contour over uj. The points indicate the positions of Green 
function poles; the crosses denote the positions of the poles of the polarization operator. 



FIG. 8. Insertion of the hole-particle interaction into the second order correlation correction. 
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FIG. 9. Hole-particle interaction in the polarization operator. 



FIG. 10. Renormalization of the Coulomb line due to the hole-particle interaction and screening. 
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FIG. 11. The electron self-energy operator with screening and hole-particle interaction included. 
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FIG. 12. Chaining of the self-energy operator. 



FIG. 13. Third-order diagrams of the interaction of a hole and particle from the loop with an 
external electron. 





FIG. 14. Correlation corrections to occupied orbitals of closed shells. 





FIG. 15. External field (denoted by a cross) in the external electron lines. 



99 



O & 'CP 

FIG. 16. Structural radiation in diagram 2(a) of the self-energy operator. 




FIG. 17. Lowest order many-body perturbation corrections to the exchange potential; the cross 
denotes the external field. 





FIG. 18. Brueckner-type correlation corrections to the PNC El transition amplitude in first 
order in the weak interaction; the crosses denote the weak interaction and the dashed lines denote 
the electromagnetic interaction. 
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5E 

FIG. 19. External field inside the correlation potential. In diagrams (a) the weak interaction 
is inside the correlation potential (<5£ denotes the change in E due to the weak interaction); this is 
known as the weak correlation potential. Diagrams (b,c) represent the structural radiation (photon 
field inside the correlation potential). In diagram (b) the weak interaction occurs in the external 
lines; in diagram (c) the weak interaction is included in the electromagnetic vertex. 
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FIG. 20. Diagram showing the anapole moment, a, the toroidal current that produces it, j, and 
the magnetic field that the current creates, H. 



d M m 



FIG. 21. Diagram showing the spin helix that occurs due to the parity violating nucleon-nucleus 
interaction. The degree of spin rotation is proportional to the distance from the origin and the 
strength of the weak interaction. 
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FIG. 22. Flow diagram of CP-violation mechanisms at different levels that induce neutron, 
atomic, and molecular EDMs. 




FIG. 23. Constant electric field E inside the nucleus produced by the P, T-odd interaction 
(Schiff moment field). E is directed along the nuclear spin I. 
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